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ABSTRACT 


The present on-going study Is concerned with experimental and theoreti- 
cal research on 2-0 axl symmetric geometries under low speed* nonreacting* 
turbulent, swirling flow conditions. The flow enters the test section and 
proceeds Into a larger chamber [the expansion ratio D/d - 2] via a sudden 
or gradual expansion. Objectives Include recirculation zone characteriza- 
tion, turbulence simulation In swirling recirculating flow* gradual expan- 
sion flowflelds* and further complexities and parameter Influences. 

The present document details progress made during the second half of 
the first year of the study. During this time period, six flowfleld con- 
figurations have been Investigated with sidewall angles a » 90 and 45 deg. 
and swirl vane angles ^ » 0, 45 and 70 deg. Photography of neutral ly- 
buoyant helium-filled soap bubbles, tufts, and Injected smoke helps to 
characterize the time-mean streamlines* recirculation zones and regions of 
highly turbulent flow. Five-hole pitot probe pressure measurements allow 
the determination of time-mean velocities u, v and w. An advanced computer 
code equipped with standard two-equation k-e turbulence model Is used to 
predict corresponding flow situations and to compare results with the experi- 
mental data. 
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1. INTRODUCTION 
1.1 The Problem 

Both experimental and theoretical studies assist In the design and 
development of gas turbine combustion chambers. ^ Up to now designers rely 
heavily on experimental evidence to produce empirical formulae. However* 
traditional design methods are now being supplemented by analytical methods 
(numerical solution of the appropriate governing partial differential equa- 
tions). Computer modeling of combustion processes Is now an established 
fact* but Improvements and new developments (both experimental and theoreti- 
cal) can and should be made* theoretical modeling being aided by carefully 
chosen experiments. 

The accuracy of currently available prediction codes for swirling re- 
circulating confined flows Is In doubt because of questionable turbulence 
models and lack of an experimental data base. A prerequisite to the pre- 
diction of more complex turbulent reacting flows Is the development of 
suitable turbulence models and computer programs for flow prediction under 
nonreacting conditions, with which the present study Is concerned. 

Progress needs to be made on computational methods. Typically* a com- 
puter code of the TEACH-T type Is appropriately modified to Include swirl 
and geometric variations and the two-equation k-e turbulence model simulates 
the mixing characteristics. One problem Is that a turbulence model whose 
basis and parameters are adequate for simple flow situations Is not adequate 
for the more complicated swirling recirculating flow situation. There Is a 
need to develop logical extensions of currently available turbulence models 
to swirling recirculating flows. In the form of modifying their parameters. 
Inclusion of nonIsotropIc effects, and/or more complex modifications. 


1.2 Objectives 


The present research encompasses steady turbulent flow In 2-f) axisym- 
metric geometries* under low speed and nonreacting conditions. The parti> 
cular problem being addressed is concerned with turbulent flow of a given 
turbulence distribution in a round pipe entering an expansion into another 
round pipe, as illustrated in Fig. 1 of Appendix A. The in-coming flow 
may possess a swirl component of velocity via passage through swirl vanes 
at an angle and the sidewall may slope at an angle a, to the main flow 
direction. The resulting main flowfield domain may possess a central 
toroidal recirculation zone CTRZ in the middle of the region on the axis, 
in addition to the possibility of a corner recirculation zone CRZ near the 
upper corner provoked by the rather sudden enlargement of the cross-sectional 
area. Of vital concern is the characterization of flows of this type in 
terms of the effects of side-wall angle a, degree of swirl 4», turbulence 
intensity k^^ of the inlet stream and expansion ratio D/d on the resulting 
flowfield in terms of its time-mean and turbulence quantities. The general 
goal of the on-going research is to perform experiments and complementary 
computations with the idea of doing t.he necessary type of research that will 
yield improved calculation capability. This involves performing experiments 
where time-mean and turbulence quantities are measured, and taking input 
conditions and running an existing prediction code for a variety of test cases 
so as to compare predictions against experiment. Hence the validity of tur- 
bulence model modification can be assessed. In fact they will also be 
deduced directly from the measured stresses and velocity gradients. The 
goals of the on-going research include: 

1. Measurements of mean flow patterns and corner and central 
toroidal recirculation zones. 
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2. Flowfleld mapping of time-moan velocity, normal and shear Reynolds 
stresses, and hence turbulent viscosities In each (1,J) -direction. 

3. Development of a computer program based on the Imperial College 
TEACH-T program for two-dimensional ax 1 symmetric, swirling, con- 
fined jet flows with recirculation regions. 

4. Advances In turbulence models for general swirling recirculating 
flows. Including nonIsotropIc simulation. 

In the experimental portion of the research, a logical sequence of 
experiments Is proposed to establish the effects on the resulting flowfleld 
of: 

1. Swirl strength 

2. Expansion angle 

3. Downstream blockage 

4. Expansion ratio 

5. Upstream turbulence level. 

Complementary computations will be performed corresponding to all the boundary 
conditions of the experiments so as to provide a thorough evaluation of state- 
of-the-art predictive capability with advanced turbulence models. 

Of course, advances attained on the proposed research may be readily 
and usefully assimilated later Into the more complex flowfleld simulation 
schemes. A more complete survey of the present research program, and Its 
significance. Is given In the original proposal, see Ref. 6. 
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The present research project started under NASA Sponsorship on July 
1980 via NASA Grant NAG 3-74 [see Ref. 6]. A discussion of accomplishments 
during the first six months [«'uly 1, 1980 to Dec. 31, 1980] was given In 
the First Semi-Annual Status Report.^ There, emphasis was on Task 1 [Basics] 
and Task 2 [Time-Mean Flow Characterization] of the original schedule of the 
Proposal.^ Discussion was Included on: 

1. Test Facility - the design and construction of three expansion 
blocks [o » 90, 70 and 45 deg.] and of the confined jet facility 
shown schematically In Fig. 4 of Appendix A, Including the design 
and construction of a variable angle swirler. 

2. Experimental Techniques - the ;.et up of flow visualization cap- 
ability with neutral ly-buoyant helium-filled soap bubbles, smoke 
tracers and tufts. Including slit lighting and photographic 
techniques. 

3. Computational Code Development - the formulation of an advanced 
version* of the Imperial College TEACH-T computer program,* In- 
cluding swirl flwW and a stairstep approximation to the sloping 
side wall. 

4. Flowfield Characterization - with emphasis on the flow character 

and predictive capability, as evidenced by streamline patterns 

and recirculation zones from the bubble photography technique. 

Results presented at that time [Dec. 30, 1980] were of a preliminary nature \ 

, • 1 

and were Incorported In a draft paper. The final form of that paper^ Is , 

I 

Included herein as Appendix A. Complete details about the experimental 
facility and techniques may be found in Section 4 of Appendix D, which con- 
sists of a recently-completed PhD Thesis.* In sumnary form, some of the 

' 

details also appear in a recent AIAA paper,* Included herein as Appendix C. j 

4 ^ * 
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1.4 The Present Contribution 

j 

The present report Is concerned with results and publications ensuing 
during the second half of the first year of the project [Jan. 1, 1981 to 
July 31, 1981]. Recent progress Is discussed In Section 2 where more 
realistic flowfleld details, from the extensive flow visualization and five- 
hole pitot probe surveys, are presented. Predictions from the computer code 
are Included and compared with the data when appropriate. The Closure 
[Section 3] summarizes the achievements, clarifies the status of current 
activity and Includes a list of publications. Copies of four recent pub- 
11cat1ons‘~‘ are appended, while the user's guide to the finally-developed 
computer code' Is still being processed by NASA. 



2. RECENT PROGRESS 


The original proposal calls for progress on Task 3 [Turbulent Fluctua- 
tion Characterization] and Task 4 [Time-Mean Flow Characterization for the 
strongly swirling flow case with swirl vane angle ^ = 70 deg.]. Because of 
some delay In Implementation of the hot-wire techniques, discussion of that 
progress Is being reserved for the next Semi-Annual Report. The emphasis 
now Is on flow visualization, pitot probe time-mean velocity measurements 
and flowfleld predictions In six distinct flowfleld configurations. These 
have side-wall angles a ■ 90 and 45 deg. and swirl vane angles (> > 0, 45 
and 70 deg. Copies of four recent publications^*’ are supplied as Appendices 
A through 0, to which reference and further details about developments may 
be found. The discussion concentrates on the more thorough documentation of 
Appendix D [a recently-completed PhD Thesis'*], but It may be noted that more 
brief details may be found In the research papers of Appendices A through C. 

2.1 Flow Visualization 

Section 4.3 of Appendix 0 describes the flow visualization techniques 
used In the study. The discussion centers on: 

1. Illumination anu Photography, 

2. Neutral 1y-6uoyant Helium-Filled Soap Bubbles, 

3. Smoke-Wire, 

4. Tufts. 

Recirculation zones are Important to combustor designers because most 
of the burning occurs In these regions and they exhibit the highest tempera- 
tures. The size and location of these regions In the present Isothermal 
flows are readily deduced from flow visualization photographs of tufts, 
smoke, and bubbles responding to the experimental flowfleld patterns. Re- 
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suiting dividing streamline sketches as well as selected photographs of the 
visualization experiments are presented and discussed In Section 5 of 
Appendix 0. Main results Include: 

1. Artistic Impressions of Streamline Patterns 

Photographs of each of the six flowflelds resulting from ^ « 0, 

45, and 70 degrees with o > 90 and 45 degrees have been examined 
In detail for each of the three flow visualization methods employed. 
The characteristics of the overall flowfleld are Illustrated and 
discussed via the resulting t1me>mean dividing streamline patterns. 
These are sketched In Figs. 20 and 21 [of Section 5.1 of Appendix 0] 
from Information obtained from the entire collection of flow 
visualization photographs. Results from the smoke-wire experiment 
are utilized near the Inlet, whereas tuft and bubble data are used 
In approximating the size and shape of the recirculation zones 
downstream. Also, bubble flow patterns reveal the existence of a 
processing vortex core, which occurs downstream of the central 
region. In the swirl flow cases Investigated. This, discussed at 
length by Syred and Beer,^° Is observed near the centerline extend- 
ing from the end of the central region to the test section exit. 

The axial location, where the vortex core begins fluctuates, 
ranges approximately from x/D = 1.25 to 1.75. This vortex core Is 
essentially a three-dimensional time-dependent phenomenon, which 
occurs as a swirling region of negligible axial velocity whose 
center winds around the test section centerline. 

2. Flow Visualization 

Tuft visualization Is very Important In that It supplies an overall 
view of local flow direction. Photographs at various shutter 
speeds were obtained. Slower speeds show more of the temporal 
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behavior, although the tufts are sometimes not distinctly visible 
In portions of the flowfleld. Some of the more noteworthy photo- 
graphs are presented here, encompassing a range of shutter speeds. 
Velocities In recirculation zones are often somewhat lower than 
In other portions of the flowfleld, and thus under such conditions 
there may be Insufficient drag on a tuft to align It accurately 
with the local flow direction. However, this Is taken Into con- 
sideration In Interpreting the photographed results. Figures 22 
through 25 Illustrate the effects which are described In Section 

5.2.1 of Appendix D. 

Local details In the nonswirling flowflelds are clearly revealed 
through the visualization of streakllnes Indicated from the 
generation of Illuminated smoke. In the swirling flow cases, 
strong mixing diffuses the smoke so that streakllnes are not 
distinguishable. However, under such conditions recirculation 
zone outlines are visible, especially In the region near the smoke- 
generation wire. A selected photograph Is exhibited and discussed 
for each of the three abrupt expansion flowflelds, see Section 

5.2.2 of Appendix 0, with Figs. 26 and 2/. 

Soap bubbles Injected Into the flow upstream of the test section 
trace pathllnes clearly when Illuminated. In relatively lower 
turbulence Intensity portions of the flowfleld mean flow directions 
can be obtained by ensemble averaging local tangents to pathllnes 
traced out by soap bubbles. This helps define the flowfleld 
geometry In terms of the outline of recirculation regions. Figures 
28 and 29 of Appendix 0 show the flow patterns tdilch we described 
In Section 5.2.3 of Appendix 0. 


3. Parametric Effects 


Streamline plots allow recirculation zones to be characterizr^ 
parametrically. The parametric effects deduced from the experi- 
mental study are illustrated in Fig. 33 of Appendix 0* which 
show the effects of a and ^ on the corner and central recircula- 
tion zone lengths. Several observations should be noted. Firstly* 
zone lengths are only slightly affected by a. Secondly* as side- 
wall angle a decreases from 90 to 45 degrees* the zone lengths 
typically decrease slightly and the inlet flow is encouraged to 
impinge on the confining walls. Thirdly* the corner recirculation 
length decreases upon increasing the swirl vane angle ^ from 0 to 
45 degrees* a parameter change which also provokes the existence of 
a central recirculation bubble and processing vortex core. In- 
creasing the vane angle to 70 degrees only slightly enlarges the 
central zone and the vortex core. These observations are discussed 
further in Section 5.4 of Appendix 0. 


2.2 F1vt~Hole Pitot Probe Measurements 


One of the simplest Instruments capable of simultaneously sensing 
both magnitude and direction of the local velocity vector Is the five- 
hole pitct probe. The particular probe employed In this study Is 
model 0C-125-12-C0 from United Sensor and Control Corp, It has a 
3.2 mm diameter sensing tip and shaft containing five tubes. The 
sensing head Is hook-shaped to allow probe shaft rotation without 
altering the probe tip location. The Instrumentation system. In 
addition to the five-hole pitot probe, consists of a manual traverse 
mechanism, two five-way ball valves, a very sensitive pressure 
transducer, a power supply, and an Integrating voltmeter. The 
differential pressure transducer Is model 5900 from Oatametrics, Inc. 

The output Is read as the d.c. signal from a TSI model 1076 Integrating 
voltmeter. The relevant part of Appendix 0 Is Section 4.4 (with 
Figures 13 through 19) which Includes Information on: 

1. Instrumentation, 

2. Measurement Technique, 

3. Data Reduction, 

4. Calibration. 

Section 6.1 of Appendix 0 presents and discusses the five-hole 
pitot probe measurements that have been completed on both nonswirling and 
swirling confined Jet flows In the test facility. It Is anticipated 
that these measurements will provide the most accurate data on time- 
mean velocities, and their spatial gradients. The results are being 
used In evaluation of the predictive capability of the computer simula- 
tion and In future evaluation and possible modification of the turbulence 
model. 
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Figures 34, 35, and 36 of Appendix D show the axial and swirl 
time-mean velocity profiles, for swirl vane angles (^ « 0, 45 and 70 
deg. with a sudden expansion [side-wall angle a > 90 deg.]. Measure 
ments are obtained with the five-hole pitot probe In the manner just 
discussed. Note that In these and subsequent figures, different 
scales are used for the normalized axial and swirl velocities. 

Radial velocities are consistently much smaller than the axial and 
swirl components, and they need not be presented here. Figures 
37, 38, and 39 of Appendix 0 exhibit velocities for the same 
sequence of flowflelds with side-wall angle a « 45 degrees. The 
Inlet profiles were not measured In this geometry because of the 
presente of the expansion block Interferes with probe positioning. 
Effects of swirl vane angle <t> on velocities, similar to those 
found for the sudden expansion cases, are found in the flowfleld 
sequence for this test section geometry. The major difference Is 
that the sloping wall encourages the Inlet flow to accelerate 
near the top wall. Also, It tends to shorten or obliterate the 
corner recirculation region. 

In the absence of pertinent measurements, previous predictions 
of gas turbine combustor flowflelds anployed an approximation for 
Inlet velocity boundary conditions. It has generally been assumed 
that V = 0 and that both u and w exhibit flat velocity profiles. 
However, the present measurements Indicate this is an unrealistic 
estimate, with sharply peaked u and w profiles. The difference 
results from the use of a swirler with ten flat blades with pitch/ 
chord ratio of unity, imperfect blade effic1enc:y, the existence 
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of a hub, and the fact that the downstream edge of the swirl vanes of the 1 

test facility Is actually located approximately 4 cm upstream of the 1 

flow expansion corner, where x/0 « 0. This swirler location allows j 

I 

the central recirculation zone to begin upstream of x/0 » o, thereby | 

changing the velocity profiles there. i 

2.3 Flowfleld Predictions 1 

A primitive pressure-velocity variable finite difference computer | 

\ 

code has been developed to predict swirling recirculating Inert turbu- 
lent flows In axisymmetric combustors In general, and for application 
to the present specific Idealized combustion chamber with sudden or gradual 
expansion. The technique Involves a staggered grid system for axial and 
radial velocities, a line relaxation procedure for efficient solution 

! 

of the equations, a two-equation k-e turbulence model, a stairstep i 

boundary representation of the expansion flow, and realistic accommodation 
of swirl effects. The development Is based on the 1974 Imperial College 
TEACH-T computer code.^ The finally developed computer program (written 
In Fortran 4) Is code-named STARPIC (mnemonic for swirling turbulent 
axisymmetric recirculating flows In £ract1cal Isothermal combustor 
geometries). A complete report* has been amalgamated and delivered to 
the NASA technical monitor, together with a copy of the card deck Illustrating 
a sample computation of air flow through a 45 deg. expansion (a » 45 deg.) 
from an Inlet pipe to a larger pipe (0/d » 2). The delivered code and 
description provide computer runs through a range of seven Inlet swirl 
vane angles ( ■ arctan ^Rual to 0, 45, 55, 60, 65, 68 and 

70 deg. The extensive document presents details of the computational 
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solution procedure. It serves as a user's manual and deals with the 
computational problemi showing how the mathematical basis and computational 
scheme may be translated Into a computer program. A flow chart, Fortran 4 
listing, notes about various subroutines and a user's guide are supplied 
as an aid to prospective users of the code. Only a brief description of 
the technique need be Included, see Section 3 of Appendix 0. 

Application of this STARPIC computer code has been made during the 
first year to the simulation of Isothermal airflow In the axisymmetric 
test facility with the following conditions: diameter expansion ratio 
0/d Is 2.0, Inlet Reynolds number Re^ * 1.26 x 10®, wall expansion angle 
o * 90, 70, and 45 deg. and swirl vane angle varies from 0 to 70 deg. 

All results are obtained via a nonuniform grid system because of the 
resulting enhancement of solution accuracy. The Inlet profiles of axial 
velocity u and swirl velocity w are Idealized as "flat" (that Is constant- 
valued). Section 6.2 of Appendix 0 describes predicted u and w velocity 
profiles for the six flowfield configurations, with the details In Figs. 

40 through 45. Corresponding streamline predictions appear In Figs. 30 
and 31, and deduced parametric effects on recirculation zone lengths In 
Fig. 32 (see Appendix 0). The general trends are similar to the experimental 
findings. 

Note that. In these predictions, the Inlet flow angle Is taken to be 
the swirl vane angle (|>, whereas In reality the blades are not 100% efficient. 
The experimental Inlet data were not available at the time the predictions 
were made, and the flat Inlet profile assumption Is now known to be In- 
adequate for precise simulation of the present study. In which the Inlet 
flow conditions are achieved via a ten-blade swirler with pitch/chord 
ratio of unity. As documented In the experimental results, the present 


Inlet profiles are highly nonunlform. Because of this, a test case was 
Investigated In which experimental Inlet profiles of u and w were taken 
as Inlet boundary conditions, with zero radial velocity also specified. 

It was for the case of swirl vane angle » 45 deg. with sudden expansion 
a » 90 deg., and the details appear in Section 6.3 of Appendix 0, with 
Fig. 46 comparing predicted and measured profiles. The discrepancy near 
the Inlet may result partially from backflow at the inlet boundary [which 
could not be measured because of transient effects], the lack of probe 
sensitivity to small velocities near the centerline, and/or the use of 
the standard two-equation k-e turbulence model. In this connection. It 
may be noted that a turbulence model whose basis and parameters are 
adequate for simple flow situations Is not adequate to handle the more 
complicated swirling recirculating flow situation. Nevertheless, good 
predictions are available for a less difficult test case [co-swirl and 
counter-swirl flow In a pipe using the data of Kef. 11], see Section 
6.5 of Appendix 0. This turbulence simulation problem In complex flow- 
fields Is clearly an area of current research Interest. Currently, pre- 
dictions are being made for a variety of more realistic Inlet station 
velocity and turbulence specifications, deduced from measurements being 
made very close to the upstream swirler. The encouraging results will 
be available very shortly. 
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3. 


CLOSURE 


3.1 Summary 

The main objectives of the research program are to determine the 
effects of swirl and combustor geometry on Isothermal flowfleld patterns, 
time-mean velocities and turbulence quantities, and to establish an Im- 
proved simulation In the form of a computer prediction code equipped with 
a suitable turbulence model. The first year's effort has focused on the 
identification of recirculation regions from flow visualization, mean 
velocity measurements with a five-hole pitot probe, and the comparison of 
experimental results with flowfleld predictions using an advanced computer 
code equipped with standard two-equation turbulence model. Discrepancy 
between experiment and calculation establishes the need for Improvements 
In the modeling of swirling confined jets, especially as applied to 
combustor geometry flows. 

3.2 Current and Future Work 

The experiments being performed are designed to provide the Informa- 
tion necessary for turbulence modeling development which is a major part cf 
the goals of the on-going research. Single normal hot-wire measurements 
In nonswirling flow have established that the experimental techniques are 
producing reliable results In comparison with past researchers (or. nonswirling 
flows) such as Chaturvedi Six-orientation hot-wire measuremc.nts In the 
facility are In progress from which estimates of the time-mean velocity 
components and the Reynolds stresses can be made.'*’*' A major portion of 
these measurements are now complete for the nonswirling flowfleld. They are 
Incorporated In an MS Thesis which Is currently being finalized.'* 

Although the single hot-wire (six-oricntatlon) measurements are 
convenient and provide a great deal of Information on the time-mean velocity 


and on the kinettc energy of turbulence, they have a shortcoming In that 
measurements of the Reynolds shear stresses are less accurate than probably 
will be necessary to meet the goals of our research.^’ Consequently, 
crossed hot-wire capability Is also being developed to measure turbulence 
properties, and specifically the turbulent shear stresses. This technique 
Is now being used for the nonswirling flow and another student's report on 
this effort is nearly complete. This work is currently being extended to a 
three-wire hot-wire technique for swirling flows. These measurements are 
most accurately accomplished by simultaneously digitizing the analog 
signals from the three hot-wires of the probe and reducing the data on a 
digital computer providing adequate storage Is available In the computer. 

This Is the technique being followed presently. 

The second year's work extends measurements of turbulence quantities 
to confined jet flowflelds with swirl and with a realistic downstream 
blockage. Later work will employ different expansion ratios and inlet 
turbulence levels. In addition. Computational studies are now in progress 
which Illuminate further on the need to specify Inlet values very accurately, 
this clearly being a pre-requisite to turbulence model development. Measurements 
of turbulence and mean velocity components will provide the necessary experi- 
mental data with which to evaluate the standard k-c turbulence model and to 
recommend advances In this model specifically to more realistically accommodate 
swirling flows. 

3.3 Publications 

The following list of publications covers research conducted with the 
support of this research grant: 

1. Rhode, D. L., Lllley, 0. G., and McLaughlin, 0. K., "On the 

Prediction of Swirling Flowflelds Found In Axisymmetric Combustor 
Geometries", ASME Symposlun on fluid Mechanics of Combustion 
Systems, Boulder, CO, June 2?-24, 1981, pp. 257-266. . 
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2. L111ey« D. G., Rhode, D. L., and Samples, J. U., Prediction of 
Swirling Reacting Flow in Ramjet Combustors, AIAA Paper No. 

81- 1485, Colorado Springs, CO, July 27-29, 1981. 

3. Rhode, 0. L., “Predictions and Measurements of Isothermal Flowfields 
in Axi symmetric Combustor Geometries, Ph.D. Thesis, Oklahoma State 
University, Stillwater, OK, Dec., 1981. 

4. Lilley, D. G. , and Rhode, 0. L., STARPIC: A Computer Code for 
Swirling Turbulent Axi symmetric Recirculation Flows in Practical 
Isothermal Combustor Geometries, NASA CR-3442, 1981. 

5. Rhode, D. L., Lilley, 0. G., and McLaughlin, D. K. , Mean Flowfields 
in Axi symmetric Combustor Geometries with Swirl, AIAA Paper No. 

82- 0177, Orlando, Florida, Jan. 11-14, 1982. 
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ABSTRACT C TurhutonC V 4 lu« 


Combuncor Ok^dcllns hat rtjchot! th« «Cag** whore Che 
ao«C useful research accivicies are likely co be on 
specific sub-problems of che general three- Jtoenslonel 
turbulent reacting flow problem. The present study is 
concerned with a tlirely fluid dynamic research task of 
interest to the combustor aodeling cofrrunity. Nuiretical 
computations have been undertaken for a basic tvo-dlmen* 
sional axisyiBReCrlc flowfleld which is similar to that 
found in a conventional gas turbine combustor. A swirl- 
ing nonrearcing flow enters a latter chamber via a sud- 
den or gradual expansion. The calculation tsethod in- 
cludes a stairstep boundary representation of the ex- 
pansion flow, a conventional k-C turbulence model and 
realistic accomr.odat ion of swirl effects. The results 
include recirculation zone characterization and oredict- 
ed mean streamline patterns. In addition, an e<^ erimental 
evaluation using flow visualization of neutrcl '> -buoyant 
hel ium- f i I led soap bubbles is yielding ver/ pr lising 
results. Successful outcomes of the work can be incor- 
porated into the mure combustion- and hardware-oriented 
activities of gas turbine engine manufacturers , includ- 
ing incorporating Che modeling aspects into already 
existing comprehensive numerical solution procedures. 

NOMlNCl.ArURE 

a Coupling coefficient 

c Constant 

0 Chamber dl are ter 

d Nozzle diameter 

G Axial flux of mocentua. k-generatlon term 

l,J Nesh point 

k Kinetic energy of turbulence 

P Time-*r*fin pressure 

R Residual source 

S Swirl number • 2C^/(G|^d), source term 

(with subscript) 

Sp, Components of linearized source term 

V • (u,v,w) Time-mean velocity (in x, r, 0 direction) 
x,r.0 

y 
r 

c 

u 

d 
3 

r 

1 


Subscri pts 

in 

i 

m 

N.a.e.w 

0 

P.N.S.C.W 


Axial, radial, azimuthal cylindrical polar 
co-ordl lates 

Distance normal to a wall 
Turbulent exchange coefficient 
Turbulence energy dissipation rate 
Effective viscosity 
Time-mean density 
Frandt l-Schmidt nur.ber 
Wall shear stress 

Swirl vane angle ( tan* ‘ fw^j^/u^jj) | , general 
dependent variable 


Inlet conditions 

Laminar value 

Station maximum value 

North, south, east, west faces of cell 

Value at orifice of jet or inlet to flowfleld 

Point, north, south, east, west neighbors 


Supe rscripts 


Prelisiinary u,v and p field based on esti- 
skated pressure field p* 

Correction value to u*, v*, p* to get u,v,p 


INTROOI.'CTION 

The combustor of the g*s turbine engine contains 
high intensity combustion and, as far as possible, must 
bum fuel completely, caunc little pressure drop, pro- 
duce gases of nearly uniform* temperature, occupy sm.ill 
volume, and maintain stable combustion over a wide rang** 
of operating conditioned)- In design situations, the 
engineer has to seek an optimum path between irrecoc.cil- 
able altemativen of, for example, efTiciency and pollu- 
tion. The general aim of most rese.irch invest ig«t ions 
is to provide information which is useful to designers 
by 'characterizing* or 'modeling* certain features of 
the phenomenon in questlori?). Up to now designers rely 
heavily on experimental evidence to produce empirical 
formulae. However, traditional design methods are now 
being supplemented by analytical methods (numerical 
solution of the appropriate governing partial differen- 
tial equations). Computer modeling of combustion pro- 
cesses is now an established fact, bur improvements and 
new developfrents (both experimental and theoretical) can 
and should be made, theoretical as^deling being aided by 
carefully chosen expe rimentsC 1) . 

Ihe present paper addresses research that is re- 
stricted to steady turbulent flow in axisytaretric geo- 
metries, under low speed and nonreacting conditions - 
a study area highlighted at a recent workshop (1) as a 
fundamental research requirement in cembustion siodeling. 
The particular problem is concerned with turbulent flow 
of a given turbulence distribution in a round pipe en- 
tering an expansion into another round pipe, as illus- 
trated in Fig. 1. The in-coming flow may possess a 
swirl component of velocity via passage throuxh swirl 
vanes at angle 0 [equal approxlmafely to tan*^(w^n/u^;i)| , 
and the side-wall may slope at an angle i, to the suin 
flow direction. The resulting flowfleld doenain may 
possess a central toroidal rec irculat ion zone CTR2 in 
the sUddle of the region on che axis, in addition to the 
possibility of a comer recirculation zone CRZ near the 
upper comer provoked by the rather sudden enlargement 
of the cross-sect ional area. Of vital concern is the 
characterlzat ion of flows of this type in terms of the 
e^f facts of side -wall angle 3, degree of swirl 0, turbu- 
lence intensity k^u of the inlet stream and expansion 
ratio 0/d on the resulting flowfleld in terms of its 
time-mean and turbulence quantities. Such problesu have 
received little attention, yet there is a definite need 
for work in this area(4-6). 

Prediction work in the form of computing the flow- 
field in the conf iguration Just described is being com- 
plemented by an associated experimental study. Consi- 
deration is given to recent work in the finite differetk>* 
solution, via a primitive variable code, of axisvmmetric 
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•vlrl flow in th* combuiitor g*oMti 7 of Fig. 1, whort 
th9 Inltt txpantlnn tidcwall imv tlopt obliqutly to th« 
cenlral axis. Thua a ayatamatic paranatrlc Invaatlga- 
tlon «ay ba contaaplatad on tha affact of aldavall angla 
a and dagraa of awirl d cm tha raaulting flowfiald pro- 
ducad. Baalc guvamlng aquatlona ara praaantad, togath* 
ar with a brlaf daacrlptlon of tha alanilatlon and aolu* 
lion tachnlqua. Tha aquationa ara alllptic in charactar 
and ara aolvad via an advancad varalon of tha Inparlal 
Collaga TtACH-T coaiputar prograa (7). 


SWIRL VANE ancle 

inlet the 



CVS 






ki/2 


la) SCHEMATIC OF TEST SCCTiON 
GEOMETRY 



' (b) EXPECTED RECIRCULATION ZONES 


Fig. 1. Tha flowflald baing invaat Igatad. 

Coapariann with avalliblt axperlmantal turbulant 
flow BiaasuraTnantt aaaiata In conf inning tha final pra- 
dictiva capability. One auch axparioant, undarway at 
Oklanoaa Stata Univaraity, la concamad with maaaurlng 
tha af facta of awirl and aida-wall angla on tha at r<>as^ 
linaa, a^an flow and turbulanca paramatara in non* 
raacting flow. Tha facility and axparinantal datalla ara 
includad. For tha praaant papar a praliminarv avaluatioo 
of tha accuracy of tha computad flowfialda ia accoopliah- 
ad by coaptriaon with flow viaualizat ion uaing nautrally- 
buoyant haHu8»-f illad aoap bubblaa aa tracar partlclaa. 
Fhotographa of tha bubblaa can ba intarpratad t^ yiald 
tima-»aan flowtiald napa which dafina approximataly tha 
boundariaa of racirculation ragicna and ragiona of high- 
ly turbulant flow. Major faaturaa of atrongly awlrling 
flow characta. irat ion ara than praaantad ao aa to axam- 
pllfy tha cu-reni pradictiva capability in taraa of valo- 
city profllaa. atraamlina pattama and racirculation 
zona charactari zat ion. Tlia final cloaura auunaritaa tha 
acliiavananta. 

THEOXETICAL MODEL 

T lowfi a ld Sjaula t i on F rob 1 a« 

Numarical procaduraa for predicting ronbuation 
cl.ambar flowa (which in practica ara complax thrac-diibcn- 
aional turbulant reacting flowa) raat on two four i ona: 
mathematical modala of phyairal proccaaaa (rurbujtnc , 
radiation, combuation and two-phaaa af facta) and coitputar 
programa for aolving tha appropriata Irval of aophiatl- 
cation of tha raaulting dlffarential aquationa. In tha 
modallng and prediction of laotharmal rriabuator flow- 
fialda, tha flow ia aimulatad by aipultanaoua nonlinear 
partial iiffarantial aquationa. Theac may ba parabolic 
(bojndary layer type) but ara more often allipric (re- 
circulating ty'^a) and tha aolutlon achena diffara accord- 
ing to tha ca«agory. Fxparting atrongly curved atrram- 
linaa and recirculation tonaa, tha problrv ia contrmad 
with tha latter. T)<a turbulant Ra/nolde aquationa for 


conaarvation of maaa, aomantuN (with x, r« t velocity 
coorponantt u, v, w), turbultnct kinatic anargy k, and 
turbulanca diaaipation rata C govam thia two-dlaanaion- 
al axiaynaiatric ataady flow. Thay laay ba axpraaaad in 
ganaral foni (7-10) 


“ (Our ^ <tT 


JL 

br 




A (,r ^)j • g 

Sr Sr^^ > 


( 1 ) 


whara i rapraaanta any of tha dwpandant variablaa and 
corraaponding valuat of ara indicatad in Tabla I« 
whara implicit uaa ia mada of tha t«#o-aquation k-C tur- 
bulanca modal, and othar quant it iaa ara dafinad aa 
follow!! 
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Tabla 1. Sourca Tarn in the Canaral Equation 
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C - CqOC 
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(C|CC - cipc*)/k 


T>*raa tranaport aquationa ara all aimilar and con- 
tain tarns for the convection and dlffuaion (via tur- 
bulant flux tarna) and **acurct** Sg of a gaf<aral varlstla 
t (which containa term describing the ganaratiun 
(creation] and consumption (diaaipation] of d). Kotico 
tha usual turbulant dif fuaioo-flux (atreas-rata of 
strain typa) laws. In the first three of *haa# aqua- 
tions (for 4 • u, V and w]Tg stands for the affective 
viac^aity u (although latar nonlaotropic variant# of tha 
turbulanca uodal may not treat all such value# at a 
point to ba equal, but rather to depend also o® tha 
direction aitu orientation of tha appropriate ahaarj# 

In othar aquationa Tg stands lor ii divided by aa 
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ropriat* Prandtl or Schaidt nuaibor, in tbit cat# 0^ 

I 9£t Constanta appearing In tha aquation tat ara 
/an valuaa racoonandad In Raf. 7t Cy • 0.09« 

• l.OO, Cl • 1.44, C 2 • 1.92, • l.OO and9£*1.21. 

Thaaa aqua t Iona hava to ba aoSvad for tiaa aaan 
aaaura p and valoctty coaponantt u, a and v. Than 
bar uaaful daalgnar Inforantion, IlUa atraaallna plota, 
aakavay and raattachaant points, racircularion tonaa 
d atagn^rion polnta, tor axampla, aay ba raadlly da- 
cad. Tba lltaratura (7*10) aay ba raliad to pro- 

da full dataila on poaalbla turbulanca aodata and 
nutation and aolution tacbniquai, and only blghJfghta 

* tha priaittva variable approach naad Ua given bara 

• tba context of tba specific problem being invaatiga- 
id. 


>0 Solution Tech*: • qua 

Tba praaaut prograai ia baaed on concepts from tha 
EACH coda (7) because of, aaong other things, Ita con- 
eptual ainpllcity and ita conputational efficiency, 
ularian finite difference approxisMtions to tha steady 
lliptic governing aquations ara aolvad implicitly for 
ba priaitiva pressure and velocity variables, via tha 
IKPLC [aani-^apllcit natbod for £raosura-Unkad aqua- 
ions] technique using tna luhA i w I-UIagcr.al 
.Igoritbm] In coluasis. Tba Incorporation of several 
ipacial features enhances tbt accuracy and rapidity of 
:onvarganca of tba finally developed computer program* 
Vdvancas to tba original formulation (7) ara described 
aith applications to swlri flow coad>ustor problems alsa- 
ahare (8-13). An earlier simulation and solution taebni- 
)ua using the derived stream function and vortlcity 
variables (14) has also bean extended and applied to 
similar problems (19-17). 

The mesh system of staggered grids Is shown in 
Fig. 2. All variables except u- and v-veloclties are 
stored at g.id point F, while these velocities are 
stored at locations corresponding to the arrows shown 
midway between grid points. Thus three staggered mesh 
systems produce a triad of points encompassed by a boom- 
erang-shaped envelope associated with each grid point 
F. Hence u^j is Che axial velocity associated with 
location (i,J), although its true location is (i-l/2,j). 



Fig. 2. Staggered grid system and nctation. 

Inuodar^es are posit lotied midway between grid points so 
that they coincide with normal velocities as in Fig. 3, 
which shows an example of grid spocif Ication for the 
geometry urJer cons iderr* ion. A recent ax>di flea' Ion (18) 
has been to position the u- aid v-ceil boundaries to be 
either exactly on or exactly halfway between the grid 
intersection points F. This improves the lecuracy of 
calculating pressure gradient effects on ihe u and v 
velocities (as pressure values nrv lie directly on ap- 
propriate cell facts) and of calculating cerraln cell 
boundary fluxes. 
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wall 



Fig. 3. An example of a coarse grid system being 
employed to fit the flow domain. 

The finite difference equations for each 9 are 
obtained by integrating Eq. <i> ov«i lUm •pprepri*'* 
control volume (centered about Che location of I) and 
expressing the result in tenna of neighboring grid 
point values. The convection and diffusion terms be- 
come surface integrals of thi convection and diffusive 
fluxes. For these, a hybrid scheme is used (7), which 
is a combination of the so-called central and upwind 
differencing approaches. The source term is linearised 
in a form which is such as to encuursge convergence of 
the iterative solution procedure. 

loundery Co nd itions 

The flowfield is covered with a nonunlfom rectan- 
gular grid fvstem; typically the boundary of the solu- 
tion domain falls halfway between its imsiedlate nearby 
parallel gridlines; and clearly specification of the x 
and r coordinates of the gridilnes, together with in- 
fonxtition concerned with the piisitioo of the sloping 
sidewall boundary, ia aufficienc to determine the fl<n#- 
ficld uf interest. The aloplng boundary is simulated 
by means of a atalrstep approach, and Interior points 
adjacent to the boundary must 'feel* the boundary in 
the usual way. Usually the standard coupling coeffl- 
clsnt of Intsreac is set to ssro and the wall effect 
given by way of a false source. Velocities or ptsssurcs 
are typically known along the boundsiiss. Implemenra- 
Clon of the former la atraightforvaiU. velocity values 
are known and 3p*/3n • 0 at the boundary. In the latter 
case pressure is sn«cified and one acta p * p* and p'* C 
at the boundarv. For the outflow boundary, velocitlew 
may be adjuatad to enforca global mats continuity. At 
the inflow boundary^ variablaa are given definite fixed 
values associated with known conditions. 

A final observation concerns valocities, k and C 
at wall boundaries. To avoid the need for detailed 
calculations in tha nar.rvall raglons, aquationa called 
*wall functions* ara introduced and used in the finite 
difference e. .culations at naar-wall points. They oc- 
cur In the momenta equaiiont and k-generxtion terms, 
and their Implementation is discussed elsevhcri (19-20), 
together with appropriate near-wall r specif ication. 

The e.fect of swirl on wall function specification 
la handled as follows.' The previous Ideas are extended 
to find total tangential wall thaar stress near bound- 
*.rlea (involving x and 9 directions for an r • constant 
wall, and r and 9 directions for an x • constant wall). 
Then appropriate componan ta are deduced direct Iv (for 
the u and w velocities which are tangential to ar r • 
constant wall, and v and w velocities which ate tongen- 



till to an X * constant wall]. The effects on u, v and 
V momenttitD equations are incorporated via the usual 
I linearized source technique (18), 

Itera tioft Sch eine 

The numerical schene, described at length 
elsewhere (7-10), considers variables at r^oints along a 
vertical line as unknowns (that is, N, P and S for each 
point P), whereas the latest values of each E and V 
neighbor are treated as knovm. The tridlaponal matrix 
algorithm (TDMA) is employed to sequentially solve the 
system of algebraic equations resulting for each verti- 
cal line, sweeping from left to right. Some underrela- 
xation is used and final convergence is attained when 
tlie residual source becomes smaller than some reference 
value for each variable at each point in the flow domain. 
Several techniques are incorporated to enhap/'.e accuracy 
and efficiency, and a complete description of the final- 
ly developed computer program is currently being pre- 
pared (18). 

EXPERIMENTAL APPROACH 

P i ous Re lev ant Experimental Studies 

There have been several experiments performed with 
nonreacting flows in expansion geometries that have been 
reported in the literature (21-30). References 21 and 
22 also include flovfield predictions, made with ver- 
sions of the TEACH-T computer program (7). These ex- 
periments include tir<r-mear. velocity measurements [with 
hot-wire and ,;<tot probes and laser Doppler anemometry), 
turbulence measurements [with hot-wires and laser ane- 
mometers] .»nd flow visualization. The majority of the 
measurements were made in nonswirling flows (24-30), 
however some noteworthy experiments were made in swirl- 
ing confined jets (21,23). Direct comparison between 
the results of the cited experiments and the present 
experimental results is gene»'allv not possible because 
of differences in geometry. However, in the nonswirling 
jet comparisons were possible with experiments of 
Chaturvedi (30), who measured mean and turbulent flow 
quantities downstream of a sudden expansion of diameter 
ratio 2.0 and various expansion sidc'w*all angles a. 
Measurements of near, velocity in regions of high turbu- 
lenrje intensity and where the direction of the velocity 
vector is unknc»wr. were made with a pitot tube. Mean 
velocity was also measured with a constant temperature 
hot-wire anemometer using a single wire. In addition, 
a cross-wire was used to measure all the Reynolds 
St resses. 

T est Faci 1 ity 

The present experiments have been conducted in the 
confined jet facility shown schematically in Fig. 4. 



Fig. 4. Schematic of flow facility. 

The facility has an axial flow fan whose speed can be 
^changrJ by altering the drive pulley combination. Nu- 
merous fine screens and straws produce flow in the set- 
tling chamber of relatively low turbulence intensity, 
rhf c^'ntractlon section leading to the lest section hat 


been designed by the method of horel (31) to produce 
a Tninimum adverse pressure gradient on uhe boundary 
layer and thus avoid unsteady problems associated with 
local separation regions. The sudden expansion centists 
of a 15 CD diameter circular jet nozzle, exiting abrupt- 
ly into a 30 cm diameter test section as sho%m in 
Fig. 1, The suhscanrial size of this test model pro- 
vides excellent probe resolution for hot-wire measure* 
ments which are currently underway and will be reported 
in a later paper. The test section is constructed of 
plexiglass to facilitate flow visualization. The side- 
wall angle a and swirl vane angle C ^re variable. The 
sidewall angle is set by inserting one of three blocks 
with a sidewall angle a of 90, 70 or 45 deg. Tlie swirl 
vane angle is continuously variable with discrete values 
of ^ ■ 0, 45, 60 and 70 deg. being chosen f-»r the pre- 
sent study, with Re^ ■ 1.26 x 10*. 

The Expe rimental Te chniq ue 

The basic technique for the experiment discussed 
in this paper involves the visualization of individual 
neutrally-buoyant helium-filled soap bubbles. The 
bubbles which are approximately 0.5 to 1 am in diameter 
are produced by a Sage Action, Inc. generator-injector. 
The bubbles are produced in an injector device which 
can be located at various positions in the flowfield. 
Because of the substantial size of the injector itself 
(1.5 cm in diameter and 6 cm long) it is not inserted 
directly into regions of flow interest. Instead it is 
either mounted upstream in the stilling chamber, or 
flush mounted to the wall of the large diameter pipe of 
the test section, to inject bubbles directly into the 
comer recirculation zone. The bubbles are illuminated 
by a beam of light from a high pewer 35 rm slide pro- 
jector which is located downstream of the test facility. 

I n t e rpr e t a_t i on o f^ H ow Visualization 

Photographs of illuminated bubbles are taken with 
various shutter speeds and camera positions. At rela- 
tively long exposure times (such as 1/6 sec) a series 
of streaks are visible in the field of view correspond- 
ing to pathlines of individual bubbles. A 35 mm single 
lens reflex carv'ia with a 45 ram lens was used. Tri-X 
Fan black and white film with a normal ASA of 400 was 
used and developed with a special process which pushed 
the ASA to 5000. In relatively lower turbulence inten- 
sity portions of the flowfield mean flow directions can 
be obtained by ensemble averaging local tangents to path- 
lines traced out by soap bubbles. This helps define the 
flovfield geometry in terms of the outline of recircu- 
lation regions. The outlines of these recirculation 
regions constitutes the best check we have on the reli- 
ability of the data. 

RESL0.TS AND DISCUSSION 

Tlie predictions presented here are computed simula- 
tions of the Isothermal airflow in axisymroet ri c combus- 
tor geometries. As noted earlier diameter expansion 
ratio D/d is 2.0, inlet Reynolds number Re^j ■ 1.26x10*, 
wall expansion angle a • 90, 70, and 45 deg. and swirl 
vane angle C varies from 0 to 70 deg. All results are 
obtained via a nonuniform grid system due to the result- 
ing enhancement of solution accuracy. Nineteen cells 
are employed in the r-dlrecilon, and they are clustered 
near the shear layer region and along the wall and 
centerline. Cells in the x-direclion are gradually ex- 
panding in the streamvise direction, and from 18 to 31 
are employed as required to produce the desired a. The 
inlet profiles of axial velocity u and swirl velocity w 
are idealized as "flat" (i.e., constant-valued). 

Axial and swirl velocity profiles are presented for 
the three a values and four C values, clearly showing 
details of these parameter influences. Then, streamline 
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Fig. 5. Predicted axial velocity profiles for * 1.26 x 10* showing 
the effect of swirl vane angle ^ for wall expa.ision angles a: 
(a) 90* , (b) 70*. and (c) 45*. 

[ 0—0 > Experiment (30) with ^ » 0 and * 2 x 10*1 


plots for each value of t, calculated and drawTi by com- 
puter for each i, allow comparison of the shape and size 
of recirculation zones. Further, the relationship be- 
tween the length of the comer recirculation zone CR2 
and 3 is seen, as well as that for the central toroidal 
recirculation zone CTRZ. 

Discussion is pr^narllv al3u*d at guiding designers 
in judlclouslv choosing where experimental emphasis 
should be placed and'or in interpolat Ing results from a 
limit *d amount of e.xperlmental data. For coinparlson 
with other results, it should be noted (6) that swirl 
number S and 0 are related by S ■ 2/3 tan so that 
vane ingles 45, 60 and 70 deg, for example, correspond 
to S Values of 0.67, 1.15 and 1.33, respectively. 

oclty Prof H « 

Predicted mean axial velocity profiles for the a ■ 
90, 70, and 45 deg. geometries are shown in Figures 5a, 

5b and 5c above. The nonswirling case (J - 0) exhibits 
good qualitative agreement with measurements of Chatur- 
vedl (30) in a geometrical Iv similar facility. 

The Influence cf 5 Is most drasvitlc near the com- 
bustor inlet, where the 5-0 profile in Fig. 5a shows 
a Urge comer recirculation region, provoked by the 
sudden increase in fh>w area. A very large value of 
3u/3r occurs near the InUt, which Is indicative of high 


turbulence energy generation In a strong shear layer. 
Further, the nonswirling centerline v«»or'fry exhibits 
little change In the streaciwise direc**^^" 

With ^ - 45 deg. the mean axial velocity profile 
is dramatically changed. Near the inlet a central toroi- 
dal recirculation zone appears and the comer recircula- 
tion zone is considerably smaller. Also, a maximum 
velocity value occurs in an annular fashion near r/D ■ 
0.25, although a more flattened shape quickly develops 
before x/D « 1.0. It should also be noted that the 
boundary layer on the outer sidewall Is relatively thin 
and is not easily seen on the figures using the chosen 
scale. 

The strong swirl case of • 70 deg. shows a much 
larger central recirculation region at x/D « 0.3, which 
is caused by strong centrifugal effects. This promotes 
a very high forward velocity near the wall rather than a 
comer recirculation region. This tendency has been 
qualitatively observed by combustor designers at high 
degrees of swirl, but little quantitative data are yet 
available to precisely substantiate this phenomenon. 

The radial extent of the central recirculation zone and 
the velocity near the vail quickly diminish downstream 
as swirl strength Is dissipated. 

The effect of a la also most visible near the inlet, 
where only the strong swirl profile exhibits a noticeable 
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difference between the a ■ 90 deg. and 70 deg. geoiae- 
triea. For a • 70 deg. this profile is more flat In the 
central region. However, comparison of velocity profiles 
for a ■ 45 deg. with those for a ■ 90 deg. reveals that 
the 45 deg. vjne prediction is also influenced, .,s no 
comer recirculation zone is present for this swirl case. 
Also the off-axis velocity maximum is located closer to 
the centerline. The 70 deg. vane case is considerably 
affected in that the center recirculation bubble has 
less radial extent and the velocity near the wall is 
greatly reduced. It should be noted, however, that only 
a slight effect of a (for the range considered here) re- 
mains beyond x/D values of about l.O. 

Figures 6a through 6c show swirl velocity profiles 
for the same three geometries with vane angles of 45, 60, 
and 70 deg. All of these profiles show solid-body-rota- 
tion behavior near the centerline, even near the inlet 
where a flat profile is a specified inlet conditior. 

The radial location of the station maximum for w tends 
to Increase with x/D in Fig. 6a. Also irregularities 
of the profiles at x/D ■ 0.32 disappears with increasing 
X. Hence, swirl as well as axial velocity profiles ap- 
pear to approach those corresponding to swirling flow in 
a rlp« (32) as X increases. 

Comparing Fig. 6b with 6a, it is the weaker swirl 
cases which show the rest appreciable effect of a at 
x/D ■ 0.32, and again this diminishes with increasing x. 


The 70 deg. vane profiles for Q ■ 90 and 70 deg. at 
x/D •• 0.37 qualitatively similar, although the maxi- 
mum velocity for <1 ■ 90 deg. is 0.9 w/uq as opposed to 
0.8 v/U(, for a - 70 deg. In Fig. 6c it is the strong 
swirl case which is most influenced by a, with a maximum 
value of 1.0 for w/uq at x/D ■ 0.32. The profiles at 
x/D • 1.68 exhibit similarity for each geometry, and if 
w/wg^ is plotted against r/D the w-curves will collapse 
to a single characteristic curve. 

S t r^a mline Plots 

Figures 7a th.'ough 7d are calculated and plotted by 
computer to show the combustor designer the sequence of 
predicted streamline patterns he should expect upon in- 
creasing the vane angle for the ^ ■ 90 deg. combustor. 

The recirculation zones are important to designers (10) 
because most of the burning occurs in these regions, 
and they exhibit the highest temperatures (11). The non- 
swirling case shown in Fig. 7a exhibits a large comer 
recirculation region as indicated also in Fig. 5a. As 
swirl is introduced, a central recirculation region ap- 
pears in conjunction with a decrease In size of the 
comer recirculation region, as seen in the 45 deg. 
swirl angle case. Further increases in swirl vane angle 
result in continued enlargement of the center zone. 
Similarly, the comer recirculation rone is gradually 
reduced in axial extent until it disappears by C * 70 deg. 
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Fig. 6. Frrdlcied swirl velocity profiles for Re^ ■ 1.26 x 10* shoving 
•he effect of swirl vanf ai*gle for wall expansion angles a: 
(a) 90*, (b) ?C*, and (c) 45' 


762 


os - IS 




S 0.50 


0.25 



c) > = 60" 



Axial Position x/D 

Fig. 7. Predicted streanllne plots for a wall ex- 
pansion angle 1 • 90* and Re^ - I 26 x 10* 
for various swirl vane angles 
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Fig. 3. Predicted streamline plots for a wall ex- 
pansion angle i * 70* and Re^j * 1.26 x 10* 
for various swirl v.ine angles 


The same series of streamline patterns is displayed 
in Figs. 8a through 8d for the 70 deg. expansion geome- 
try. The same trend Is found as the vane angle Is In- 
creased, except comparison with corresponding figures 
for the 90 deg. expansion reveals that the center recir- 
culation bubble Is slightly larger In both radial and 
axial extent. The effect of vane angle Is again the 
same for the deg. expansion combustor as shown in 
Fig. 9 parts (a) - (d) . In this case the center reclr- 
culatlon rone Is generally slightly smaller in both 
directions than the 90 deg. geometry. The combustor 
designer mav obtain further insight by observing a simi- 
lar series of streamline plots predicted by Novlck et al 
(10) for an isothermal dump combustor flowfield with the 
following differences: a 90 deg. expan.«ion, an inlet 

hub and a constricted exit. 

An indication of the validity of the computations 
can be obtained bv comparing predicted mean streamline 
patterns with the pathlines traced out by the soap bub- 


bles in the flow visualization experiments. A sample 
flow visualization photograph is presented in Fig. lOa 
corresponding to the zero swirl, 90 deg. expansion angle 
flowfield. The photograph, taken with a relatively long 
time exposure (1/8 sec.) where the flow is from left to 
right clearly shows a great number of Individual path- 
lines. Photographs of this type can be used to distin- 
guish regions of highly turbulent flow from smoother 
regions (for example near the centerline of the flow- 
field) which have smoother, stralghter pathlines. In 
addition the outline of the corner recirculation region 
can be estimated from Fig. IGa (and numerous additional 
photographs taken at the Identical run condition). For 
this geometry the mean stagnation point defining the 
end 01 the recirculation zone appears to be at about 2 
chamber diameters (8 step heights) downstream of the 
sudden enlargement. Examination of Fig. 7a shows that 
Chit Is in quite good agreement with the numerical pre- 
diction. 
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of Che centr^il recirculation region, an Inaccuracy which 
has also been encountered by others in red rculat ion 
zone length prediction with the standard k-c turbulence 
mcdel ^ especially under svirl flow conditions (21,27). 
There are a number of possible reasons for this, but 
Che Dost likely one is chat a turbulence model whose 
ba;Ls and parameters are adequate for simple flow sicua- 
dons is not adequate to handle the tnore complicated 
swirling recirculating flow situation. 

CLOSURE 

The present research Is concerned with specific 
sub-problems of the conplete three-dimensional problem 
and restricted initially to axi symmetric geometries un- 
der low speed and nonreacting conditions. These are 
areas in which there is a need for more fundamental re- 
search (with practical applications) and in which the 
prime gas turbine engine manufacturers would like to 
see developments. 

Many factors affect the existence, size and shape 
of the corner recirculation zone and central toroidal 
recirculation zone (10,11). A major outcome of the 
current work is the ability to characterize and predict 
more realistically than previously the existence, size 
and shape of the comer and central recirculation zones 
as a function of the angle of the sloping wall, the 
degree of swirl imparted to the Incoming flow, the ex- 
pansion ratio, and other swirler and geometric para- 
meters. In order to accomplish these goals, it was ne- 
cessary to include the development of a suitable solu- 
tion procedure in the form of a computer program. Com- 
puter code embellishments include a satisfactory techni- 
que for handling the irregular boundary. A few para- 
meter variations were Investigated computationally in 
order to make combustor designer information available 
in a directly usable form. 

Comparisons of the calculated results yielded en- 
couraging results as the gross features of the flow- 
fields were quite well predicted. A problem in swirl- 
ing flows is the accuracy with which the sizes and 
shapes of the recirculation zones may be predicted. 

This may be partially attributed to the quality of the 
turbulence model. Further research should emphasize 
turbulence model development for swirling recirculating 
f I ows . 
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Abstract 

Numerical compuCacions have been undertaken 
for a basic 2-D aKisymmetric flowfield which is 
similar to that foisid in conventional gas turbine 
and ramjet combustors. A swirling flow enters a 
larger chamber via a sudden or gradual expansion. 
The calculation method involves a staggered grid 
system for axial and radial velocities, a line 
relaxation procedure for efficient solution of the 
equations, a two-equ.it ion k-C turbulence model, a 
stairstep boundary representation of the expansion 
flow, and realistic accommodation of swirl effects. 
The results include recirculation zone characteri- 
zation and predicted mean streamline patterns. 
Predictions with and without chemical reaction are 
obtained. An associated isothermal experimental 
flow study is providing a useful data base. Suc- 
cessful outcomes of the work can be incorporated 
into the more combustion- and hardware-oriented 
activities of industrial concerns. 

K In t roduc t i on 

U_L . The P rob lem 

The basic ramjet combustor configuration is 
the so-called sudden expansion dump combustor with 
corner recirculation zones. In this type of com- 
bustor, liquid fuel is sprayed into the ram a^'. 
upstream of the dump station, although it may also 
be injected directly Into the chamber via side-wall 
Inlets. Primary flank> stabilization is provided by 
the flow recirculation regions, which may be sup- 
plemented, at the expense of total pressure loss, 
with mechanical flameholding devices at the air 
inlet/combustor interface and/or the presence of 
inlet air svdrl, obtained by the use of tangential 
injection or swirl vanes. The flow throughout is 
multiphase, subsonic, turbulent and involves large- 
scale comer recirculation zones. With strong 
swirl in the inlet flow a central toroidal recir- 
culation zone (a recirculation bubble in the mid- 
dle of the chamber near the inlet) also presents 
itself. Even gross features of the flow are not 
known quantitatively with certainty: for example, 

factors affecting the existence, size and shape of 
the recirculation zones.* 

The designer has a formidable problem in aero- 
thermochemistry, and the modeling task is to pro- 
vide a route which leads to the accomplishment of 
design objectives more quickly and less expensively 
than current practice permits. Some combustor 
modeling problems are: 

I. Physical processes - turbulence, radia- 
tion, combustion and multiphase effects. 

* Aas&^clate Professor, School of Mechanical and 
Aerospace Engineering, Associate Fellow AIAA 
** Graduate Student, School of Mechanical and 
Aerospace Engineering 
*** instructor. Department of Mechanics 
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2. Computer programs - 0-, 1-, 2- and 3-D ap- 
proaches in steady state and transient 
cases. 

3. Unresolved problems - 

Irregular boundary representation 

effect of swirl, recirculation and 
wall proximity on turbulence 

turbulence - reaction interaction 

multiphase simulation 

Current combustor design and development problems, 
the needs of the combustion engineer in practice 
and proposed research tasks which will assist in the 
attainment of design objectives are becoming clear. 
Improvements and new developments (both experimental 
and theoretical) can and should be made, theoretical 
modeling being aided by specific carefully chosen 
experiments being performed. 

In combustion chamber development, designers 
are aided by experiments but, as a suppleoient to 
them, economical design and operation can be greatly 
facilitated by the availability of prior predictions 
of the flowfield. These may be obtained by use of a 
mathematical model incorporating a numerical finite 
difference prediction procedure. This work combines 
the rapidly developing fields of theoretical combus- 
tion aerodynamics and computational fluid dynamics, 
and its improvement and use will significantly in- 
crease understanding and reduce the time and cost of 
development.* * 

The present research work is concerned with 
complementary experimental and theoretical studies 
and is described in a recent paper. ^ Th«> problem 
being investigated experimentally is cjnceried witr 
steady turbulent flow in axi symmetric geometries, 
isider low speed and nonreacting conditions - a study 
area highlighted recently" ** as a fundau;‘'ntal re- 
search requirement in combustion modeling. Predic- 
tions are being made both for nonreacting and react- 
ing flows. The particular problem is concerned with 
turbulent flow of a given turbulence distribution in 
a round pipe entering an expansion into another 
round pipe, as illustrated in Fig. 1. The in-coming 
flow nay possess a swirl component of velocity via 
Pfssage through swirl vanes at angle t (equal ap- 
pijximately to l"i** (w^„/u< „) 1 , and the side-wall 
may slope at an angle a, to the main flow direction. 
The resulting flowfield domain nay possess a cen- 
tral toroidal recirculation zone CTRZ in the middle 
of the region on the axis, in addition to the pos- 
sibility of a comer recirculation zone CRZ near 
the upper comer provoked by the rather sudden en- 
largement of the cross-sectional area. Of vital 
concern is the character! zat ion of flows of this 
type in tersis of the e f fects of side-wall angle a, 

degree of swirl turbulence intensity k. of the 

in 


J 


lnl#t stream and expansion ratio D/d on the result- 
ing flowfield in terms of its time-mean and turbul- 
ence quantities. Such problems have received 
little attention, yet there is a definite need for 
work in this area even under nonreacting flow con- 
ditions. 
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being complemented by an associated experimental 
study. Consideration is given to recent work in 
the finite difference solution, via a primitive 
variable code, of axlsyimer ric swirl flow in the 
combustor geometry of Fig. I, where the inlet ex- 
pansion sidewall may slope obliquely to the central 
axis. Thus a systematic parametric investigation 
may be contemplated on the effect of sidewall angle 
Cl and degree of swirl on the resulting flowfield 
produced. Comparison with available experimental 
turbulent flow measurements assists in confirming 
the final predictive capability. One such experi- 
ment, ifliderway at Oklahoma State University, is 
concerned with measuring the effects of swirl and 
side-wall angle on the streamlines, mean flow and 
turbulence parameters in nonreacting flow. The 
facility and experimental details are described in 
a recent paper. Use is made of the confined jet 
facility shown schem.iC ically in Figure 2. 
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Fig. 1. The flowfield being investigated. 

1.2 Pj*v_iou8 Work 

A mathematical solution of the flowfield of 
interest should provide results, if possible more 
cheaply, quickly, and correctly than possible by 
other means (for example, experiments on real-life 
systems or models). In order to achieve this, the 
model should simulate the flow in all its impor- 
tant respects (geometry, boundary conditions, phys- 
ical properties of gases, turbulence, etc.) and 
provide a means %^ereby the governing equations may 
be solved. Mathematical models of steadily in- 
creasing realism and refinement are now being de- 
veloped, both in the dimensionality of the model 
(together with the computational procedures) and 
in problems associated with the simulation of the 
phy-icc. processes occurring. Clearly there are 
two areas of difficulty: the simtilation and solu- 

tion. 

Several previous publications discuss these 
problems at length, in terms of practical applica- 
tion,*** turbulence modeling''’* and numerical 
solution of 2-D axisymmetric problems via the 
stream function - vorticity or primitive pressure - 
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velocity approach. 


Whereas the former ap- 


proach, used in the 1968 computer program from 
Imperial College for cx.tmple'*, reduces by one the 
number of equations to be solved and eliminates the 
troublesixne pressure (at the expense of trouble 
with the vorticity equation), the preferred ap- 
proach now is SIMPLE (mnemonic for semi-^mp 1 icit 
Mthod for pressure jinked e^quations) which focuses 
attention directly on the latter variables. Be- 
cause it possesses many advantages, the present 
work has been developed immediately on this new 
technique, the basic ideas of which have been em- 
bodied into the 1974 Imperial College TEACH (teach- 
ing elliptic axisymmetric ^aracter istics heuristi- 
cally) computer program.'^ Advances to this code 
have been made elsewhere and applied to other 


2-0‘- ” and 3-D'* 
combustor designer. 


problems of interest to the 


1 . 3 Asso cia ted Ex perimental Study 

Prediction work in the form of computing the 
flowfield in the configuration Just described is 


1 1 ^ n — j 
d |j ! i i|ii! r-5; 
1^- 


Fig. 2 Schematic of flow facility. 

The facility has an axial flow fan whose speed 
can be changed by altering the drive pulley combina- 
tion. Numerous fine screens and straws produce flow 
in the settling chamber of relatively low turbulence 
intensity. The contraction section leading to the 
test section has been designed to produce a minimisa 
adverse pressure gradient on the boundary layer and 
thus avoid unsteady problems associated with local 
separation regions. The sudden exp.msion consists 
of a 15 cm diameter circular Jet nozzle, exiting 
abruptly into a 10 cm diameter test section as shown 
In Fig. 1. The substantial size of this test model 
provides excellent probe resolution for hot-wire 
measurements which are currently isiderway and will 
be reported in a later paper. The test section is 
constructed of plexiglass to facilitate flow visual- 
ization. The sidewall angle a and swirl v.^ne angle 
$ are variable. The sidrvall angle la set by 
inserting one of three blocks with a sidewall angle 
cx of 90, 70 or 45 deg. The swirl vane angle is con- 
tinuously variable with discrete values of t ” 0, 

45, 60 and 70 deg. being chosen for the present 
study, with Re^ • 1.26 x 10*. Presently, experi- 
ments are underway with nonreacting flow as follows: 

1. Photography of neut rally-buoyant heliuB- 
filled soap bubbles, wool tufts, and in- 
jected smoke, so as to characterize the 
time-mean streamlines, recirculation 
zones and regions of highly turbulent flow. 

2. Five-hole pitot probe pleasure measure- 
ments, so as to determine tlM-mean 
velocities u, v and w. 
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3* Singl»-wlrc ilx-oritfritaCion and doubl«* 
wirc hot-wirv mcaturcm^ntf, so aa to 
daduc* tlm«-B»an wlocitias, u, v and w, 
normal ttrcttas u'|», and 

and shear sCrassss vTv'^. uV" and 


Extansiva rasults will ba raportad in a latar papar. 
1.4 Out Una of t ha Papa r 

Tha prasanc papar discusses Cha pradlctlon 
Cachnlqua and axhlblts both Isotharmal and react* 
ing flowflald predictions. Reacting flow predic- 
tions are made on Cha basis of a simple homogeneous 
stoichiometric pramlxed mixture entering and react- 
ing In the chamber. Tha demonstration rasults are 
assessed in terms of their use in the combustor 
design situation. 

2j The Ma t h em a t i c aj P rob I am 


2j I _ The Gove rn Equa 1 1 on s 


The turbulent Reynolds aquations for conser- 
vation of mass, momentum (In x, r and ^ directions), 
turbulence energy k, turbulence dissipation rate C, 
stagnation enthalpy h, fuel mass fraction m^ and 
mixture fraction f ■ (where s Is the 

stoichiometric oxygen/fuel ratio) which govern the 
2-D axis ymme trie, swirling, steady flow may be 
taken as previously.**"^' The transport equations 
are all similar and contain terms for convection 
and diffusion (via turbulent flux terms) and source 
$4 of a general variable (which contains terms de- 
scribing the generation (creation) and cons< spt ion 
(dissipation) of ^) . Introducing turbulent ex- 
change coefficients and the usual turbulent dif- 
fusion-flux (stress - rate of strain type) laws. 

It c^ be shown that the similarity between the 
differential equations and their diffusion rela- 
tions allows them all to be put in the common form: 


Ti li’ * T? fl’' ■ 


( 1 ) 


where t is the general dependent variable. The 

forms of the source term S. are given in earlier 
. 2 1 * 
wore. 


loq)liclt here is the use of the two-equation 
k-C turbulence model****' and constant Prandtl- 
Schaldt miBbers from which exchange coefficients 
are calculated: 



"fu ■ 
"fu ■ 



The value of the latter model Is that It tries to 
assert the effect of turbulence on the reaction, the 
time to heat up the premixed mixture by eddies of 
hot combustion product being related to the rate of 
dissipation of the fuel concentration. In many 
cases, this Is the limiting factor In controllirg 
the reaction, and not the final stages of molecular 
processes governed by an Instantaneous Arrhenius 
expression. The variable g (mean square fluctuating 
component of fuel concentrat ion) may be obtained 
from Its governing differential equation or, with 
the assumption that generation equals dissipation. 
Its reduced algebraic equation. In the present re- 
acting flow calculations, its algebraic form is used 
with reaction rate specified as the average of the 
values given by Eq. (3). Two important problem 
areas currently being Investigated are: 


1. More realistic simulation of complex 
chemistry and pollutant emission in 
turbulent reacting flows. 


2. Inclusion of radiation heat transfer 
effects. 

These equations have to be solved for the tlaie- 
mean pressure p and velocity comptinents u, v, and 
w, etc. Then other useful designer Information, 
like streamline plots, breakaway and rear tachsK*nt 
points, recirculation zones and stagnation points, 
for example, may be readily deduced. So also may 
temperature, flame zones and regions of ^burned 
fuel be deduced. Earlier publications***'* provide 
details of the present simulation and solution tech- 
nique, and only highlights of the primitive-variable 
approach need be given here in the context of the 
specific problem being Investigated. 

2 . 2 The F In i t e Pi f feren ce Fjj rau be t ion 

Solution may be via the strejn funct lon-vortlc- 
Ity, or primitive pressure-velocity approach. 

Whereas t.ie former approach, used in .he 1968 com- 
puter program from Imperial College for example**, 
reduces by one the number of equatizms to be solved 
and eliminates the troublesome pressure (at the ex- 
pense of trouble with the vortlclty equation), the 
preferred approach now Is SIMPLE (r«>emonlc for Mmi- 
i.sipliclt Mthod for £ressure linked equations) which 
fccuses attention directly on the latter variables. 


g • C^ok'/e ♦ 

r, ■ u/o^ 


( 2 ) 


Implicit also is the use of a simple exothermic 
three -component premixed chemical reaction charac- 
terized by equations for h, m , and f, %diich Is 
a conserved property. Fuel afU oxidant may coex- 
ist at a point In the flow, and the consimiptlon 
rate of fuel calculated depending on 

local species concent rat ions, temperature, and 
turbulence levels. For this, one can use either 
the time-averaged Arrhenius model or eddy-breakup 
reaction model: 


Because It possesses m^«ny advantages, the pre- 
sent work has been developed immediately on this 
new techniques, the basic Idea^ of which had been 
embodied into the 19 74 lsq«crlal College TEAiTi 
(Reaching eliptlc Mlsymmetrlc characterl stlc Heur- 
Istlcally) computer program**. A primitive pres- 
sure-velocity variable, finite difference cosipuier 
code his been developed to allow the computation of 
nonreacting turbulent swirling flow aerodynamics 
within the second stage nozzle section of the axl- 
sysssetrlc combustor. The method and program Involve 
a staggered grid system for axial and radial veloc- 
ities, and a line relaxation technique for efficient 
solution of the equation** '*. Turbulence s/mula- 
tlon Is by way of two-enuetlon k-€ model. The pre- 
sent analysis has been developed based on Ideas 
which had been embodied into the TEACH cos^uter pro- 
gram. Several techniques are Incorporated to en- 
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hancc #ccuracy and cfflc'.mc/, and a cowplate da- 
acrlDtlon of tha finally davalopad coda fa avail* 
•bU'* 

Tha flnlta diffaranca aquations ara aolvad on 
a cooplax Besh aystaa, ahovn In fig. 3. Tha intar* 
aactiona, the point P for axsopla, of tha aolld 
llnaa aark tha grid nodca tdiara all varlablaa ax- 
capt tha u and v valoclty coaponanta ara atorad. 

Tha lattar ara atorad at points ifhlch ara danotad 
by arrows (and labalad w and a raapactivaly) lo- 
catad aildvay batwaan tha grid Iniarsact Iona • Da- 
tails of tha spaclal narlta of this sta|aarad grid 
systan hava barn raportad pravloualy. * ^ Tha 
diffarant control volmcs C, U and V which ara ap- 
propriata for tha P, w and a locations raapactivaly 
ara also glvan in tha flgura. Thus tha u- and v- 
call boundaries lla althar exactly on or exactly 
halfway batwean tha grid Intaraactlon points P. 

This laprovaa tha accuracy of calculating pressure 
gradient affects on tha u and v valocltlas (as 
pressure values now lie directly on appropriate 
call faces j and of calculating certain call bound- 
ary fluxes. Tha flowficld dooialn bomdarlas ara 
positional aidway between grid points so that they 
coincide with normal velocities as in Fig. 4, which 
shows an axanpla of grid specification for tha 
gaoaetry under consideration. 


a. V tre 
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Fig. 3 The three control volumes associated 
with points of tha three grids. 



Fig. 4 An axasg>la of a coarse grid system 

L'aing employed to fit the flow domain 

In obtaining tha PDFs in tha primitive ap- 
proach it is commim to uea a hybrid formulation, 
which is a combination of tha central and upwind 


differencing schema; it is applied to tha represen- 
tation of the convection and diffusion terms and 
provides, whan needed the required degree of up- 
stream differencing. In this manner tha PDF for a 
typical variable becomes 


♦j ^ *5 


(*) 


whe re 


• r, •* - s* 
p J J p 


• sun over N. S. E and U neighbors 


for each variable ^ at each point P of a grid sys- 
tem covering the flow domain. 


2.3 The Solution Procedure 


The flowfield is covered with a nonuniform 
rectangular grid system; typically the boundary of 
tha solution domain falls halfway between its im- 
mediate nearby parallel gridlines; and clearly 
specification of the x and r coordinates of the 
gridlines, together with information concerned with 
the position of the upper and curved boundary, is 
sufficient to determine the flowfield of interest. 
The ct *’ved or sloping boundaries are simulated by 
means of a stairstep approach with sloping I oun- 
dary segments where appropriate. 

finite difference equations like Eq. (4) hold 
for all variables at all internal points of the 
mesh system, except those adjacent to wall tcun- 
daries. At these points, the correct boisidary in- 
fluence is inserted by a linearized source tech- 
nique. using wall functions so as to avoid detailed 
calculations in near-wall region. Tlie effect of 
swirl on wall function specification is handled as 
follows. The previous ideas are extended to find 
tota_l tangential wall shear stress near boundaries 
TTnvolving X and 9 directions for an r • constant 
wall, and r and 9 directions fot an x ■ constant 
wall). Then appropriate c omponents are deduced 
directly (for the u and w velocities which are tan- 
geritial to an r • constant wall, and v and w velo- 
cities which are tangential to an x • constant 
wall). The effects on u. v and w momentum equations 
are incorporated via the usual linearized source 
technique. ** 

The finite difference equations and boundary 
conditions constitute a system of strongly-coupled 
simultaneous algebraic equations, though they ap- 
pear linear they are not since the coefficients and 
source terms are themselves fiaictions of some of 
the variables, and the velocity equations are 
strongly linked through the pressure. In the solu- 
tion procedure algebraic equations like Eq. (4) are 
solved many times, coefficient and source updating 
being carried out prior to each occasion. The 
practice uaed here la to m.ske uae of tiia well-known 
tri-diagonal matrix algorithm (TOHA). wiiereby a sat 
of equal iona. each with exactly thiee unknowns in 
a particular order axcept the first and last which 
have exactly two unknowris. may be solved sequenti- 
ally. In the 2-0 problem one conaloers the values 
of gridpolnta along a vertical gridline to be un- 
known (values at f. M and S for each point P) . hut 
take as known, most racenr. valuaa being used, the 
values at each C and U neighbor. The TOKA la then 
applied to this vertical gridline. In th^t manner 
one can trawerae along all Unas in the vertical 
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dlrvction ■«qu#nctally fron l«ft to right of who 
intcgrotlon donutn, using « cortsln iSogroo of 
undor rolsRst ioii to onhanco convorgonco. 

Tho »aln o«phastt new is on thr application of 
th# ccMsputsr coda to Idoalifod ar isysmiotr i£ cosi' 
bustor gooo-rtrios to assist in design and develop* 
■lent studies. Soth nunreacting and stoichiosierr ic 
preoixed »echane/air reacting flows are considwied, 
with idealised wall cooling to lOOOK in the latter 
case, so as to present calculations in increasing 
order of cosipleKity. Velocity profiles, streamline 
plots and temperature profiles are exhibited. 
Finally, recirculat ion tones are discussed In terms 
of the effects of swtrl strength, sidewall angle 
and chemical reaction on their sites. Discussion 
is primarily aimed at guiding designers in judf* 
cioiisly deciding where exper iisenta 1 emphasis should 
be placed and/or in interpolat ing results from a 
limited amount of experimental data. 


Gimputat ions shi«wn include geometric and vane*swirl 
angle variations. For comparison purposes. It le 
worth noting that swirl number S and # are related 
by S • ;/} tan g, so that vane angles 45, 60 and 
70 deg. for example, correspond to S values of 0.67, 
1.15 an«< 1.8.*, respectively. 

As rioted earlier diameter expansion ratio 
D/d is 2.0, inlet Reynolds number Re^ • 1.26 x 10^, 
wall expansion angle a • 90, 70, and 45 deg. and 
swirl vane angle 6 varies from 0 to 70 deg. All 
results are obtained via a nonunifom grid system 
due to the resulting enhancement of solution ac* 
curacy. Nineteen cells are employed in the r* 
directicu, and they are clustered near the shear 
layer region and along the wall and centerline. 
CelH in the x*directlon are gradually expanding 
in the streamwise dirw»ctiun, end from 18 to 31 are 
employed «a requf red to produce the desired ci. The 
inlet profUes of axial velocity u and swirl velo* 
city w are idealised as **flst** (i.e., constant* 
va lued ) . 




(a) 6 • 0 dpq 


(a) • 0 deg 



0.00 1.00 

(b) 6 • 4S deg 



0.00 1.00 

(b) ■ 4S dtrg 




u.oo i.eo 

(c) ♦ • 60 «le 9 


O.OC 1.30 2.00 3.00 

(c) ♦ • 50 dpq 



0.00 1.00 2.00 3.00 
rtXIRL POSITION X/D 

(d) ♦ • 70 de. 


0.00 1.00 2.00 3.00 
flXJRL POSITION X/0 

(d) • ■ 70 d«9 


Fig. 5. Predicted streamline plots for a wall ex* 
pension angle 3 • 90* for varioua awirl 
vane angles |, nonreacting tlowfield. 


Fig. 6. Predicted streamline plots for a wall es- 
panslon angle 3 ■ 45* for varioua swirl 
vsne angles t, nonreacting flowfield. 



3. 1 Nonreacting Flow 

The predictive capability of the code la now 
demonstrated for the flowfleld of Fig. 1 under non- 
reacting conditions. Figures 5a through 5d are 
calculated and plotted by computer to show the com- 
bustor designer the sequence of predicted stream- 
line patterns he should expect upon Increasing the 
vane angle for the a « 90 deg. combustor. The 
recirculation zones are Important to designers^ ^ 
because most cf the burning occurs In these regions, 
and they exhibit the highest temperatures. The 
nonswlrllng case <«hrv.\ ^n Fig. 5a exhibits a large 
corner recirculation region. As swirl Is Intro- 
duced, a central recirculation region appears In 
conjunction with a decrease In size of the corner 
recirculation region, as seen .a the 45 deg. swirl 
an'^le case. Further Increase** in swirl vane angle 
result In continued enlargement of the center zone. 
Similarly, the corner recirculation zone Is grad- 
ually reduced In axial extent until It disappears 
by ^ ■ 70 deg. 

The same series of streamlli. . atterns Is 
displayed In Figs. 6a through bd tor tie 5 deg. 
expansion geometry. The sa.ns trend Is found as 
the vane angle Is increased except that In this 
case the center recirculation zone is generally 
slightly smaller In both directions. The combustor 

•/■ • 0.1 


designer may obtain further Insight by observing 
a similar series of streamline plots predicted by 
Novlck et al^^ for an Isothermal dump combustor 
flowfleld with the following differences: a 90 deg. 

expansion, an Inlet hub and a constricted exit. 

An Indication of the validity of the computa- 
tions can be obtained by comparing predicted mean 
streamline patterns with associated experimental 
studies In progress at Oklahoma State University 
as discussed In Section 1.3. Extensive details 
will be reported later. Suffice It to say here 
that the agreement between predictions and avail- 
able flow visualization Is excellent.^ Currently, 
velocity and turbulence measurements are In pro- 
gress. 

Predicted mean axial velocity profiles for the 
a • 90, 70, and 45 deg. geometries are sho%ni In 
Figs. 7a, 7b and 7c. The nonswlrllng case ■ 0) 
exhibits good qualitative agreement with measure- 
ments of Chaturvedl^^ In a geometrically similar 
facility. The Influence of ^ Is most dramatic near 
the combustor Inlet, where the ^ ■ 0 profile in 
Fig. 7a shows a la**ge corner recirculation region, 
provoked by the sudden Increase In flow area. A 
very large value of 9u/3r occurs near the Inlet, 
which Is Indicative of high turbulence energy gen- 
eration In a strong shear layer. Further, the 

•/• • l.t I/I • I.M 
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Fig. 7. Predicted axial velocity profiles showing the effect of swirl vane angle ^ for wall 
expansion ai.gles a: (a) 90*, (b) 70*, and (c) 45*, nonreactlng flowfleld 

(o e- o Experiment’* with ^ • 0) 
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nonswlrllng centerline velocity exhibits little 
change in the stre.imwise direction. With • 45 
deg. the aean axial velocity profile is dramati- 
cally changed. Near tht inlet a central toroidal 
recirculation cone appears and the corner recircu- 
lation tone is considerably smaller. Also, a suixi- 
■urn velocity value occurs in an annular fashion 
near r/D • 0.25, although a more flattened shape 
quickly develops before x/D • l.O. It should also 
be noted that the boundary layer on the outer side- 
wall is relatively thin and is not easily seen on 
the figures using the chosen scale. The strong 
swirl case of ^ ■ 70 deg. shows a much larger 
central recirculation region at x/D • 0.3, which 
is caused by strong centiifugal effects. This 
promotes a very high forward velocity near the 
wall rather than a corner recirculation region. 

This tendency has been qualitatively bserved by 
combustor designers at high degrees of swirl, but 
little quantitative data are yet available to pre- 
cisely substantiate this phenomenon. The radial 
extent of the central recirculation zone and the 
velocity near the wall quickly diminish downstream 
as swirl strength is dissipated. 


The effect of a is also most visible near the 
inlet, where only the strong swirl profile exhibits 
a noticeable difference between the a ■ 90 deg. and 
70 deg. geometries. For a • 70 deg. this profile 
is mi>re flat in the central region. However, 
comparison of velocity profiles for a • 45 deg. 
with those for a • 90 deg. reveals that the 45 deg. 
vane prediction is also influenced, as no corner 
recirculation zone is present for this swirl case. 
Also the off-axis velocity maximum is located 
closer to the centerline. The 70 deg. vane case is 
considerably affected in that the center recircula- 
tion bubble has less radial extent and the velocity 
near the wall is greatly reduced. It should be 
noted, however, that only a slight effect of a (for 
the range considered here) remains beyond x/D 
values of about l.O. 

Swirl velocity profiles for the same three 
geemetries with vane angles of 45, 60, and 70 deg. 
are shown in Figs. 8a through 8c. All of these 
profiles show sol id-body-rotat ion behavior near the 
centerline, even near the inlet where a flat pro- 
file is a specified inlet condition. The radial 
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Pig. 8. Predicted swirl velocity profiles showing the effect of swirl vane angle ^ for 
wall expansion angles (a) 90*, (b) 70*, and (c) 45*, nonreacting flowfield. 
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lOiMtlon of chft* scat ion for w tends Co in- 

cre.ise with x/D In Fig* 8s. Also Irregular it lex 
of Che profiles at x/D ■ 0.12 disappears with In- 
creasing X. Hence, twirl as well as axial velocity 
profiles appear to approach Chose corresponding to 
swirling flow in a pipe*’ as x increases. Compar- 
ing Fig. 8b with 8a, it is the weaker swirl cases 
which show the ohysC appreciable effect of a at 
x/D • 0.12, and again this diminishes with increas- 
ing X. The 70 deg. vane profiles for a • 90 and 70 
deg. at x/D • 0.32 are qualitatively similar, al- 
though the maximtim velocity for a • 90 deg. is 0.9 
w/u^j as opposed to 0.8 w/u^, for a • 70 deg. In 
Fig. 8c it la the strong swirl case which is most 
influenced by it, with a maxiimim value of l.O for 
w/Ujj at x/D • 0.32. The profiles at x/D • 1.68 
exhibit similar I tv for each geometry, and if w/Wq}^, 
is plotted against r/D the w-curves will collapse 
to a single charac ter 1st Ic curve. 

^._2 ^Re.ictlng Flow 

Ceneial predictions are now given for the 
effect of swirl and sidewall .'rngle on a correspond- 
ing flow with chemical reaction. Figures 9a through 
9c show what happens to the axial velocity profiles 


we • «.i 



when the inlet flow is replaced by a premixed 
stoichiometric mixture (methane gas and air) and 
combustion la allowed to proceed downstream of the 
inlet. This figure may be compared with and con* 
trasted to the corresponding figure without chemi- 
cal reaction. Fig. 7. Both fU>ws have the same 
general velocity profiles. Notice now that under 
reacting conditions the radial locations of maximum 
values occur at approximately the same values as 
before. This is especially true at locations 
further do%nistream from the inlet. The presence 
of chemical reaction decreases the site of both the 
corner and central rec lrcul.it ion tones, with higher 
recirculating velocities. The presence of reaction, 
gas exp.insion, and subsequent increased axial velo- 
city tends to shorten the recircul.it ion tones; alow 
moving gas is accelerated toward the exit. 

Corresponding temperature calculations provide 
further evidence about the reaction tones, heat 
transfer and velocities effects inside the flow 
domain. Figures 10 and 1 1 show temper.iture profll«i 
at axial stations x/D ■ 0. I and 0.75 for different 
inlet swirl strengths with the sidewall angle a • 

45 deg. only. The top w.sll and sloping sidewall 
are cooled in a oianner simulating the effects of 
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Figure 9. Predicted axial velocity profiles showing the effect of swirl vane angle ^ for 
wall expansion angles i: (a) 90*, (b) 70", and (c) 45), reacting f. vfield. 
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Fig. 10. Temperature profiles at x/D ■ 0.1. 
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Fig. ll. Temperature profiles at x/D • 0.75. 

cooling by way of heat transfer, so restricting 
wall temperatures to at most lOOOK. Notice that 
the CR2 and the CTRZ have the highest temperatures 
since most of the burning is taking place in these 
regions. The recirculation of hot combustion pro- 
ducts has a preheating effect in reverse flow zone\ 
a condition not enjoyed in the high velocity region 
emanating from the expansion corner. Hence, a cool 
zone 1s present there and extends downstream from 
this corner. 

Figure 10 illustrates this effect for the 
I ■ 70 deg. case, while regions of no reverse flow, 
in Che 1*0 and 45 deg. cases especially, have 
cool legions. Further downstream, the temperature 
profiles of Fig. 11 show that all burning has been 
completed in the swirling cases, but chat Che cool 
central region of Che nonrecirculating nonswirling 
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case still persists. Here, of course, only Che 
gross effects of combustion are simulated; no 
attempt has been siade to correspond exactly to any 
particular case. There is no droplet spray of fuel, 
for example, while geometric and boundary details 
are '‘.mplified as compared with an actual combustor. 

3.3 Recirculation Zones 


Interpretation of an extensive range of stream- 
line plots allows recirculation zones to be chara- 
terized and parametrized. Figures 12 and 13 illu- 



and swirl vane angle ^ on corner re- "j 

circulation zone CRZ length. 
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Fig. 13. Predicted effect of sidewall angle 
a and swirl vane angle > on central 
toroidal recirculation zone CTRZ 
length. 

(— - nonreacting, reacting] 
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straCe the predicted effects of sidewall angle a 
and swirl vane angle ^ on the corner CRZ and cen- 
tral CTRZ recirculation zone lengths, respectively, 

under both nonreacting ( 1 and reacting ( ) 

flow conditions. Four points oiay be perceived. 
Firstly, recirculation zone lengths are only 
slightly affected by the sidewall angle, as found 
previously** for the CRZ under isothermal nonswirl- 
Ing conditions. Also at these high Reynolds num- 
bers, the CRZ length is Refolds number independent 
at high Reynolds numbers.* **^ As the sidewall 
angle decreases from 90 to 43 deg. the zone lengths 
decrease slightly: the inlet flow is encouraged to 

'stick* to the confining walls. Secondly, the CRZ 
length decreases and, thirdly, the CTRZ length in- 
creases with swirl strength. Both these effects 
are well-known to combustion engineers who strive 
to utilize the recirculation of hot combustion pro- 
ducts and the bluff body effects of the zones as an 
aid to the combustion processes. Finally, the 
effect of combustion is to decrease the lengths of 
both corner and central recirculation zones. This 
has also been observed by others** with the dump 
expansion case onJy being considered, where shorter 
zone lengths, higher recirculating velocities and 
less mass flow recirculated. The presence of re- 
action, gas expansion, and increased axial velocity 
tends to shorten the zones: slow moving gas near 
the downstream end is accelerated toward the exit. 

4. Concl usions 

Fundamental theoretica4 studies with applica- 
tion to practical geometries provide useful simpli- 
fications of the complete three-dimensional problem. 
Many factors affect the existence, size and shape 
of the corner recirculation zone and central toroi- 
dal recirculation zone. A major outcome of the 
current work is the ability to characterize and pre- 
dict more realistically than previously the exis- 
tence, size and shape of the corner and central re- 
circulation zones as a function of the angle of the 
sloping wall, Che degree of swirl imparted to the 
incoming flow, and the presence of chemical reac- 
tion. Numerical computations have been undertaken 
for a basic 2-D axisymmetric flowfield which is 
similar to that found in conventional gas turbine 
and ramjet combustors. A swirling flow enters a 
larger chamber via a sudden or gradual expansion. 

The calculation method involves a staggered grid 
system for axial and radial velocities, a line 
relaxation procedure for efficient solution of the 
equations, a t«M>-equat ion k-C turbulence model, a 
stairstep boundary representation of the expansion 
flow, and realistic accomnodat ion of swirl effects. 
An on-going experimental program is currently pro- 
viding a useful data base for Che nonreacting flow- 
field. Further development and application is pro- 
viding a valuable supplementary technique for de- 
signers of practical combustion equipment. 
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constant 

chamber diameter 

nozzle diameter 

Arrhenius constants 

fuel, oxidant, product 

ousmentum flux, k generation term 

stagnation enthalpy 

kinetic energy of turbulence 

molecular weight 

time-mean chemical mass fraction 

cell Pec let number 

time->rean pressure 

universal gas constant, mass rate of 
creation per unit volume, residual 
source, radius 

swirl number • 2Cq/ (0^6) ^ source term 
(with subscript) 
stoichiometric ratio 
time-mean temperature 

time-mean velocity (in x,r,0 directions) 
axial, radial, polar coordinates 
turbulent exchange coefficient 
turbulence energy dissipation rate 
effective viscosity 
time-mean density 
Prandtl-Schmidt number 
swirl vane angle ( tan“ * (w£|,/u£,^) 1 , gen- 
eral dependent variable 


Subscripts 


EBU eddy-break-up model 

P,N,S,E,W point, north, south, east, vest neighbors 
p constant pressure 
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Abetrect 


Subecrlpte 


A swirling nonreectlng flow enters e lerger 
cbenber vie a sudden or gradual expansion. Six 
flowfield configurations are Investigated with side- 
wall angles a ■ 90 and 43 deg. and swirl vane angles 
^ • 0, 43 and 70 deg. Photography of neutrally- 
buoyant hellun-f Hied soap bubbles* tufts* and in- 
jected snoke helps to characterise the tlae-nean 
streamlines* recirculation zones and regions of 
highly turbulent flow. From the photographic evi- 
dence, It Is found that central recirculation zones 
occur for the swirling flow cases Investigated. 

They extend from the inlet to approximately x/D ■ 
1.7* after which a relatively narrow precessing 
vortex core exists i^ar the axis. Five-hole pitot 
probe pressure measurements allow the determination 
of time-mean velocities u* v and w. At the inlet* 
they exhibit maximum values at approximately r/D ■ 
0.2 In a sharply-peaked annular fashion. An ad- 
vanced computer code equipped with standard two- 
equation k-c turbulence model Is used to predict 
corresponding flow situations a*.id to compare results 
with the experimental data. Although excellent 
agreement is found for the nonswlrllng flow* poor 
agreement in detail occurs for the swirling flow 
cases* especially near the inlet. This discrepancy 
may be attributed to the lack of realism in the k-€ 
turbulence model and/or to inaccurate specification 
of time-mean velocity and turbulence energy distri- 
butions at the Inlet. 
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o value St inlet to flowfield 

1. Introduction 


In gas turbine and ramjet combustion chamber 
development* designers are aided by both experi- 
mental and theoretical studies.^ The present re- 
search work is concerned %rlth such complementary 
studies. The problem being Investigated Is con- 
cerned with steady turbulent flow in axlsysnetric 
geometries* under low speed and nonreacting condi- 
tions - a study area which is fundamental to combus- 
tor modeling. The particular problem discussed In 
this paper is the flow in a round pipe entering an 
expansion into another round pipe as illustrated in 
Fig. 1. The in-coming flow may possess a swirl 
component of velocity via passage through swirl 
vanes at angle ^ from the streamwiae direction* and 
the side wall may slope at an angle a* also from the 
streamwise direction. The resulting flowfield do- 
main may possess a central toroidal recirculation 
zone in the middle of the region on the axis* in 
addition to the possibility of a comer recircula- 
tion zone near the upper comer provoked by the 
rather sudden enlargement of the cross-sectional 
area. Of vital importance is the characterization 
of flows of this type in terms of the effects of 
side-wall angle a* degree of s%/lrl ^* turbulence 
intensity of the inlet stream and expansion ratio 
D/d on the resulting flowfield in terms of Itg time- 
mean and turbulence quantities. Such problems have 
received little attention* although they are funda- 
mental to the physical processes occurring In air- 
craft combustors. 



Fig. 1 The flowfield being investigated. 

A systematic paraomtric investigation Is being 
undertaken on the effect of side-wall angle a arid a 
swirl vane angle # on the resulting flowfield pro- 
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duced. Conparlson with availablt •xptriwntal tur- 
bulent flow meaeurereentt aaaiata in avaluatlrr tur- 
bulence nodcla and the final predictive capability. 
One such experlnent, underway at Oklahoma State 
Univeraity, la concerned with neeauring the effecta 
of awirl and side-wall angle on the atreanlinaa* 
mean flow and turbulence paranetera in nonreacting 
flow. This paper focuses on the Identification of 
recirculation regions froa flow visualitation, mean 
velocity measurements with a five-hole pitot probe • 
and the comparison of experimental results with 
flowfield predictions using an advanced computer 
code equipped with standard two-equation turbulence 
model. 

2j E xperimental Approach 

2.1 Previous Relevant Experimental Studies 


The awirl vane assembly contains tan vanes which 
are individually adjustable for any vane angle 
The diameter is located at the center of the awirler 
with a atreamlinad nose facing upstream. The down- 
stream end is simply a flat face, simulating the 
geometric shape of a typical fuel spray noitle. 

The idealized combustion chamber model is com- 
posed of an expansion block and a long plexiglass 
tuba. The expansion block is a 30 cm diameter disk 
with a 15 cm diameter hole centered on its axis 
through which the air entering the model flows 
smoothly. Tha downstream face of the expansion 
block has been shaped to provide the desired flow 
expansion angle a which is shown in Fig. 1. There 
are currently three interchangeable expansion blocks 
and the appropriate choice gives a ■ 90, 70, or A5 
degrees. 


There have been several experiments performed 
with nonreacting flows in expansion geometries that 
have been reported in the literature, examples of 
which are contained in Refs. 5-8. References 5-7 
also include flowfield predictions, made with ver- 
sions of the TEACH-T computer program.* These ex- 
periments include time-mean velocity measurements 
[with hot-wire and pitot probes and laser Doppler 
anemometry], turbulence measurements [%rlth hot-wires 
and laser anemometers] and flow visualization. The 
majority of the earlier measurements were made in 
nonswirling flows; however some noteworthy experi- 
ments have been made in swirling confined Jets.*** 
Direct comparison between the results of the cited 
experiments and the present experimental results is 
generally not possible because of differences in 
geometry. However, in the nonswirling Jet compari- 
sons were possible with experiments of Chaturvedl,* 
vrho measured mean and turbulent flow quantities 
downstream of a sudden expansion of diameter ratio 
D/d * 2 and various expansion side-wall angles a. 

2.2 T est Facility 

The present experiments have been conducted in 
the confined Jet facility shenm hematically in 
Fig. 2. The facility has ar axis', flow fan whose 
speed can be changed by alteriug the varldrlve 
mechanism. Numerous fine screens and straws produce 
flow in the settling chamber of relatively low tur- 
bulence intensity. The contraction section leading 
to the test section has been designed by the method 
of Morel to produce a minimum adverse pressure 
gradient on the boundary layer and thus avoid un- 
steady problems associated with local separation 
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Fig. 2 SchesMitic of overall flow facility. 

The test section consists of a swirl vane as- 
sembly and an ideslizad combustion chamber model. 


There are no film cooling holes or dilution 
air holes in the test facility, and the chamber 
wall consists of a straight pipe. The inside dia- 
meter of the pipe is 30 cm and the length is approx- 
imately 125 cm. The substantial size of the teat 
model provides excellent resolution for five-hole 
pitot probe measurements and also for flow visuali- 
zation. 

2.3 Flow Visualization Techniques 

Illumination and Photography . A slide projec- 
tor, located downstream of the test facility, serv- 
ed as the light source for illumniation of the 
bubbles. A vertical sheet of light 4 cm thick was 
produced to illuminate the rx-planc of the bubble 
flow pattern. This was provided by using a slide 
which is opaque except for a thin slit cut out for 
light passage to the test section. • Trl-X Pan, a 
very light-sensitive photographic film rated at ASA 
400, along with a large camera aperture of f 2.0, 
was employed in order to obtain photographs of ac- 
ceptable contrast. Morevoer, the film was sub- 
stantially overdeveloped to compensate for underex- 
posure; this is equivalent to using an extremely 
light-sensitive film rated at approximately ASA 
6,000. Photographs of the bubble stresklines were 
typically taken at shutter speeds ranging from 1/60 
to 1/4 second. 

Neutrally-Buoyant Helium-Filled Soap Bubbles . 

Flow visualization experiments using three 
different techniques were conducted during this in- 
vestigation and each is described in the following 
paragraphs. An attractive flow visualization tech- 
nique for both detailed as well as overall charac- 
terization of strongly turbulent flows employs 
neutraiiy-buoyant helium-fi J ied soap hubbies. Al- 
though snich larger than particles of smoke, soap 
bubblea are almost Ideal because they trace parti- 
cle paths and have the average density of air. Of 
all the previous work with soap bubbles, that most 
related to the present study was conducted by Owen 
et al.,*^ who employed them to study the flowfield 
in a cylindrical vortex tube. 

The bubble generator is manufactured by Sage 
Action, Inc.^* It conaiata of the injector, in 
which the bubbles are actually formed, and a con- 
sole that supplies the constituents to the injectot. 
Only the injector is situated in the tunnel. It is 
located far enough upstream of the test section to 
avoid any disturbanca to tha flow pattern under 
study. The bubble film solution ia a special soap 
mixture, which is commercially available from Sage 
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Actioiiv Inc. After foriMtlon, 0.3 to I.O wm dlo- 
■oter bubble! ere Injected into the eurroundlng 
freeetreea, which cerrlee then through the teet 
eection. It le eetlneted^* fron photogrephe that 
thie bubble production rate le aa high aa 300 bub- 
ble! per lecond. 

Snoke-Wlre . Although enoke par tic lee nay not 
be observed Individually, enell end dlacrete etreak- 
llnes can be located within the flowfleld ao that 
snail acale details nay be studied. For exanple, 
sero velocity points irlthln a large recirculation 
zone nay bO' Identified. This la acconpllahed with 
the Mmok€~wirm tmehniqum. alnllar to that used by 
Naglb and his covorkera. ' * * 

The nethod consist a of a fine wire positioned 
in the flowfield, coated with oil, and heated by 
passing an electrical current through the wire. Aa 
the wire is coated with the oil, small beads of the 
oil form on the wire and at each of these beads a 
smoke filament originates when the wire ia heated. 
Cornell^* ^:udied the sister iala used in generating 
'*snoke" for flow visualization, and explained that 
particles created with thia technique are actually 
a vapor-condensation aerosol. That ia, they are 
very small liquid particles near 1.0 x 10*’ oub in 
diameter rather than products of a combustion pro- 
cess. However, in keeping with the literature, 
these particles are referred to as smoke. For the 
mechanical details of the system, the reader is re- 
ferred to Rhode. ^ 

Tufts . Perhaps one of the simplest ways of 
visualizing a flowfield ia by inserting a wire grid 
with a tuft attached at each intersection point of 
the mesh. Besides pointing in the local flow direc- 
tion. their vascillating motion indicates the degree 
of turbulence in steady mean turbulent flow situa- 
tions. Very large tuft fluctuations, however, can 
be p’^oduced by the combined effects of Man flow 
unst. idiness and turbulence. Tuft BcremnM have been 
employed to visualize the vortex shedding behind 
swept wings. They have also been used to observe 
the interaction of an air jet and a cross wind.^^ 
Again. Rhode^ provides details of the tuft appara- 
tus as used in the present study. 

2. A Mean Velocity Heaaurenent 

Many instruments arc used for separately ms- 
surina the magnitude and direction of fluid velo- 
city.^* However, there are only a few instruments 
capable of simultaneously sensing both magnitude 
and direction. One of the simplest of tht.ae ia the 
five-hole pitoi probe which has been developed and 
used by various investigators.***** The particular 
probe employed in this study ia model DC-125- 12-0 
from United Sensor and Control Corp. It is shown 
schematically in Fig. 3 and has a 3.2 am diameter 
sensing tip and shaft containing five tubes. The 
sensing head is hook-shaped to allow probe shaft 
rotation without altering the probe tip location. 

The instruMntat I'm system, in addition to the 
five-hole pitot probe, consists of a manual traverse 
Mchanism. two five-way ball valves, a very sensi- 
tive pressure transducer, a power supply, and an 
integrating voltMter. The differential pressure 
transducer is model 590D from DatsMtrics, Inc. 

The output is read as the d.c. signal from a TSI 
model 1076 integrating voltmeter. 

Prior to production MasureMnts, the five-hole 



Fig. 3 Five-hole pitot probe. 
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Fig. 4 Velocity components and flow direction 
angles associated with five-hole pitot 
measurements (yaw angle 8 in the hori- 
zontal plane and pitch angle 6 in the 
vertical plane. 

pitot probe ia aarodynamically zeroed for yaw, which 
is in the horizontal plane of Fig. 4, ao that x and 
6 axes of the Maaurement coordinate frame coincide 
with those of the test section. The measurement 
procedure for each location within a traverse begins 
irlth aarodynamically nulling the yaw, and determin- 
ing the yaw angle 0. Thie is indicated by a zero 
reading for Py * Pgt tfhere the pressures are identi- 






fled In Fig. 3. Th»n tht fivt-viy nvltchlng valvct 
MV9 set MO that - Pc la aanaad by tha tranaducar. 
Finally, the rtading of - py la aloillarly ob- 
tained. 


The data reduction eaiploya tvo calibration 
curvet wliich were obtained froa a aingla calibration 
velocity. The underlying principle la that tha 
calibration ia independent of probe Raynolda nuabar 
Ke bated on probe tip dlaaetar. Careful calibra- 
tion experimenta reveal that thia condition exiata 
for Rep 1090 correaponding to a local velocity of 
5.4 a/a. Hence meatureaenta of auch low velocitiea 
auffer froa a necetaary calibration eiror. However, 
thia error affecta the velocity meatureaenta typi- 
cally by leaa than ^ percent for Rtp > 400, corre- 
aponding to a local velocity greater than 2 a/a. 


The data at each aeaaurcment location are re- 
duced uaing a coapu^r prograa by firat calculating 
the pitch coetficient (pj^ - Pg)/(P(^ - Py) • Froa 
these values an Interpolation technique la uaed to 
obtain the pitch angle 6 In the vertical plane froa 
the appropriate calibration characterlatic. The 
resulting value of 6 la utilized to deteralne the 
velocity coefficient pV*/( 2 (p 0 - py)] uaing the 
correapondifig calibration characterlatic. 


Valuca for V aa well aa the axial, radial, and 
swirl velocity component a u, v, and w, ahown in 
Fig. 4, are calculated from the velocity coefficient, 
pitch angle 6 , and yaw angle which la In the 
horizontal plane. The magnitude of the velocity 
•ector it given by 


V 


2 f V* , 
P [ 2 (Pc - P„) 



1/2 


fella below epproxisMtely 2.0 ■/• becauae the probe 
ia not very aenaltive to auch low dynamic preaaurea. 
In addition, the calibration ia not independent of 
probe Reynolda number. The preaaure tranaducer ia 
eatimated to have a meaaurement uncertainty of ap- 
proximately 0.25 percent for a typical velocity 
magnitude of 5 m/a and 0.75 percent for 2 m/a. 

An initial experiment waa performed to inveatl- 
gate the aaymptotic invariance of Reynolda number 
Re^. Thia waa needed to enaure that the data will 
be applicable to combuator-like flowa which operate 
at higher Re^ valuea than thoae desired for defini- 
tive flow viauallsation experiments. Velocity 
meaaurement traverses were repeated at the same con- 
ditions for increasing Re^j for each of the flow- 
fields at x/D ■ 2.5. Thia axial station was chosen 
becauae a probe location slightly downstream of a 
recirculation zone was :onaidered moat sensitive 
to Re^. The resulting Re^ values at which the 
experiments were conducted for nonswirling and 
s%rlrllng flews were 1.05 x 10* and 7.0 x 10**, re- 
spectively. 

Also an experiment to determine the extent to 
%rhlch the five-hole probe disturbs the flow was 
conducted. A dummy probe of approximately the 
same size and shape wa^ Inserted through the lo%#cr 
wall of the teat section, erectly opposite to the 
five-hole probe entering th jgh the upper wall. 
Comparison of the axial and swirl velocity profiles 
with and without the dummy probe, which ia traversed 
exactly like the measuring probe, were obtained. It 
ia clear that the presence of the dummy probe has 
negligible effect. 

3. Flowfleld Computations 



find the velocity components are obtained from thia 
maKnltude and the pitch and yaw angles. 

Proper calibration of the five-hole pressure 
probe is extremely critical to the accuracy of the 
experiment. Consequently, considerable care was 
exercised in the construction and use of the cali- 
bration equipment. The apparatus consists of a 
.small air Jet, a rotary table, a probe mounting 
bracket, and the Instrumentation system previously 
described. The rotary table has a rotary vernier 
wliich is readable within ^.5 minute. The probe 
mounting bracket Is secured to the rotary table, 
and supports the probe which rests In a cylindrical 
Steel collet. 

The motion of the rotary table orients the 
probe at the desired pitch angle, whereas the yaw 
is aerodynamical ly nulled. The probe aenalng tip 
remains at the centerline within the potential core 
of the jet and leas than one throat diameter down- 
stream of the nozzle discharge plane. The pitch 
uigle Is zeroed by very carefully aligning the probe 
shaft In a plant parallel to that of the nozzle 
throat. Differential pressures P|n * P 5 xnd * Py 
srr obtained at 5 degree Increments In 6 over the 
range -55*i6<55*. 

2 . 5_ Reliability 

Little inlormatlon ia available concerning the 
effect of turbulence on pressure probes In swirling 
flows.*’ Hrrwever, It Is presumed the the five-hole 
pitot probe Is sccursts irlthln spproxlmstely 5 per- 
cent for most of th# messurements. This vslue may 
Increase to 10 percent sc the velocity magnitude 


Equations representing conservation of mass, 
momentum, turbulent kinetic energy and its dissipa- 
tion rate are represented In the general form for 
axisynnetrlc flows vis 

i (pur^) 4 ^ (pvr^) - ^ (rr^ |J) 

• ^ s 

where 9 represents any of the dependent variables 
u, V, w, k and C. These equations are the Reynolds 
equations for axisynnetrlc flow written In diver- 
gence form. The standard two-equstlon k-€ turbu- 
lence model Is employed, whereby the exchange coef- 
ficients Fa may be specified In each equation. The 
equations differ primarily In their final source 
terms S^. The corresponding finite difference equa- 
tions are solved vis an advanced version of the 
TEACH computer code,* seml-lmpllclt llne-by-line 
method for values at points of s variable size rec- 
tangular grid, with variable under. .laxstlon. A 
complete description of the finally developed com- 
puter program Is available. 

Current messurements for mean velocities at 
the inlet are employed as boundary conditions, while 
the inlet values for k and c srt presumed. Further, 
previous svirilng flow messurements are incorporat- 
ed aa wall functions to avoid the expense of compu- 
ting within the boundary layer. Zero velocities on 
all walla are assumed, vlth symswtry conditions for 
most variables along the centerline. 
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^ ^ low Vliyllf tlon 

Rtcirculatlon tones are laportent to coabuetor 
designers baceuse aost of the burning occurs In 
these regions nnd they exhibit the highest teapera- 
turcs.** The site end location of these regions In 
the present isothemal flows are deduced froa flow 
/ visuslitation photographs of tufts, saoke, and bub- 
bles responding to the experiaental flowfield 
patterns. Resulting dividing streamline sketches 
as «#ell as selected photographs of the visualisation 
experiments are now presented and discussed. 

A.l A rtis tic Impressions of Streamline Patterns 

i Photogiaphs of each of the six flowfields re- 

sulting froa ^ • 0, 45, and 70 degrees with a • 90 
and 45 degrees have been examined in detail for each 
of the three flow visualization aethods currently 
employed. The characteristics of the overall flow- 
field are illustrated and discussed via the result- 
ing time-mean dividing streamline patterns. These 
are sketched in Figs. 5 and 6 from information ob- 
tained from the entire collection of flow visualiza- 
tion photographs. Results from the smoke-wire ex- 
periment are utilized near the inlet, whereas tuft 
and bubble data are used in approximating the size 
and shape of the recirculation zones downstream. 
Also, bubble flow patterns reveal the existence of 
a precciising vortex core, %fhich occurs downstream 
of the central region. 
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Fig. 6 Artistic impressions of dividing stream- 
lines with wall expansion angle a ■ 45* 
for swirl v.Tne ingles: (a) ^ • 0*. 

(b) 9 - 45*, and (c) ♦ - 70*. 


.0 51 1 

0 1.0 2 ^ 3J0 

U) • • s* 



The resulting streamlines for the three cases 
of the sudden expansion a ■ 90 degree geometry arc 
shown in Fig. 5. The nonswirling flow sketch in 
part (a) exhibits a large comer recirculation zone 
which is in excellent agreement with the correspond- 
ing streamlines froa the measuresMnts of Ghaturvedi.* 
Appropriate prior measurements for the present swirl 
ing flows have not been found. For the moderate 
swirl vane angle case of 0 • 45 degrees a central 
recirculation region appears in conjunction with a 
decrease in size of the comer zone. A thin pre- 
ceasing vortex core, discussed at length by Syred 
and Beer,’* is observed rear the centerline extend- 
ing froa the end of the central region to the teat 
section exit. The axial location where the vortex 
core begins fluctuates, ranging approximately froa 
x/D • 1.25 to 1.75. This vortex core is essential- 
ly a three-dimensional time -dependent phenomenon, 
which occurs as a swirling region of negligible 
axial velocity whose center winds around the test 
section centerline. A further Increase in vane 
angle to • 70 degrees results in slight enlarge- 
ment of the central xone; however, the comer bub- 
ble is essentially unaffected. For this flowfield 
the vortex core is slightly expanded in the radial 
direction. 


Fig. 5 Artistic impressions of dividing stream- The corresponding sequence of dividing stream- 

lines with wall expansion angle a - 90* lines is found in Fig. 6 for the gradual flow ex- 
for swirl vane angles: (a) 0-0*, pension case with a • 45 degrees. As with Fig. 5, 

(b) 0 • 45*, and (c) 0 • 70*. Che nonswirling flowfield exhibits excellent agree- 


a 


\ 


i 


fMHt with corresponding •eaturewonto of Chatunrodi.* 
Ho evidence of e corner tone le found froa exealne- 
tlon of Che entire photogreph collector for the In- 
temedlate vene angle caae. although the flow pat* 
tern la othervlae very alallar to the correapondlng 
flow for a • 90 degreea. The flowfleld with A • 70 
degreea producea exactly the aape effecta aa in the 
abrupt expansion caae of Fig, 5. 

4.^ Flow Vlaualltation 

Tuft vlaualltation la very laportant In that It 
auppliea an overall view of local flow direction, 
fva dlacuaacd in Section 2.3, photographs at various 
shutter speeds were obtained. Slower speeds show 
•ore of the temporal behavior* although the tufta 
ore somerimea not distinctly visible in portions of 
the flowfleld. Some of the 0 M>re noteworthy photo- 
graphs are ptesented here* encompassing a range of 
shutter speeds. Velocities in recirculation sor.es 
are often somewhat lower than in other portions of 
the flowfleld, and thus under such conditions there 
may be insufficient drag on a tuft to align it ac- 
curately with the Ic^al flow direct lot!. However* 
this is taken irto consideration in interpreting the 
photographed reaults. 

The three flowfietds for a • 90 degrece are 
characterized in Fig. 7 aa photographs of the rx- 
tcrecn tufta. Fig. 7 (a) identifies tha cornar zone 
reat tachi»»nt point for the nonswirling flowfleld. 

It la found from the length ecale as approximately 
x/D ■ 2.0* wherein tufta showing no air motion ara 
considered stagnation points. Further* an indica- 
tion of the turbulence level is depictsd as ths 
rather slow shutter speed of 1/15 second was used. 

For swirl vane angle 0 • 45 degrees. Fig. 7(b) 
is a photograph taken at 1/60 second. The swirl is 
In the counterclockv'lse direction %^cn viewed from 
downstream. Thare is no evidsnee of s comer zone 
and the central region apparently extends to x/D 
1.5. Finally, the A * ^0 degree csac shown In Fig. 
7(c) exhibits a central region extending downstream 
as far as approximately x/D • 1.85 for 1/125 second 
shutter speed. 

Local details in the nonswirling f lowf ields arc 
clearly revealed through the visualization of streak- 
lines indicated from ths generation of illuminated 
*.js(.>ka. In the swirling flow cases* strong mixing 
diffuses the ^isMke so chat streaxlines are not dls- 
clnguiahable. However* under such conditions rs- 
circulatlun zone outlines arc visible* especially 
In the region near the smoke-generstlon wire. A 
selected photograph is sxhlbitcd and discussed for 
earh of the three abrupt expansion flowfields. 

The corner recirculation bubble in the non- 
wir'lng caae with wAl expansion angle a • 90 de- 
grees is revealed in Tig. 8 (a) using 1/!K) second 
shutter speed. Also* the radial location of the 
zero velocity point within tha uppar and lower cor- 
ner bubble Is estimated to b# approximately 0.15 D 
frcmi the respective walls of the test facility. 

This afrses with the velocity messurementa presenctd 
In Section 5. 

The modsrsts swirl vane sngls flow with a 1/30 
sscond exposure is shown in Fig. 8 (b). The short- 
ened comer tone is easily identified because tha 
tdjacent flow rontaina no amok# near tha inlat. 

Thla bubbla la seen to extend to approximately x/D* 
0.45. The tuft photographa for this caaa Indicata 



(a) ♦ - 0* 



(b) A - 45* 
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(c) 9 • 70* 

Fig. 7 Flow viaualitation photographa of tufts 
in the rx-plane with wall expansion 
angle a - 90* for swirl vane anglea: 

(a) A - 0*. (b) A - A5** and (c) A - 70* 

a allghtly shortar zona ending at approximately 
x/D • 0.4. The upstream portion cf the central 
zona la alao claarly aaen* as low vslocity fluid 
carrlas a denaa mass of smoks which slowly moves 
upjtresm of c«is Inlst. further* the precesslng 
vortex core ia aeen to contain only slight smoke. 
Since the core exhibits negligible axial velocity* 
the 8SK>ke la aaaentially cairlad around it by the 
high velocity fluid outside the core. 

Figure 8 (c) is a photograph at 1/8 aecond of 
tha A * 70 degree caae wherein thla core Is not ss 
dlfetlnct. In this caae eome asMke haa diffused 
Into tha core due to a slightly longer delay before 
activating tha camera abut tar. Obaarva that both 
tha corner and central zonae near tha Inlet reveal 
that the only change from those for A * 45 degrees 
la a slightly wider central region. 
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Fig. 8 Flow visualization photographs of smoke- 
wire streaklines with wall expansion 
angle a » 90* for swirl vane angles: 

(a) i|) - 0\ (b) (|) - 45*, and (c) (|)-70*. 

Soap bubbles injected into the flow upstream of 
the test section trace pathlines clearly when illu- 
minated. In relatiyely lower turbulence Intensity 
portions of the flowfield mean flow directions can 
be obtained by ensemble averaging local tangents to 
pathlines traced out by soap bubbles. This helps 
define the flowfield geometry in terms of the out- 
line of recirculation regions. 

I 

A sample flow vdsualization photograph is pre- 
sented in Fig. 9 (a) corresponding to the zero svlrl, 
90 degree expansion angle flowfield. The pi ^tograplv 
taken with a relatively long time exposure (1/8 
second), clearly shows a great number of individual 
pathlines. Photographs of this type can be used to 
indicate regions of weakly turbulent flow such as 
that near the centerline of the flowfield which ex- 
hibits relatively straight pathlines. In addition, 
the outline of the comer recirculation region can 
be estimated from Fig. 9 (a) (and numerous additional 


Fig. 9 


Flow visualization photographs of path- 
lines indicated by Illuminated neutral- 
ly-buoyant soap bubbles for wall expan- 
sion angle a > 90** and swirl vane 
angles: (a) ^ - 0®, (b) ^ - 45®, and 

(c) 4) « 70®. 


photographs taken at the identical run condition). 
For this geometry the mean stagnation point defin- 
ing the end of the recirculation zone appears to be 
at x/D - 2.0. 

A photograph with ^ ■ 45 degrees and 1/8 second 
shutter speed is shown in Fig. 9 (b), where the prs- 
cessing vortex core is clearly seen extending from 
x/D ■ 1.5 to the exit. Its upstream extent fluctu- 
ates randomly from approximately x/D ■ 1.25 to 1.75. 
The comer bubble is observed in both the upper and 
lower portions of the flowfield, extending to ap- 
proximately x/D ■ 0.4, which agrees almost exactly 
with the smoke flow pattern for this flowfield in 
Fig. 8 (b). A photograph using 1/8 second shutter 
speed is presented in Fig. 9 (c) for the ^ ■ 70 de- 
gree flowfield where a thicker vortex core is ob- 
vious. The comer zone is faintly visible here. 





and 1C8 axial length seems to extena to about x/D ■ 
0.4. 

^ 3 Par amet r ic Effects 

Streamline plots allow recirculation zones to 
be characterized parametrically. The parametric 
effects are lllustri.ted in Figs. 10 and 11, which 
show the effects of a and ^ on the comer and cen- 
tral recirculation zone lengths. Several observa- 
tions should be noted. Firstly, zone lengths are 
only slightly affected by a, as found previously 
for the corner bubble under isothermal nonswirling 
conditions.***^ Note that, at these high Reynolds 
numbers, the corner region length is independent 
of Reynolds nu^aber. ***** Secondly, as side-wall 
angle a decreases from to 45 degrees, the zone 
lengths typically d^iease slightly and the inlet 
flow is encouraged to impinge on the confining 
walls. Thirdly, the corner recirculation length 
decreases upon increasing the swirl vane angle ^ . 
from 0 to 45 degrees, a parameter change which also 
provokes the existence of a central recirculation 
bubble and precessing vortex core. Increasing the 
vane angle to 70 degrees only slightly enlarges the 
central zone and the vortex core. 
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Fig. 10 Experimental effect of wall expansion 
angle a and swirl vane angle ^ on cor- 
ner recirculation zone length. 

5^ V elocity Measurements 

5. 1 Sudden Expansion a ■ 90 Degrees 

Figures 12, 13, and 14 show the axial and 
iwirl time-mean velocity profiles, for swirl vane 
angles C* * 0, 45 and 70 deg. with a sudden expan- 
sion [side-wall angle a ■ 90 deg.]. Measurements 
ire obtained with the five-hole pitot probe in the 
manner discussed in Section 2.4. Note that in 
these and subsequent figures, different scales are 
used for the normalized axial and swirl velocities. 
Radial velocities are consistently much small than 
Che axial and swirl components, and they are indi- 
cated in the dividing streamline patterns deduced 
from flow visualization experiments depicted in 
Figs. 3 and 6. They need not be presented here. 
Figure 12 shews the comer recirculation region for 
the nonswirling flow on the verge of reattachment 
near x/0 ■ 2.0. The corresponding visualization 



Fig. 11 Experimental effect of wall expansion 
angle a and swirl vane angle 0 on the 
central toroidal recirculation zone 
length. 

results shown in Fig. 5 (a) give this location as 
x/D ’■2.1 and the corresponding measurements of 
Chaturvedi* give a value of x/D ■ 2.3. 

In the absence of pertinent measurements, pre- 
vious oredictions of gas turbine combustor flow- 
fields^*******® employed an approximation for inlet 
velocity boundary conditions. It has generally 
been assumed that v ■ 0 and that both u and w exhi- 
bit a flat velocity profiles. However, the present 
measurements indicate this is an unrealistic esti- 
mate, with sharply peaked u and w profiles, as 
shown in Figs. 13 and 14. The difference results 
from the use of a swirler with ten flat blades with 
pitch/chord ratio of unity, imperfect blade effi- 
ciency, the existence of a hub, and the fact that 
the downstream edge of swirl vanes of the test fa- 
cility is actually located approximately 4 cm up- 
stream of the flow expansion comer, where x/D • 0. 
This swirler location allows the central recircula- 
tion zone to begin upstream of x/D ■ 0, thereby 
changing the velocity profiles there. 

Figure 13 also reveals that a 45 degree swirl 
vane angle produces a maximum swirl flow angle 
{tan*' (w/u)} of 30 degrees at the test chamber in- 
let. The maximum swirl flow angle for the 70 de- 
gree vane angle case shown in Fig. 14 is only 34 
degrees. Thus there is only a slight increase of 
swirl flow angle, although these two flows are dif- 
ferent in that the inlet profiles are considerably 
more sharply peaked for the latter case. Figures 
13 and 14 show zero u and w velocity values near 
the axis at the inlet, but actually the probe was 
insensitive to the very low velocities there. This 
is consistent with earlier flow visualization re- 
sults that the central bubble extends upstream of 
the inlet. The measurements shown in Fig. 13 (a) 
provide no evidence regarding the existence of a 
corner zone. This is expected because flow visual- 
ization reveals that the region only extends to 
x/D • 0.4. However, there is clear evidence of a 
rather large central zone whose length is similar 
to that shown in Fig. 3 (b). Although the axial 
velocity profiles are beginning to flatten in the 




Fig. 13 Measured velocity profiles for vail expansion 
angle a • 90* and swirl vane angle ^ ■ 45*. 
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Fig. 14 Measured velocity profiles for wall expansion 
angle a • 90* and swirl vane angle ^ • 70*. 









downutrean direction, they retain a zero velocity 
value on the axis, which la conalatcnt with soap 
bubble flow patterns as seen in Fig. 9 (b). The 
early erratic behavior shown by the swirl velocity 
profile In Fig. 13 (b) quickly transforms, exhibi- 
ting a solid-body-rotation core with a rather flat 
profile outside this region. 

Figure 14 also reveals the discrepancy regard- 
ing inlet velocity profiles. The large central re- 
circulation region causes the downstream flow to 
be accelerated near the top wall. The correspond- 
ing swirl velocity component at the x/D *^0.5 axial 
station typically shows local minima where the 
mean axial velocity Ir zero. The accuracy of the 
velocity measurements at these locations Is poor 
as discussed in Section 2.5, so that any physical 
interpretation here is suspect. The early erratic 
behavior found In Che swirl velocity profile at 
x/D * 0.5 quickly develops into a shape similar to 
that seen in Fig. 13 (b). The precessing vortex 
core motion discussed in Section 4.2 results in 
poor measurement repeatability which promotes the 
irregular behavior within this core region. 

5 . 2 Gra dual Expansion a ■ 45 Degrees 

Figures 15, 16, and 17 exhibit velocities for 
th€ same sequence of flowflelds with side-wall 
angle a * 45 degrees. The inlet profiles were not 
measured in this geometry because the presence of 
the expansion block interferes with probe position- 
ing. Effects of swirl vane angle 0 on velocities, 
similar to chose found for the sudden expansion 
cases, are found in the flowfield sequence for this 
test section geometry. The major difference Is 
that the sloping wall encourages the inlet flow to 
accelerate near the top wall. Also, it ends to 
shorten or obliterate the corner recirculation 
region. 

^ 3 Comparison of Measurements and Predictions 

Flowfield computations have been made^"** for a 
variety of side-wall and swirl vane angles, using 
the technique discussed briefly in Section 3 and at 
length elsewhere.^** These and other earlier pre- 
dictions have generally idealized the inlet flow as 
a plug-flow axial velocity profile with a flat or 
solid body rotation swirl velocity profile. Though 
this may be adequate for those applications, it is 
clearly inadequate for the present study, in which 
the inlet flow conditions are achieved via a ten- 
blade svirler with pitch/chord ratio of unity. As 
documented in Section 5.1 and 5.2, the present in- 
let profiles are highly nonunlfonn. 

As an example cf the type of predictions re- 
sulting from application of the STARPIC code*** to 
the flowflelds studied experimentally. Fig. 18 
shows predicted and measured velocity profiles for 
the case of the sudden expansion flow (a ■ 90 deg.] 
with swirl vane angle 0 ■ 45 deg. Inlet profiles 
to this prediction are taken as the measured pro- 
files at x/D • 0. The prediction shows that the 
large peak of Inlet axial velocity at x/D - 0.5 Is 
located at r/D ■ 0.25 rather than 0.40 as measured. 
Also, a much smaller central zone is predicted. 
However, no effect of this discrepancy is seen be- 
yond x/D ■ 1.5. The measured inlet swirl velocity 
peak troadens very quickly by x/D ■ 0.5. Also, a 
solid-body-rotation core is exhibited there. Com- 
parison with measurements shows considerably larger 
disagreement at 'x/D • 0.4 than downstream. This 
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discrepancy near cne liuet may result partially from 
backflow at the inlet boundary [which could not be 
measured becauat of transient effects], the lack of 
probe sensitivity to small velocities near the cen- 
terline, and/or the use of the standard two-equation 
k-C turbulence model. In this connection, it may be 
nctad that a turbulence model whose basis and para- 
metara are adequate for simple flow situations is 
not adequats to handle the more complicated swirling 
recirculaing flow situation. Nevertheless, good 
predictions are available** for the use of a less 
difficult test case (co-swirl and counter-swirl flow 
in a pipe using the data of Ref. 31]. This turbu- 
1 l..C€ simulation problem is complex flowflelds is 
clearly an area of current research interest. 


6? Closure 


A major outcome of the current study is the ex- 
perimental characterization of corner and central 
recirculation zones in six basic flowfield configu- 
rations of an axisynmetric expansion with side>wall 
angle oi * 90 and 45 deg. and swirl vane angle « 0, 
45 and 70 deg. The size and shape of the recircula- 
tion bubbles for each flowfield is illustrated via 
an artistic impression deduced from a collection of 
flow visualization photographs of tufts, smoke, and 
neutrally-buoyant soap bubbles responding to the 
flow. Increasing swirl vane angle ^ from 0 to 45 
degrees produces a shortened corner region and the 
appearance of a central bubble typically extending 
downstream to approximately x/D ■ 1.7, after which 
a precessing vortex core exists near the axis reach- 
ing CO the exit of the test section. A further in- 
crease in ^ to 70 degrees enlarges the central zone 
and vortex core with negligible effect on the corner 
region in those flowflelds where it occurs. The 
effect of side-wall angle a on the nonswirling flows 
Is negligible. However, a decrease from 90 to 45 
degrees apparently eliminares the comer bubble in 
the swirling flow cases investigated. This decrease 
in a also causes the inlet flow to impinge more 
severely on the top wall, where larger axial velo- 
cities occur. 

A more detailed experiment consists of the 
measurement of time-mean velocity components in the 
axial, radial, and azimuthal directions using a 
five-hole pitot probe. These measurements generally 
agree with the flow visualization results and pro- 
vide a more complete understanding of each flowfield 
At the inlet, the axial and swirl velocity profiles 
exhibit maximum values at approximately r/D ■ 0.2 
In a sharply-peaked annular fashion. This nonuni- 
formity arises for several reasons: the use of a 

ten blade swirler with pitch/chord ratio of one, 
combustor flowfield predictions are unrealistic. 

This arises from blade inefficiency, the presence 
of a hub, and the fact that the swirl vane exit 
station Is typically located upstream of the expan- 
sion station. This allows the central recirculation 
zone to begin upstream of the expansion where 
x/D - 0. 

An advanced computer code equipped with stan- 
dard two-equation k-C turbulence model is used to 
predict corresponding flow situations and to com- 
pare results with the experimental data. Although 
excellent agreement is found for the nonswlrllng 
flow, poor agreement in detail occurs for the swirl- 
ing flow cases, especially near the inlet. This 
discrepancy may be attributed to the lack of realism 
in the k-C turbulence model and/or to Inaccurate 
specification of time-mean velocity and turbulence 
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Fig. 15 Measured velocity profiles for wall expansion 
angle a ■ 45* and swirl vane angle d * 0*. 
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Fig. 16 Measured velocity profiles for wall expansion 
angle a ■ 45* and swirl vane angle 0 ■ 45*. 
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Fig. 17 Measured velocity profiles for wall expansion 
angle a ■ 45* and swirl vane angle ^ - 70*. 








Fig. 18 Predicted and measured velocity profiles for a specific 
test case with wall expansion angle a ■ 90* and swirl 
vane angle ^ • 45*, 

( predictions, A-A-A measurements] 


energy distributions at the inlet. Clearly, the 
time-mean velocity measurements constitute a seri- 
ous Iv needed data base for the validation of com- 
puter prediction codes and the development of tur- 
bulco e models for their simulation. 
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CHAPTER I 


INTRODUCTION 

1.1 Combustor Design and Development 

The combustor of the gas turbine engine illustrated in Figure 1 con- 
tains high intensity combustion and, as far as possible, must burn fuel 
completely, cause little pressure drop, produce gases of nearly uniform 
temperature, occupy small volume, and maintain stable combustion over a \ 

wide range of operating conditions (1). Efforts are continually being | 

expended to produce efficient clean combustion. The designer has a for- j 

I 

midable problem in aerothermochemi st ry , and more accurate analysis pro- 
cedures can provide a route which leads to the accomplishment of design 
objectives more quickly and less expensively than current practice per- 
mits (2, 3). 

In design situations, the engineer has to seek an optimum path be- 
tween alternatives of, for example, efficiency and pollution. The gener- 
al aim of most research investigation is to provide information which is 
useful to designers by "characterizing" or "modeling" certain features 
of the phenomenon in question. Investigations may be theoretical or ex- 
perimental; the two approaches are complementary. Until recently, design- 
ers have relied heavily on the experimental approach, but traditional 
design methods (empirical formulas and experimentation) are now being sup- 
plemented by analytical methods (numerical solution oF the appropriate 
governing partial differential equations). The mathematical modeling 
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approach is now finding favor and is be ng used to suppleinent existing 
design procedures. Current combustor design and development problems, 
the needs of the combustion engineer in practice, and proposed research 
tasks which will assist in the attainment of design objectives are becom 
ing clear. Improvements and new developments (both experimental and the 
oretical) can and should be made, theoretical modeling being aided by 
specific carefully chosen experiments (k) . 

1.2 Practical Needs and the Present Research 


Many details of combustor phenomena ren^ain poorly understood at the 
macroscopic as well as the microscopic level. Knowledge of these will 
undoubtedly reduce combustor development time and cost substantial ly, and 
lead to increases in performance. Some of these research needs are: 

1. Physical processes--turbulence, radiation, combustion, and mul- 
t i phase effects . 

2. Computer programs--0- , 1-, 2-, and 3“bimensional approaches in 
steady state and transient cases. 

3. Unresolved prob lems--ef feet of swirl, recirculation, and wall 
proximity on turbulence; turbulence-react ion interaction; and multiphase 
simulation. 

The present research work is concerned with complementary experimen- 
tal and theoretical studies on a specific subproblem of this general prob- 
lem and is described in a recent paper (5). The specific problem bein^ 
investigated experimentally is concerned with steady turbulent flow in 
axisymmetric geometries, under low speed and nonreacting cond i t i ons--a 
study area highlighted recently as a fundamental research requirement in 

1 

combustion modeling (2, 6, 7). The particular problem is concerned with I 
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turbulent flow of a given turbulence distribution in a round pipe enter- 
ing an expansion into another round pipe» as illustrated in Figure 2. 

The in-coming flow may possess a swirl component of velocity via passage 
through swirl vanes at angle (equal approximately to tan ' ) * 

and the side-wall may slope at an angle a, to the main flow direction. 

The resulting flowfield domain may possess a central toroidal recircula- 
tion zone in the middle of the region on the axis, in addition to the 
possibility of a corner reci rculat ion zone near the upper corner provoked 
by the rather sudden enlargement of the cross-sectional area. Of vital 
concern is the characterization of flows of this type in terms of the 

effects of side-wall angle a, swirl vane angle , turbulence intensity k. 

I n 

of the inlet stream, and expansion ratio D/d on the resulting flowfield 
in terms of its time-mean and turbulence quantities. Such problems have 
received little attention, yet there is a definite need for work in this 
area even under nonreacting flow conditions. 

1.3 Theoretical Investigation 


1 


I 


A mathematical solution of the flowfield of interest should provide 
results more economically, quickly, and correctly than possible by other 
means (for example, experiments on real-life systems or models). In order 
to achieve this, the model should simulate the flow in all its important 
respects (geomet ry , boundary cond i t ions , phys i ca 1 properties of gases, tur- 
bulence, etc.) and provide a means whereby the governing equations may be 
solved. Mathematical models of steadily increasing realism and refine- 
ment are now being developed, both in the dimensionality of the model 
(together with the computational procedures) and in problems associated 
with the simulation of the physical processes occurring. Clearly there 





are two areas of difficulty: the simulation and the solution. Recent 


publications discuss these difficulties in the context of combustor flow- 


field predictions (4, 5, 8, 9t 10). 


An advanced prediction computer program, which solves the governing 


partial differential equations in finite difference form, has been devel* 


oped. A detailed discussion of this program, including a user-oriented 


section, has been prepared recently (11). Therefore, extensive discus- 


sion of this program in Chapter I I I of this document is not necessary and 


only highlights of the technique are included. 


1.^ Experimental Investigation 


Prediction work in the form of computing the flowfield is being com- 


plemented by an associated experimental study which is used to evaluate 


the final predictive capability. The experiments, under way at Oklahoma 


State University, are concerned with measuring the effects of swirl and 


side-wall angle on streamlines as well as mean flow and turbulence para- 


meters in nonreacting flow. Use is made of the confined jet facility 


shown schematically in Figure 3. 


The facility has an axial flow fan whose speed can be changed by 


means of a varidrive mechanism. Numerous fine screens and straws produce 


flow in the settling chamber of relatively low turbulence intensity. The 


contraction section leading to the test section has been designed to pro- 


duce a minimufr. adverse pressure gradient on the boundary layer and thus 


avoid unsteady problems associated with local separation regions. The 


sudden expansion consists of a 15 cm diameter circular jet nozzle, exit- 


ing abruptly into a 30 cm diameter test section, as shown in Figure 3. 




r 


I 


I. 


mm 
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The substantial size of this test model provides excellent probe resolu- 
tiori for five-hole pitot meacurcmenis . The test section is constructed 
of plexiglass to facilitate flow visualization. The side-wall angle 
and swirl vane angle ^ are variable. The side-wall angle is set by in- 
serting one of three blocks with a side-wall angle a of 90, 70, or ^5 de- 
grees. The swirl vane angle is continuously variable with discrete val- 
ues of ♦ ■ 0, ^5» and 70 degrees being chosen for the present study. 

The experiments are concerned with nonreacting flow, and the following 
two aspects of the study are addressed in the present thesis: 

1. Photography of neutral ly-buoyant helium-filled soap bubbles, 
tufts, and smoke, so as to characterize the time-mean streamlines, recir- 
culation zones, and regions of highly turbulent flow. 

2. Five-hole pitot probe pressure measurements to determine time- 
mean velocities u, v, and w. 

Complete details of these experimental studies appear in Chapter IV. 

1.5 Swirl Flow Concepts 

Most practical combustion equipment includes the phenomenon of swirl 
and it is convenient to recall here some of the basic concepts of its gen- 
eration and characterizat ion (12). Swirling flows result from the appli- 
cation of a spiraling motion, with a swirl velocity component (also known 
as a tangential or azimuthal velocity component) being imparted to the 
flow via the use of swirl vanes, in an axial-plus-tangential entry swirl 
generator or by direct tangential entry into the combustion chamber. Ex- 
perimental studies show that swirl has large scale effects or flowficlds; 
jet growth, entrainment, and decay (for inert Jets) and flame size, shape. 
Stability, and combustion intensity (for reacting ^lows) are affected by 


6 


the degree of swirl imparted to the flow. This degree of swirl usually 
is characterized by the swirl number S, which is a nondimensional number 
representing axial flux of swirl momentum divided by axial flux of axial 
momentum times equivalent nozzle radius. That is 


S 


G d/2 

X 


( 1 . 1 ) 


where 


J (puw + pu'w' ) 


,^r 


is the axial flux of swirl momentum, including the x6*directIon 
turbulent shear stress term; 

r® 2 2 

G ■ 1 (pu '♦•pu* -f (p - p )) r dr 
X o * 


is the axial flux of axial momentum, including the x-direction 
turbulent normal stress term and a pressure term; d/2 Is the nozzle 
radius; and u,v,w are velocity components in (x,r,G) cylindrical 
polar coordinate directions. 

U . free j« in 5t.9n.nt surroundings, .nd G, .re const.nts, th.t is. 
invariants of the Jet, as may be deduced from radial integration of the 
appropriate governing Reynolds equation (13). 

It is especially convenient to be able to relate the angle of a 
swirl vane pack to the swirl number produced. In this context, for com- 
parison purposes it is worth noting that swirl vane angle ^ and swirl 
number S are related approximately by 


S 


2 

3 


I - (d^/d) 
1 - (d^/d) 



tanp 


( 1 . 2 ) 












nssk 


7 


where d and d. are nozzle and vane pack hub diameters, respectively. This 
h 

relationship follows from assumptions of plug flow axial velocity u^ in 
the annular region, and very thin vanes at constant angle 4> to the main 
direction so imparting a constant swirl velocity = u^ tan4> to the flow. 
From the 6^ and definitions, integration over the range (= to 

R (= d/2) gives 



(1.3a) 

(1.3b) 


from which the quoted relationship follows. In the case of a hubless 

swirler (or one for which d,/d is very sj^'dll) the expression simplifies 

h 

to 

S = y tani{) ( 1 .^) 

so that vane angles of 15, 30 , 45, 60, 70, and 80 degrees, for example, 
correspond to S values of approximately 0.2, 0.4, 0.7, 1-2, 2.0, and 4.0. 
Here 100 percent efficiency also is assumed for the swirl vanes, which 
deteriorates as the vane angle ii creases. In fact, flat blades are not 
perfect, rnd indeed the assumption of flat axial and swirl velocity radial 
profiles immediately after tne swirler is not correct, as may be seen from 
experimental measurements in Chapter VI of this tnesis and elsewhere (’2, 
13, 14). 


1.6 The Present Contribution 


The basic objectives of the present research Include the development 
of a computer code for the computation of swirling, axisymmet n c , conf i ned 
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jet flows using state-of-the-art computational procedures and the k-e 
turbulence model. Further objectives entail mean flowfield measurements 
and flow visualization to characterize this particular flow and to evalu- 
ate the developed prediction procedure. The computations and experiments 
have been performed for six basic confined jet flowfields, all of which 
conform to the schematic diagram of Figure 2 and have an expansion ratio 
D/d =■ 2. The parameter combinations which define these flowfields are 
side-wall angle a = 90 and ^5 degrees and swirl vane angle = 0,^5, and 
70 degrees. The computational techniques are discussed in Chapter III. 

The experimental program consisted of performing mean velocity mea- 
surements wi th a f i ve-hole pi tot probe and flow visualization of neutral ly- 
buoyant soap bubbles, tufts, and smoke. Details of the experimental 
methods are presented in Chapter IV. The pitot probe measurements re- 
sulted in maps of axial, radial, and swirl velocity components while the 
flow visualization provided streamline patterns and recirculation zone 
geometries. Results of the overall flowfield characteristics highlighted 
by the reci rculat icn region geometries, both predicted and visualized, 
are presented in Chapter V. Chapter VI contains results of predictions 
and measurements of time-mean velocity components. Finally, Chapter VII 
summarizes the conclusions of the present investigation and gives recom- 
mendations for further work. 


CHAPTER II 


\ 



REVIEW OF PREVIOUS AXISYMMETRIC 
RECIRCULATING FLOW STUDIES 

Previous experimental and theoretical work in axisymmet ric recircu- 
lating flows provides an important background for this investigation. 

Information regarding both swirling and nonswirling, nonreacting, turbu- 
lent flows has been extensively surveyed. Noteworthy results are includ- 
ed in the following summary where the emphasis is on nonreacting turbu- 
lent swirling flows . 

2.1 Experimental Work in Nonswirling Flows 

The effect of expansion angle a was investigated by Chaturvedi (15) 
who measured mean and turbulent flow quantities behind expansion half- 
angles of 15, 30, 45, and 90 degrees for air flow in a pipe. The expan- 
sion diameter ratio was 2.0. Measurements of velocity head in regions of 
high turbulence intensity and where the direction of the velocity vector 
was unknown were made with a 2.5 mm diameter pitot tube. Mean velocity ^ 

was also measured with a constant temperature hot-wire anemometer using a 
single wire. A cross-wire was used to measure all the Reynolds stresses. 

Krall and Sparrow (16) utilized an orifice in an electrically heated 
tube to create flow separation. It was presumed that the peak Nusselt 
numoer occurred at the reat taenmont point. Their results agreed with 
those of Phaneuf and Net.^er (17) in that reattachment length wa^ unaffected 
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by Reynolds number in the turbulent range. Also, they found that peak 
Nusselt number became spread out into a broad zone for strong separations 
caused by large step heights. This transition from a point to a broad 
zone occurred between h/D = 0.16 and h/D = 0.25, where D is the pipe dia- 
meter and h is the step height. Also, it was found that increasing the 
step height moved the peak point slightly downstream. 

Back and Roschke (18, 19) investigated the reattachment length be- 
hind a pipe step by visual observation in water flows. The diameter ex- 
pansion ratio was 2.6 and a conical contraction section just upstream 
reduced the boundary layer thickness. Small 0.?6 mm diameter holes were 
located along the larger diameter tube one step height apart. Reattach- 
ment locations were determined by slowly metering dye through these holes 
and observing whether it moved upstream or downstream. 

Wall static pressure was measured by Phaneuf and Netzer (17) as well 
as mean axial velocity behind abrupt pipe steps wherein the step heights 
were 0.303 and 0.195 outlet diameters. They found that the reattachment 
zone spreads out with both increasing Reynolds number (for a given step 
height) and increasing step height provided the Reynolds number is in the 
turbulence transition range. It was further concluded that the point of 
maximum heat transfer does not coincide with the zone of reattachment, 
and furthermore, it is dependent upon the inlet mass flux. 

Ha Minh and Chassaing (20) studied the rest ructurat ion of a turbu- 
lent pipe flow past three types of geometric perturbations, one of which 
was the abrupt pipe expansion. Agreement between hot-wire and pitot 
probe measurements of time-mean axial velocity was exhibited. Owing to 
probe resolution cons iderat ions , only normal or inclined single wire 
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probes were used. All Reynolds stresses were measured, and the rotated 
inclined hot-wire technique was utilized. 

Both a one-color and two-color laser-Doppler system were utilized by 
Moon and Rudinger (21) to measure mean axial velocities behind an abrupt 
pipe step. The diameter expansion ratio D/d was I.A 3 . Velocity reverses 
in the recirculation zone were easily detectable through use of a Bragg 
cell to shift the frequency of one probe beam relative to the other. 

An e .perinoental technique was developed by Johnson and Bennett (22) 
for measuring mean and fluctuating velocity and concentration distribu- 
tions along with their cross correlations. Measurements were obtained 
using a laser-Doppler velocimeter to measure velocities, and a laser- 
induced fluorescence technique to measure the concentration of fluorescent 
tracer particles. Such measurements were displayed for coaxial jets dis- 
charging into a larger axisymmetric duct with D/d = 2.07. 

2.2 Experimental Work in Swirling Flows 

Both an aircraft and an industrial type combustor were investigated 
by Hiett and Powell (23) under reacting and non'^eacting conditions. A 
water-cooled three-dimensional pitot probe was developed. Mean velocity 
profiles, recirculation boundaries, and temperature contours were present- 
ed. Also, cotton tufts were employed to visualize the nonreacting flow. 

Extensive studies of free and enclosed swirling jets from vane swirl- 
ers were undertaken by Mathur and MacCallum (24). In the enclosed case, 
the expansion ratio was approximately D/d = 2.5 into a square-sect ioned 
test section. An evaluation of swirler efficiency was conducted for both 
annular (with hub) and huoless swirlers. Extensive presentation is given 

^ _ j 


of velocity profiles and reci rculat ion zone sizes as functions of the de- 


gree of inlet swi rl . 


Laser velocimeter measurements were made by Owen (25) in the initial 


mixing region of a confined turbulent diffusion flame burner. The experi* 


mental facility consisted of an axisymmetric combustor in which a central 


gaseous fuel stream mixed with a high speed coaxial annular air stream. 


The expansion diameter ratio was approximately 1.3 and swirl numbers of 


0, 0.3, and 0.6 were chosen to contrast swirling with nonswirling combust- 


ing flowfields. Measurements were presented for axial and tangential mean 


velocity profiles along with the rms and probability density distributions 


of velocity fluctuations. 


Three mean velocity components and normal stresses were measured with 


a laser anemometer by Baker et al. (26) in an experimental facility of a 


small scale axisymmetric furnace. The test chamber has an expansion dia- 


meter ratio of 3-33 and contains a constricted outlet. Measurements were 


obtained in both an isothermal air flow and a combusting mixture of natu- 


ral gas and air, with exit swirl numbers of 0 and 0.52. The boundary and 


field measurements of three components of time-mean velocity and corre- 


sponding normal stresses were designed to be of good use in the evaluation 


of turbulent flow prediction procedures. It was clear from the results 


that regions o? recirculation under reacting flow conditions differed sub- 


stantially from those under nonreacting conditions. Surprisingly, they 


were found to be in general larger under reacting flow conditions. The 


turbulence was found to be far from isotropic over most of the flow field. 


Beltaqui and MacCallum (27» 28) studiea premixed vane swirled flames 




in furnaces with expansion ratios D/d = 2.5 ana 5.0. where D and d are 
chamber and nozzle diameters, respectively. They found that a given vane 
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swirler can produce different velocity profiles as D is varied, and that 
the size and shape of the central reel rcul at ion zone are primarily func- 
tions of D rather than d, A modified swirl number using D rather than d 
was suggested as being more characteristic. 

Mean velocity measurements were obtained by Vu and Gouldin (29) for 
coaxial jet mixing in a pipe under coswirl and counterswirl conditions. 
Both the central and annular flows had a swirl number of approximately 
0 . 5 . In the counterswirl case, a recirculation zone formed on the axis 
with very steep velocity gradients. In the coswirl case, no recircula- 
tion was observed. 

Syred et al. (30) concentrated on turbulence measurements in strong- 
ly swirling flows via a single-wire six-orientation hot-wire technique. 
Local turbulence intensities were found to be extremely high in and near 
the central reel rculat ion zone. Measurements of all six turbulent stress 
components showed strong variations of absolute turbulent kinetic energy 
levels and strong nonisotropy of the stresses and associated turbulent 
viscos i ty . 

Further work of less importance to the present study is readily 
available (31 . 32, 33, 3^) . 

2.3 Theoretical Work in Nonswirling Flows 

Laminar flow over a pipe step has been computed for the diameter ex- 
pansion ratio 2.0 by Macagno and Hung (35). Computations for Reynolds 
numbers up to 200 were undertaken using stream function and vorticity 
variables. Flow visualization results were in good agreement with pre- 


dicted values. 
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Runchal et al. (36) also predicted the flow downstream of a pipe 1 

step, solving transport equations for stream function, vorticity, enthal- I 

py, and turbulence energy. The turbulence length scale was empirically | 

determined to ensure that the reattachment point was accurately predicted. j 

The effects of the outlet boundary condition are significant for only a 
short distance upstream, for example, one pipe diameter. The Stanton num- 
ber was correctly predicted to vary as the -0.38 power of the Reynolds | 

number. j 

As previously mentioned, Ha Mi nh and Chassaing (20) investigated the j 

flow behind a pipe step. They obtained good agreement between measured I 

and predicted values consir'e'* the extent of the flow perturbation. 

The k-E turbulence mode! wjs emplcye«j 'or most of the computations, al- 
though the five-equation Reynolds stress transport model predicted the 
centerline k*.uetic energy more jccurately. | 

Other predictive studies on nonswirling flow are discussed in sec- | , 

ticn 2.^, v^'Siere the emphe:.i » is on swirling flows and in which the non- l 

i 

v^irliiiG cast af^^ears as a degenerate case. Selected theoretical work on | 

t .Duli'rcw nx)deling is also useful (37, 38, 39, ^0). | 

2.^ Theoretical Work in Swirling Flows 

Roberts (k\) investigated similar axi symmetric, swirling flowfields 
by developing a finite difference code incorporating the k-W turbulence 
model while using stream function and vorticity variables. Poor perform- 
ance of the k-W model led to the crude approach of specifying constant 
effective viscosities where i 3 1 “ ^0- Also, static wall 

pressures were measured and visual izet ion information recorded for a 
water nxjdel experiment which gave general agreement with predictions. 

k : 




^5 

Khalil et a1. (42) predicted the flowfield of several two-dimensional 
model furnaces for which various measured quantities are available using a 
similar finite difference algorithm. As detailed knowledge of inlet boun- 
dary conditions are rarely known, a parametric study was undertaken, and 
results indicated that the inlet velocity profile has negligible effect. 
Also, provided that total enthalpy of the inlet flow is held constant, 
changes in the radial temperature profile exhibited a very small effect. 
However, it was found that inlet distributions of k and t exert an appre- 
ciable influence, and of the three inlet k profiles tested for the flow- 

2 

field measured by Baker ec al. (26), values for k/u of 0.009 for near 

o 

wall nodes and 0.003 for central nodes gave closest agreement for axial 
velocity. Nonswirling flows are more sensitive to inlet profiles and 
thus the above numerical tests were conducted without swirl. Also, no 
significant differences could be observed for various outlet boundary con- 
ditions tested for axial velocity, enthalpy, mixture fraction, and concen- 
tration fluctuation. 

An axisymmet ric, combusting, swirling flov-/ computer code was used by 
Hutchinson et al. (A 3 ) to compare predicted and measured quantities in 
the furnace test facility studied by Baker et al. (26) except at higher 
flow rates. The furnace consisted of a coaxial burner with swirling annu- 
lar air flow and a constricted exit. Measurements of mean velocity and 
normal stress in both axial and tangential directions were generally pre- 
dicted along with mean temperature and wall heat rlux. Similar calcula- 
tions using the same physical information but with a different numerical 
scherne showed that differences of up to 15 percent can be introduced by 
the numerical arrangement, whereas the rms difference^ were less than 5 
percent. Discrepancies in axial velocity ranged up to AO percent of the 
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maximum in the appropriate plane. The largest discrepancies were in the 
burner wake and it was suggested that the two-equation turbulence model 
is inadequate there. 

A finite difference code was developed by Li 1 ley (AA) for inert and 
reacting recirculating flows with strong swirl. The effect of swirl num- 
ber on mean axial and swirl velocities was presented along with predicted 
streamlines in an axi symmetric combustor. Also, the effect of swirl vane 
angle on central recirculation zone length and flame length was exhibited. 

Scrag-Eldin ana Spalding (A5) predicted the nonaxisymmetric, swirl- 
ing, combusting flow in a can combustor with air entering through a A5* 
degree vane angle swirler in an annular tube enclosing a central jet of 
fuel. Temperature contours were measured with a bare 13 percent Rh-Pt 
versus Pt thermocouple wire of 3.8 x 10 ^ mm diameter. Predicted results 
gave the same trend as measured ones, but quantitative agreement was not 
always acceptable. 

Swirling, combusting flow in axisymmetric prechamber dump combustors 
was predicted by Novick et al. (A6) using a finite difference procedure. 
Application of the coniputer code to specific combustor designs is stress- 
ed, and consideration is given to the relative effect of certain impor- 
tant design para.7>eters : main chamber to prechamber expansion racio, 

swirler dimensions, double concentric reverse swirler, constricted out- 
let, and fuel concentration profile. The standard base configuration con- 
sists of a 70-degree swirl vane angle inlet flow (stoichiometric mixture) 
which proceeds over a central hub, enters a prechamber, followed by a sud- 
den expansion. Computed streamlines and temperature contours for this 


case are shown. 
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Details of the above model Including physical assumptions* governing 
equations* finite difference equations* and solution procedure are found 
in a recent paper by Novick et al. (^7). The predictive capability is 
demonstrated for a flowfield v/ith coaxial jets entering a sudden expan- 
sion. Computed streamline patterns are presented for various swirl 
strengths with and without combustion in a prechamber dump combustor. 

Other advances to prediction procedures have been made elsewhere and 
applied to other two-dimensional axisymmetric problems involving swirl 
flows in expansions of interest to the combustor designer. Habib and 
Whitelaw (^8) compared hot-wire and laser measurements with predictions 
in confined coaxial jets with and without swirl. 

Predictions of the swirl mixing pipe flow data of Vu and Gouldin 
(29) were made by Srinivasan and Mongia (^9) for both coswirl and counter- 
swirl cases. The measured inlet profiles were utilized as inlet boundary 
conditions with both the original k-e turbulence model and the version ex- 
tended to include Richardson number effects. The extended model gave bet- 
ter 3greement with measurements* but the recirculation zone was somev^hat 
elongated. 

Sturgess and colleagues have predicted a variety of flow types with 
and without combustion as part of a continuing validation study of turbu- 
lence and reaction models for aircraft gas turbine combustor application 
(50, 51* 52). Associated sudden expansion experiments are yielding much 
needed data for this evaluation (22). 

Lllley et al. (9) discussed the predicted effects of swirl vane angle 
(f and inlet side-wall expansion angle a on axisymmetric idealized gas tur- 
bine combustor flowfields. Also* the predicted effects of a and p on both 
the corner and central reci rculat i on zone lengths were presenteo, under 



uoth reacting and nonreacting flow conditions, so extending earlier iso- 


thermal studies, 


Other papers presented at a recent conference address turbulent re- 


acting flows in practical combustion systems (53)- The computation of 


flows with complex chemical kinetics problems, complex fluid mechanic 


problems, and complex simulation problems are all considered. 


Fully three-dimensional flowfields are not specifically addressed in 


the present thesis. However, by way of introduction, interested readers 


are referred to the work of Mongia and Reynolds (5^), Serag-EIdin and 


Spalding (45), and Swithenbank et al. (55). Mongia and Reynolds (54) 


developed and applied an advanced fully three-dimensional prediction tech- 


nique to compute the reacting flowfield as part of an analytical design 


procedure for gas turbine combustion systems. Complexities such as the 


eddy-breakup chemical reaction model, and fuel spray droplet phenomena 


were included. Good qualitative agreement was found to be quite reward- 


ing in developing a number of relatively complex combustion systems. 


Gaseous phase diffusion flames in can-type combustors have been predicted 


(45). Patterns of velocity and temperature show that qualitative trends 


may be simulated. A three-dimensional two-phase mathematical model of a 


Lycoming combustor has been ama 1 garriated by Swithenbank et al. (55). Some 


30 partial differentia) equations were solved and encouraging results 


were obtained. 






CHAPTER 111 


COHPUTATIONAL TECHNIQUES 

The economical design and operation of practical fluid flow equip- 
ment can be greatly facilitated by the availability of prior predictions 
of the flowficld. These may be obtained by use of a mathematical model 
i ncorporat ing a numerical finite difference solution procedure, and sev- 
eral textbooks discuss the approach at length (56, 57, 58). Several re- 
cent review papers also provi de useful details ( 1 ^ , 37, ^2 , , ^5 1 ^7 • 52 , 

59). These theoretical studies complement associated experimental stud- 
ies, and their improvement and use can significantly increase understand- 
ing and reduce the cime and cost of development work. In the context of 
the present study, the main goals include improved flowfield simulation 
and turbulence model development, and computer code development plays a 
significant role in this endeavor. 

The computer program STARPIC (acronym for girling ^urbulent ^xisyn- 
r.etric recirculating flow in £_ractical Jjsothermal £ombustor geometries) 
nas been developed and employed in the theoretical aspects of this inves- 
tigation. It is used in a parametric study to predict niean velocity pro- 
files and streamline patterns under various flow conditions. The develop- 
ment of STARPIC began with the nonswirling TEACH (peaching elliptic £xi- 
symmctric £haract er i St i cs £eur ist ical ly) computer program (60) as opposed 
to that of Li I ley iAA) because the latter contains extensive combustion 
calculations which are i nappropr i ate under the present objectives. 


Besides the incorporation of all swirl momentum effects, numerous special 
features were included, such as: (1) a generalized "stairstep" simulation 
of the sloping combustor wall, (2) an advanced nonuni form grid system, 

(3) momentum wall functions derived from a recent experiment, (4) use of 
advanced estimates of ini*'ial field values, and (S) streamline calcula- 
tion and plotting. 

» • 

Discussion of this finite difference computer program is minimized 
here because it is extensively described in a recent NASA report (11). 

The report deals with the complete computational problem, showing how the 
mathematical basis and computational scheme may be translated into an 
effective computer program. A flow chart, FORTRAN IV 1 i sting, notes about 
various subroutines, and a user's guide are aiso supplied. 




3.1 The Governing Equations 


The turbulent Reynolds equations for conservation of nass , rnomentum 
(in X, r, and 6 directions), turbulence energy k, and turbulence dissipa- 
tion rate c which govern the two-dimensional, axi symmet r i c, swirling, 
steady flow may be taken as previously (1*0, k2, kk, kj ^ ^8# ^9, 6l). The 
transport equations are all similar and contain terms for convection and 
diffusion (via turbulent flux terms) of a general dependent variable 
The source term S contains terms describing the generation (creation) 
and consumption (dissipation) of 0. Introducing turbulent exchange co- 
efficients and the usual turbulent diffusion-flux (stress - rate of strain 
type) laws, it can be shown that the similarity between the differential 
equations and their diffusion relations allows them all to be put in the 
common form: 
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The forms of the source term S. are given in Table I and other quant i- 
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ties are defined as follows: 
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(3. 2d) 


Implicit here is the use of the two~equation k~c turbulence model 
(62, 63) and constant Prandtl -Schmidt numbers from which exchange coeffi* 
cients are calculated: 


y * C^pk /e + y^^ 


r. - v/o. 

4> ❖ 


(3.3) 

(3.4) 


The total equation set must be solved for the time-mean pressure p, 
velocity components u, v, and w, etc. Then other useful designer inform- 
ation, like streamline plots, breakaway and reattachment points, recircu- 
lation zones and stagnation points, for example, may be readily deduced. 


3.2 The Finite Difference Formulation 


J 


Solution may be via the stream funct ion-vorticity, or primitive 
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pressure-vtiocity approach. Mtaraas tha formar approach, usad in tha 
1968 cofflputar pro9ram from Imperial Collage for example reduces by 

one the number of equations to be solved and eliminates the troublesome 
pressure (at the expense of trouble with the vorticity equation), the 
preferred approach now is SiMPLE (mnemonic for semi*hnpl icit method for 
£ressure Hnked equations) which focuses attention directly on the lat> 
ter variables. 

Because it possesses many advantages, the present work has been 
developed using this new technique, which was incorporated in the TEACH 
program (60). 

The finite difference equations are solved on a complex mesh system, 
shown in Figure The intersections, the point P for example, of the 
solid lines mark the grid nodes where all variables except the u and v 
velocity components are stored. The latter are stored at points which 
are denoted by arrows (and labeled w and s, respectively) located midway 
between the grid intersections. Details of the special merits of this 
staggered grid system have been reported previously (60, 6$). The dif* 
ferent control volumes C, U, and V, which are appropriate for the P, w, 
and s locations, respectively, are also given in the figure. Thus the 
u* and v'cell boundaries lie either exactly on or exactly halfway be- 
tween the grid intersection points P. This improves the accuracy of cal- 
culating pressure gradient effects on the u and v velocities (as pressure 
values now lie directly on appropriate cell faces) and of calculating 
cert8^in cell boundary fluxes. The flowfield domain boundaries are posi- 
tioned midway between grid points so that they coincide with normal velo- 
cities as in Figure 5, which shows an example of grid specification for 
the geometry under consideration. 
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In obtaining tho flnltt dlfforoneo equatlona using the primitive 
epproeehi it is eommon to use e hybrid formuletion, which is e combine* 
tion of the centre! end upwind differencing scheme. It Is eppifed to 
the representetion of the convection end diffusion terms end provides 
(when needed) the required degroe of upstream differencing. In this man- 
ner the finite difference equation for a typical variable becomes 

•I *P • tj *J ♦ sj (’-S' 

where 

a* ■ E a* - 

£j - sum over N, S, E, and W neighbors 

for each variable d at each point P of a grid system covering the flow 
domain. 


3.3 The Solution Procedure 

The fiowfleld is covered with a nonuniform rectangular grid system; 
typically the boundary of the solution domain falls halfway between its 
immediate nearby parallel gridlines: and clearly specification of the x 
and r coordinates of the gridlines, together with information concerned 
with the position of the boundary, is sufficient to determine the flow- 
field of interest. The curved or sloping boundaries are simulated by 
(Mans of a stairstep approach with sloping boundary segments where appro- 
priate. 

Finite difforoncQ equations like Equation ( 3 .S) hold for all vari- 
ables at all internal points of the mesh system, except those adjacent 
to wall boundaries. At these points, the correct boundary Influence is 


zk 

tnsert«d by a llnaarUad aourca tachniqua, uainq wait functions so as to 
avoid dotal lad calculations In naar-wall raqions. Tha affact of swirl on 
wall function spaciflcatlon is handled as follows. Tha previous ideas 
are extended to find total, tangential wall shear stress near boundaries 
(involving x and 6 directions for an r • constant wall, and r and 9 direc* 
tions for an X ■ constant wall). Then appropriate components are deduced 
directly (for the u and w velocities which are tangential to an r ■ con- 
stant wall, and v and w velocities which are tangential to an x "constant 
wall). The effects on u, v, and w moroantuin equations are incorporated 
via the usual linearized source technique (II). 

The finite difference equations and boundary conditions constitute a 
system of strongly-coupled simultaneous algebraic equations. Though they 
appear linear they are not, since the coefficients and source terms are 
themselves functions of some of the variables; and the velocity equations 
are strongly linked through the pressure. In the solution procedure alge- 
braic equations like Equation (3.S) are solved many times, coefficient 
and source updating being carried out prior to each occasion. The prac- 
tice used here is to make use of the well-known ^rl -diagonal Mtrix algo- 
rithm (TOMA), whereby a set of equations, each with exactly three unknowns 
in a particular order except the first and last which have exactly two un- 
knowns, may be solved sequentially. In the two-dimensional problem one 
considers the values at gridpoints along a vertical gridline to be unknown 
(values at P, N, and S for each point P), and the most recent values at 
each E and W neighbor are considered known. The TOMA is then applied to 
this vertical gridline. In this manner one can traverse along all lines 
in the vertical direction sequentially from left to right of the integration 
domain, using a certain degree of underrelaxation to enhance convergence. 


CHAPTER IV 


EXPERIMENTAL FACILITY AND TECHNIQUES 
A. I Wind Tunnel 

The Oklehome Stete University swirling confined Jet test fecility is 
used in the experimentel study. A schemetic of the overall fecility is 
shown in Figure 3> Ambient eir enters the low*speed wind tunnel through 
an air filter consisting of foam rubber. Tnen it flows through an axial 
flow fan whose speed is continuously variable over a wide range in that 
it is driven by a S h.p. varidrive motor manufactured by U.S. Electrical 
Motors, Inc. Next, the flow gradually expands with the tunnel cross sec- 
tion without separation because numerous fine mesh screens are encounter- 
ed along the way. 

The flow conditioning section is next, it consists of: a perforat- 

ed aluminum plate (2 mm diameter holes), followed by a fine mesh screen, 
a 12.7 cm iength of packed straws, and five more fine mesh screens. Next 
is the contoured nozzle leading to the test section. This axisymmetric 
nozzle was designed by the method of Morel (66) to produce a minimum ad- 
verse pressure gradient on the boundary layer to avoid the flow unsteadi- 
ness phenomenon associated with local separation regions. The area ratio 
of the flow conditioner cross section to that of the nozzle throat is 
approximately 22.5. Further details are available elsewhere regarding 
the design and construction of the nozzle, flow conditioner, and ideal- 
ized test model which is described in the next section (67, 68). 
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6.2 T«tt Section 

The test section consists of e swirl vene essenibly end en ideelited 
combustion chember model which is mounted onto the throat of the wind 
tunnel nozzle. A photograph of the general arrangement in the laboratory 
is given in Figure 6. The swirl vane assembly contains ten brass vanes 
which are individually adjustable for any vane angle The diameter of 
the flow passage within the swirier matches that of the contour nozzle 
throat which is approximately IS cm. A solid hub of 6 cm diameter is 
located at the center of the swirier with a streamlined nose facing up- 
streamt as seen in Figure 7. The downstream end is simply a flat face, 
simulating the geometric shape of a typical fuel spray nozzle. Figure 8 
shows the swirier assembly from the downstream side. A special Jig, also 
seen in the figure, is employed to set each vane at the desired angle. 

The idealized .combustion chamber model is composed of an expansion 
block and a long plexiglass tube. The expansion block is a 30 cm diameter 
disk of oak wood which is mounted to the downstream section of the swirl 
generator. The block has a IS cm diameter hole centered on its axis 
through which the air entering the model flows smoothly. The downstream 
face of the expansion block has been shaped to provide the desired flow 
expansion angle a which is shown in Figure 9. There are currently three 
interchangeable expansion blocks and the appropriate choice gives a ■ 90, 
70, or degrees. Figure 9 shows a view of the swirier with two of the 
expansion blocks. 

The actual combustor of a gas turbine engine, which is illustrated 
in Figure I, is idealized in the present study as there are no film cool* 
ing holes or dilution air holes, and the simulated chamber wall consists 
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of a straight pipa. Tha inside diamatar of tha pipa is 30 cm and tha 

iangth is approximataiy 12$ cm. Tha substantiai sixa of tha tast modal 

providas axcallant proba rasolution for fiva*hoia pitot proba maasura* 

roants, and it also anhancas flow visualization. Plaxiglass was tha 

chosan matariai fOr tha tast saction pipa in ordar to facMitata flow 

visualization. It Is supports.’ by a woodan eradla which is mountad on a 

• ' 

table. 

Tha tast saction is carefully aligned with the wind tunnel so that 
tha tast section and wind tunnel centerlines are colinear. The swirl vane 
assembly and the expansion block are always securely centered on the wind 
tunnel centerline, as they are both positioned on the discharge end of 
the contoured contraction with dowel pins. The test section inlet is 
centered since it slips snugly over the expansion block which wa** care~ 
fully machined with a lathe. The test section discharge, however, is 
aligned with a iowpower Spectra*Physics Stabiiite, model 120, laser beam 
which passes through its center. A plexiglass disk which fits tightly in 
the model dischar^ identifies this center. 

A removable length scale exhibiting axial distances downstream is 
mounted above the test section for flow visualization photography. 

k.3 flow Visualization 

Flow visualization experiments using three different techniques were 
conducted during this investigation and each is described in the follow* 
ing subsections. Such experiments serve as an important investigative 
tool by immediately providing a clear view of the flowfield. During the 
initial phase of this study, it was employed: (1) to examine the overall 

flowfield; (2) to detect unforeseen difficulties with measurements; (3) 


to check the gross cherecteristics of the computer predictions; (k) to 
locate possible regions of interest in which additional probe measure- 
ments might be desired; and (5) to help Interpret measured data. The un- 
foreseen difficulties with measurements mentioned above include time- 
dependent instability effects, regions of severe turbulence, or probe 
interference. 

Probably the most widely used flow tracer to date is smoke. Smoke 
particles are very small, typically less than 1.0 x 10 ^ mm in diameter. 
Hence they follow air motion quite well; however, they cannot be viewed 
or photographed individually. Thus one can only create a relatively 
dense cloud of particles and view the motion of its edges. The cloud 
will diffuse if the air flow is turbulent, the edges becoming indistinct 
and then disappearing as the particle concentration drops. An attractive 
flow visualization technique for both detailed as well as overall charac- 
terization in strongly turbulent flows employs neutral ly-buoyant helium- 
filled soap bubbles. Although much larger than particles of smoke, soap 
bubbles are almost ideal because they trace particle paths and have the 
average density of air (5). 

k.3.1 Neutral I y-Buoyant Soap Bubbles 

Apparently soap bubbles were first used as flow tracers by H. C. H. 
Townend, whose work was reported by Lock (69) in 1928. Coal gas was used 
to make the bubbles neutral ly-buoyant for use in a propeller flowfield 
study. Redon and Vinsonneau (70) injected air and also hydrogen into 
soap bubbles In a wind tunnel and concluded that air-filled bubbles gave 
essentially correct trajectories. Some statistical measurements of tur- 
bulence using bubbles in a wind tunnel were reported by Kampe de Feriet 
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(71) • Of all tha previous work with soap hubbies, that most related to 

the present study was conducted by Owen (72), who employed them to study 

« 

the fiowfieid in a cylindrical vortex tube. 

The bubble flow visualization system used in this study is composed 
of (i) the bubble generator, (2) the light source, and (3) two cameras 
and tripods. The arrangMient of soma of these components is schematical* 
ly Illustrated in figure 10. The bubble generator is model 3, manufactur* 
ed by Sage Action, Inc. It consists of the injector, shown In Figure 11, 
in which the bubbles are actually formed, and a console that supplies the 
constituents to the Injector. Only the injector is situated In the tun* 
nel. It is located far enough upstream of the test section to avoid any 
disturbance to the flow pattern under study. The injector currently em> 
ployed is 1.6 cm in diameter and approximately 11.5 cm long with a 1 cm 
diameter stem perpendicular to the main body. Compressed air at a gauge 
pressure of A. I bar Is supplied to the injector through this stem, while 
helium at a gauge pressure of 1.6 bar and bubble film solution are intro* 
duced through smaller tubes directly upstream of the stem. The bubble 
film solution is ^ special soap mixture, which is commercially available 
from Sage Action, Inc. After formation, 0.5 to 1.0 mm diameter bubbles 
are injected into the surrounding freestream and reach the freestream 
velocity within a short distance. They are then carried by the air stream 
through the test section. The bubbles are injected at a fairly even rate, 
it is estimated by Hale et al. (73) from photographs that this bubble pro* 
duct ion rate is as high as 500 bubbles per second. 

The bubble generator console is placed in a convenient location near 
the tunnel. Three lengths of plastic tubing carry the compressed air, 
helium, and bubble film solution from the console to the injector. The 
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console contains micromataring valvas to control tha constltuant flOM* 
rates, a prassurixa4 cylinder to hold tha bubble film solution, and Inlet 
connections for high-prassura air and helium. 

An American Optical Co. slide projector, modal 1389, served as the 
light sourM for Illumination of the bubbles. A vertical sheet of light 
k cm thick was produced to illuminate the rx'plane of the bubble flow 
pattern. This was provided by using a cardboard slide which is opaque 
except for a thin slit cut out for light passage to the tc3t section. 

Tri-X Pan, a very light-sensitive photographic film rated at ASA 
AOO, along with a large camera aperture of f 2.0, were employed in order 
to obtain photographs of acceptable contrast. Moreover, the film was 
substantially overdeveloped to compensate for underexposure; this is 
equivalent to using a severely light-sensitive film rated at approximate- 
ly ASA 6,000. 

A typical procedure begins with adjustment of the flow metering 
valves for air, helium, and bubble film solution In order to obtain a 
high generation rate of small bubbles. Then e Polaroid Land Camera, model 
l80, is used to obtain preliminary photographs since- Immediate results are 
obtainable. Once the lighting and camera orientation, camera settings, 
and glare are resolved the desired production photographs are taken with 
a 35 mm Minolta single lens reflex camera, model SRT200, with a kS mm 
lens. 

Each bubble leaves an image of its path, or a streak, on the film 
and this streak becomes a record of the air movement within this time 
frame. The record of a number of streaks on a single photograph produces 
a history of the air motion. The intensity of each streak on the film is 
not dependent on shutter speed, but rather on the aperture and bubble 
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velocity, M well e$ the sensitivity of the film Itself. The length of 
pethllnes and (for the slower shutter speeds) the number of pathlines 
are, however, strongly dependent on shutter speed. A relatively large 
number of streaks provides an overall picture of the flow pattern. A 
much smaller number though Is useful to see details of the air motion 
that are evident only by following individual streaks for an appreciable 
distance. Production photographs were typically taken a.t Shutter speeds 
ranging from 1/60 to 1/k second for each flow condition. 

k.3.2 Smoke-Wire 

Although smoke particles may not be observed individually, small and 
discrete smoke streaklines can be located within the flowfield so that 
small scale details may be studied. For example, zero velocity points 

1 within a large recirculation zone may be identified. 

This is accomplished with the smoke~wire technique, originally devel- 
oped by Raspet and Hoore in the early l9S0s. it has been improved and ex- 
tended since then to generate sheets of very fine streaklines (7^, 75, 76). 
However, much of the present smoke-wire capability results from the devel- 
opment of Nagib and his coworkers (77, 73). 

The method consists of a fine wire positioned in the flowfield, coat- 
ed with oil, and heated by passing an electrical current through the wire. 

I As the wire is coated with the oil, small beads of the oil form on the 

! 

wire and at each of these beads the smoke filaments originate when the 
wire is heated. Cornell (79) studied the materials used in generating 
"smoke" for flow visualization, and explained that particles created with 
this technique are actually a vapor-condensation aerosol. That is, they 
are very small liquid particles near 1.0 x 10*^ mm in diameter rnthcr 
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than products of a contuation procots. HovMvar* tn kaaping with tha lit- 
•ratura, thasa particlaa art rafarrad to as sinoha. 

Tha apparatus, soma of which Is shown In Figura 12, consists of a 
0.28 mm diamatar and 0.6 m long constantan wira, a woodan scaffold sup- 
port frame, a 6 cc madical syringe, a 0.45 hg weight, and a solid state 
power supply, modal TR36-4N, manufactured by Electronic 'Research Associ- 
atas, Inc. Tha support frame holds tha medical syringe oil reservoir on 
the upper horizontal member, allowing tha wire to pass through tha syringe 
needle, the test section, and finally through the lower horizontal member. 
The frame is easily relocated so that local details at virtually any 
axial station can be examined. 

The wire is kept taut as the weight is attached to its lower end and 
is free to move as the wire elongates from the/.n^l expansion. The elec- 
trical leads of the power supply are connected to each end of the constantan 
wire, and conduct about 3>S amps at 25 volts d.c. The same slide projec- 
tor described for the soap bubble experiment provides the vertical sheet 
of light to illuminate the rx-plane. 

The smoke-wire technique requires air velocities of less then approx- 
imately 10 m/si otherwise the oil droplets are blown off the wire before 
It is properly coated, and also the smoke diffuses too quickly. Hence 
the smoke experiments were conducted at lower air flow rates than other 
experiments. 

Tha procedure begins with applying a slight force on the syringe 
plunger, which causes several oil drops to run down the wire, coating it 
with small beads while the wind tunnel is not running. A uniform coating 
of fine oil beads is necessary, as each bead is a source of a distinct 
streakline. Several types of oil were tried, aittough tighiy refined 
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kerosene is utiUzed for ell production photoprephs. Then the wind tun- 
nel motor is turned on. The photoprapher waits approximately I second 
before manually activating the shutter while the required heat is suppli- 
ed to the oil. A heavy sheet of streaklines then appears, typically , for 
a duration of about O.S second. Only slight experience with this system 
is required for effective photography. Further details concerning con- 
struction of the apparatus and photography of the nonswirling flowfields 
are found elsewhere (80). 

The smoke-wire has primarily been used in considerably less compli- 
cated flow situations such as the flow over airfoils and bluff bodies. 

The present swirling flows exhibit very high levels of turbulent mixing, 
causing the smoke to diffuse very quickly. Since only very localized re- 
sults are possible under such flow conditions, this experiment was con- 
ducted at several axial stations for each flowfield under investigation. 

A large camera operture of f 2.0 was used for most production photo- 
graphs. As for soap bubbles, faster shutter speeds such as 1/60 second 
provide a more detailed observation whereas slower speeds, for example, 
1/8 second, give a more time-averaged view. Therefore, shutter speeds 
for each flow condition typically include at least l/60, 1/30, 1/15. and 
1/8 second. The same film development procedure used for soap bubble 
photographs is employed here. 

A. 3. 3 Tufts 

Perhaps one of the simplest ways of visualizing a flowfield is by 
inserting a wire grid with a tuft attached at each intersection point of 
the mesh. These tufts respond to the flow patterns and generally point 
in the local flow dire<.tion. Furthermore, their vascillating motion 
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indicates the degree of turbulence in steady mean turbulent flow situa- 
tions. Very large tuft fluctuations, however, can be produced by the 
combined effects of mean flow unsteadiness and turbulence. Such tuft 
screens have been employed to visualize the vortex shedding behind swept 
wings by Bird (8l) .and Gersten (82). They have also been used by McMahon 
et al. (83) to observe the Interaction of an air jet and a cross wind. 

Observations or photographs from a position normal to the tuft screen 
indicate the flow direction in the plane of the screen. Since the combus' 
tor flowfield entails three mean velocity components, a tuft screen was 
constructed For each of the rx- and r6 -planes. The rx-screen slides into 
the combustor model with only slight wall clearance, to minimize the flow 
disturbance. 

The overall dimensions of the rx-screen in the radial and axial direc- 
tions are 29.7 cm and 1.52 m, respectively. Two pieces of extruded alumi- 
num, with 2.5 cm and O.b cm as cross-sectional dimensions, are employed 
as axial frame members, spaced to provide a flow passage radial height of 
28.5 cm. Four 1.8 mm diameter steel rods provide radial support and are 
intermittently positioned axially. Music wire of 0.33 mm diameter tra- 
verses the frame with a uniform grid spacing of 1.9 cm radially. A uni- 
form axial spacing of 1.9 cm, 3.2 cm, and k.k cm is used for each of 


three axial sections of length 61 cm, 40 cm, ana 46 cm, respectively. 

The tufts consist of 1.1 mm diameter cotton thread and are 16 mm in 
length. In construction, the tufts are consistently tied around each 
wire intersection point and a 3*0 mm diameter rod which was then removed, 
resulting in an oversized loop around each wire intersection giving en- 
hanced freedom of movement. 
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The rO'Screen is composed of similar materials, except that its frame 
is a strip of 24 gauge steel 2.5 cm wide, wrapped into a 29.7 cm diameter 
ring. It has a uniform grid spacing of 1.9 cm in both directions. The 
tufts are identical for both screens. 

The light source previously described for soap bubble visualization 
was used as well as the same photographic equipment and technique of over' 
developing the film. Again, slower shutter speeds indicate the time- 
averaged motion, and speeds of 1/2, i/15, and 1/60 second were typically 
used. 


4.4 Mean Velocity Measurements 
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Many instruments are used for separately measuring the magnitude and 
direction of fluid velocity (84). However, there are only a few instru- \ 

ments capable of simultaneously sensing both magnitude and direction. 

The simplest of these senses pressure at certain locations on the surface 
of a sphere, hemisphere, or some other shape of probe. Spherical pitot 
probes using five holes were used by Janes (85) as early as 1915, Meyer 
and Borren (86) in 1928, and Gutsche (87) in 1931 « Pien (88) in 1958, 
further developed the technique upon realizing the fact that the three 
pressure measurements on a great circle of a sphere uniquely determine 
the velocity component in that plane. 

The five-hole pressure probe is recommended by Beer and Chigier (13) 
for mean velocity measurements in turbulent swirling flows and several 
investigators have successfully used it (23, 89, 90). The probe employed | 

in this study is model 0C-I25-12-CD from United Sensor and Control Corp. 
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It is shown schoinatlcony tn Figure 13 and has a 3«2 mm diameter steel 
sensing tip and shaft containing five tubes. The sensing head is hook- 
shaped to allow probe shaft rotation without altering the probe tip loca- 
tion. At the opposite end of the shaft* which is 30.5 cm long* is a 
manifold for the five tubes. 

There are basically three methods of operations. The first and roost 
direct is to adjust the orientation of the probe until both pairs of 
opposing* direction-sensing pressures are nulled* which means the probe 
is aligned with the local flow direction. With the second method* the 
probe is fixed and the flow direction is determined from the calibration 
reiationship between probe pressures and flow direction. The third method 
is simpiy a hybrid of the first two* wherein the yaw angle in the horizon* 
tal plane is nulled and the pitch angle in the vertical plane Is correlat- 


ed. 

The first method is avoided since complicated probe Orientation 
equipment is required. Further* this equipment would probably disturb 
the flowfieid more seriously than the other methods. The second approach 
is compietely unacceptable. This is because the probe cannot properly 
sense pitch or yaw angles relative to its axis greater than approximately 
60 degrees, and the yaw angle often exceeds this limit. This limit is 
caused by flow separation at the probe tip. Therefore* the third method 
is exclusively employed in this investigation* as pitch angles are always 
less than 60 degrees except in a very small fraction of the measurement 
locations. 

The instrumentation system* in addition to the fiverhole pitot probe, 
consists of a manual traverse mechanism* two five-way ball valves* a very 
sensitive pressure transducer* a power supply* and an integrating volt- 
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meter. The probe is mounted vertically on the test chamber* as shown in 
Figure 14, using a manual traverse mechanism, model CI000*l8 from United 
Sensor and Control Corp. This mechanism, shown in Figure 15, is made en* 
tirely of steel with a linear vernier accurately readable within ±0.25 mm; 
This allows the probe to vertically traverse across the chamber radius 
with the capability of manually rotating the probe about its shaft axis 
to null the yaw angle felt by the pressure sensing tip. This yaw angle 
is read from the rotary vernier of the traverse unit which is accurately 
readable within ±0.2 degree. 

The differential pressure transducer is model 590D from Datametrics, 

3 2 

Inc. It has a differential pressure range from 0 to 1.3x10 N/m . The 
sensing element is a high precision stable capacitive potentiometer in 
which the variable element is a thin, highly prestressed metal diaphragm. 
Power for the transducer is provided by a model 721A Hewlett Padtard 
power supply. 

The pressure transducer output is read as the d.c. signal from the 
TSI model 1076 integrating voltmeter. Use of an integrating voltmeter 
removes pressure fluctuations from the vibrating tygon tubing connecting 
the probe to the valves and transducer. The five-way ball valves are 
model SS-432F2 by Whitey Valve Co. These valves along with the pressure 
transducer are mounted on a small portable platform, shown in Figure 14, 
in order to minimize the length of tygon tubing for improved dynamic re- 
sponse of the measurement system. Thus the traverse unit and platform 
are moved together to each successive axial station. 

Finally, auxiliary equipment is used, including model 631-B strobo- 
tac from General Radio, Inc., to measure fan speed. Also, a barometer/ 
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thermoflieter unit from Cenco Corp. is used for local pressure and temper* 
ature readings. 

I|.lt.2 Measurement Technique 

The basic measurement technique entails aerodynamical ly nulling the 
yaw in the horizontal plane by rotation of the probe about its vertical 
shaft and then reading two differential pressures (P|^ - p^ and p^ * p^). 
These pressures along with the yaw angle 6 are used in data reduction to 
obtain pitch angle £ in the vertical plane and the three mean velocity 
components. The data reduction employs two calibration curves which were 
obtained for a single calibration velocity. The underlying principle is 
that the calibration is independent of probe Reynolds number Re^, which 
is based on probe tip diameter. 

Careful calibration experiments reveal that this condition exists 

for Re I 1090, or a local velocity of 5. 4 m/s. Hence measurements of 
P 

such low velocities suffer from a necessary calibration error. However, 
this error affects the velocity measurements typically by less than 6 
percent for Re^ i 400, corresponding to a local velocity gre:<ter than 2.0 
m/s. 

Prior to production measurements, the five*hole pitot rotary vernier 
must be zeroed for yaw so that the x and 6 axes of the measurement coordi* 
nate frame coincide with those of the test section, which are illustrated 
in Figure I4. This is accomplished by rotating the probe until the yaw 
is aerodynamical ly nulled near the center of the test section inlet for 
nonswirling flow. The radial axes of the probe and the test section coin- 
cide since the traverse unit base has been carefully machined to mount 
vertically on the plexiglass tube. 
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Each travarsa begins at tha eantarlina, and this location Is vartflad 

If tha proba casts Its shadoM In tha cantar of a lasar baam which passas 
** 

upstraam along tha cantar of tha tast saction. If nacassary, tha travarsa 
mechanism Is carefully shlmmad until tha above condition Is obtained. 

After tha pressure transducer Is zeroed, tha measurament procedure 
for each location within a travarsa begins with rotation of tha proba un* 
til tha yaw is aerodynamical ly nulled. This is indicated by a zero read* 
ing for where the pressures are identified in Figure I3> The 

resulting probe rotation angle is read from the rotary vernier. Then 
the five*way switching valves are set so that P|^ ~ is sensed by the 
transducer. Finally, the reading of P^ ** Py 's similarly obtained. 

Hass flow rate entering the test section must be determined for the 
calculation of inlet Reynolds number Re^ and spatial mean inlet velocity. 
This quantity is obtained for each flowfield from a measurement of dyna* 
mic pressure with a conventional pitot*static probe. The probe tip is 
located on the centerline and slightly upstream of the contour nozzle 
exit. It has been experimentally verified with a radial traverse of the 
nozzle discharge that, in this facility, the centerline exit velocity is 
extremely close to the spatial mean velocity. This results from the fact 
that the boundary layer momentum deficit is counterbalanced by a slightly 
higher*than*average momentum value Just outside the boundary layer. Thus 
the nozzle exit spatial mean velocity for each flowfield Is obtained from 
a single pitot-static reading. Finally, a continuity correction for the 
slight flow area differe.nce between the measurement and exit stations pro- 
vides the necessary value. 

The ambient pressure and teo^rature were monitored near the facility 
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for the determination of density. Also,. the wind tunnei fan speed was 
frequently checked with the strobotac. 

An Initial experimbnt was performed to investigate the asymptotic in* 
variance of Re^ on the flowfield. This is needed to ensure that the 
dimensionless measurements will be applicable to actual combustor hard- 
ware, which often operates at considerably higher Re . values than those 
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desired for the present test facility. Measurement traverses were per- 
formed at selected axial stations at increased values of Re^ for each 
flowfield. Adjustment of the varidrive fan motor provides excel ient and 
stable control of Re^. The resulting velocity profiles are discussed in 
Chapter VI. 


A. A. 3 Data Reduction 

The differential pressure readings from the five-hole pitot probe 
are utilized directly to obtain the square of the vector velocity. In 
turbulent flow conditions, the time-mean pressures give rise to the time- 
mean average of the square of the vector velocity V^. Since 

V^ - ^ (A.l) 

where V is the fluctuating portion of the velocity magnitude and the 
overbars denote time-averaging, it is slightly incorrect to infer that 
is equal to the square of the magnitude of the time-mean velocity vector 
V . However, the fluctuation term V'^ is not known and very little in- 
formation is available for the effect of turbulence in swirl flows on 
pressure probes, which is probably considerable for turbulence Intensi- 
ties greater than about 20 percent (12). Furthermore, the procedures for 
making corrections for turbulence levels are long and tedious, and even 
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then the confidence In their applicability is unknown. Therefore, no 
attempt Is presently made to incorporate such corrections, and the deduc- 

W ** 

ed velocity is taken to be the time-mean velocity magnitude V which is 
written without the overbar from here onward. 

The data at each measurement location are reduced with the FORTRAN 
computer program listed in Appendix C by first calculating the pitch co- 
efficient (Pn “ « cubic spline Interpola- 

tion technique is used to ol>tain the pitch angle 5 in the vertical plane 
from the calibration characteristic presented in Figure 16. The result- 
ing value of 6 is similarly utilized to determine the velocity coeffi- 
dent pV /[2(pj. - p^)] from the calibration characteristic given in 
Figure 17. 

Values for V as well as the axial, radial, and swirl velocity com- 
ponents u, V, and w, shown in Figure l8, are easily calculated from the 
velocity coefficient, pitch angle, and yaw angle 3, which is in the hori- 
zontal plane. The latter angle is given by 

3 • 360* - (4.2) 


where ij; !s the probe rotation angle read on the rotary vernier of the 
traverse mechanism. The magnitude of the velocity vector is given by 


V 


p 2{p^ - Py) 




and the velocity components are obtained from 
u • V cosfi cos3 
V • V slnfi 


(4.3) 


(4.4a) 

(4.4b) 


and 
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w ■ V cos6 sine 


(4.4c) 


Further calculations include that for mass flow rate and axial flux 
of angular momentum, which are respectively expressed as 


m ■ 2it / p ur dr 
o 


(4.5) 


-X 2 

4ft • j P uw r dr 


(4.6) 


These are evaluated from the velocity measurements with the data reduc- 
tion computer program. 


4.4.4 Calibration 


Proper calibration of the five-hole pressure probe is extremely 
critical to the accuracy of the experimental work. Consequently, consid- 
erable care was exercised in the construction and use of the calibration 
equipment. The calibration equipment consists of a small air jet, a ro- 
tary table, a probe mounting bracket, and the instrumentation system pre- 
viously described. The calibration Jet supply line consists of a com- 
pressed air line, which delivers the desired flow rate through a small 
pressure regulator and a Fischer and Porter model 10AI73SA rotameter. 

The jet housing consists of an effective flow management section follow- 
ed by a contoured nozzle with a 3*5 cm diameter throat. 

The rotary table is model BH-9 from Troyke Manufacturing Co. Its 
rotary vernier is readable within ±0.5 minute. As shown in Figure 19, 
the aluminum probe mounting bracket is secured to the rotary table, and 
it supports the probe which rests in a cylindrical steel collet. 


The motion of the rotai 7 table orients the probe at the desired pitch 

angle, whereas the yew is aerodynamical I y nulled. The probe sensing tip 

remains at the centerline within the potential core of the Jet and less 

than one throat diameter downstream of the nozzle discharge plane. The 

pitch angle is zeroed by very carefully aligning the probe shaft in a 

plane parallel to that of the nozzle throat. 

The calibration procedure consists of recording the voltage output 

from the pressure transducer for differential pressures P|^ ~ And 

P^'Pyt where these pressures are identified in Figure 13. These data 

are obtained at 5 degree increments in 6 over the range -SS” s 6 s 55". 

The measurement technique requires the entire calibration to be conducted 

at a constant Jet velocity. This is permitted since the dimensionless 

calibration coefficients are independent of Re^ for Re^ > 1090 (5.^ m/s) 

P P 

determined by careful calibration experiments. Measurements of much 
lower velocities suffer from a calibration error. However, this is rare- 
ly over 6 percent for Re > 400 (2.0 m/s). 

P 

Figures 16 and 17 show the calibration characteristics from which £ 

2 

and velocity coefficient pV /[2(p^-p^)] are obtained, respectively. 

Both curves exhibit considerable symmetry, as the five pressure sensing 
holes are almost symmetrical about the probe tip axis. 


CHAPTER V 


GROSS FLOWFIELD CHARACTERIZATION 

Recirculation zones are important to combustor designers because most 
of the burning occurs in these regions and they exhibit the highest tempera 
tures (47). The size and location of these regions is determined by flow 
visualization experiments and predicted by numerical computations. The 
predicted time>mean streamline patterns are compared with corresponding 
artistic impressions deduced from flow visualization photographs of tufts, 
smoke, and bubbles responding to the experimental flowfield patterns. Se> 
lected photographs of the visualization experiments are also presented and 
discussed in the subsequent sections. 

S.i Artistic Impressions of Streamline Patterns 
From Flow Visualization 

Photographs of each of the six fiowfieids resulting from ■ 0, 45, 
and 70 degrees with a * 90 and 4S degrees have been examined in detail for 
each of the three flow visualization methods currently employed. The 
characteristics of the overall flowfield are illustrated and discussed via 
the resulting time-mean dividing streamline patterns. These are sketched 
in Figures 20 and 21 from information obtained from the entire collection 
of flow visualization photographs. Results from the smoke-wire experiment 
are utilized near the inlet, whereas tuft and bubble data are used in 
approximating the size and shape of the recirculation zones downstream. 
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Bubble flow petterns also reveal the exUtence of a precasalng vortex 
core, which occurs downstream of the central region. 

The resulting streamlines for the three cases of the sudden expan- 
sion a • 90 degree geometry are shown In Figure 20. The nonswirling flow 
sketch in part (a) exhibits a large corner recirculation zone which Is In 
excellent agreement with the corresponding streamlines from the measure- 
ments of Chaturvedi (1$). Appropriate prior measurements for the present 
swirling flows have not been found. For the moderate swirl vane angle 
case of ^ As degrees a central recirculation region appears in conjunc- 
tion with a decrease in size of the corner zone. A thin precessing vor- 
tex core, discussed by Syred and Beer (91), is observed near the center- 
line extending from the end of the central region to the test section 
exit. The axial location of the upstream varies, ranging approximately 
from x/0 • 1.2S to 1.7S. This vortex core is essentially a three- 
dimensional time-dependent instability which occurs as a swirling region 
of negligible axial velocity whose center winds around the test section 
centerline. A further increase in vane angle to ^ • 70 degrees results 
in continued enlargement of the central zone; however, the corner bubble 
is essentially unaffected. For this flowfield the vortex core is expand- 
ed in the radial as well the axial direction, as its upstream edge fluc- 
tuates from approximately x/D ■ 1.0 to l.$. 

The corresponding sequence of dividing streamlines is found in Fig- 
ure 21 for the gradual flow expansion case with a • degrees. As with 
Figure 20, the nonswirling flowfield exhibits excel lent agreement with cor- 
responding •measurement of Chaturvedi. No evidence of a corner zone is found 
from examination of the photogrephs for the Intermediate vane angle case, 
although the flow pattern Is otherwise very similar to the corresponding 


ftow for 0*90 dogroot. Tho flowfiold with d*70 dogroot producot exact* 
1y the seme effects es In the abrupt expansion case of Figure 20. 

5*2 FIom Visualization Results 

5.2.1 Tufts 

Tuft visual tzetion Is very Important In that It supplies an overall 
view of zones very clearly. As discussed In section photographs 

at various shutter speeds were obtained. Slower speeds show more of the 
temporal behavior, although the tufts are sometimes not distinctly vlsl* 
Die in portions of the flowfield. Some of the more noteworthy photographs 
are presented here, encompassing a range of shutter spMds. Velocities in 
recirculation zones are often somewhat lower than in other portions of the 
flowfield, and thus under such conditions there may be insufficient drag 
on a tuft to align it accurately with the local flow direction. However, 
this is taken into consideration in interpreting the photographed results. 

The three sudden expansion flowflelds ere characterized in Figure 22 
as photographs of the rx'screen tufts. Figure 22(a) identifies the corner 
zone reattachment point for the nonswirling flowfield. It is found from 
the length scale as x/0 “ 2.1, wherein tufts showing no air motion are 
considered stagnation points. Further, an Indication of the turbulence 
level is depicted ar the rather slow shutter speed of 1/15 second was 
used. 

For swirl vane angle d ■ A5 degrees. Figure 22(b) is a photograph 
taken at 1/125 second. The swirl is in the counterclockwise direction 
when viewed from downstream. Both the corner and central recirculation 
regions are clearly visible, and their downstream stagnation points appear 
to be at x/0 ■ O.k and x/0 * 1.6, respectively. The large vane angle case 


U photo9r«ph*4 In Ft9tir« 22(c), •9«ln at t/125 Mcond shutter speed. Al- 
though the iq>per corner is not Mil illumineted in this photograph, the- 
loMer corner region Is seen to extend to x/0 ■ O.li end the central asone to 
x/0 ■ 1.75. 

Photographs of the fiowfields with d ■ AS degrMS ere shown in Figure 
23. The nonswiriing cese in Figure 23(e) is elnost identical to Figure 
22(e), with e corner zone length of x/D - 2.0. For ^ • kS degrees, Figure 
23(b) shows no evidence of e corner zone using 1/125 second shutter speed 
end the central region apparently extends to x/0 ■ 1.5. Finally, the p ■ 

70 degree cese exhibits e central region extending downstream as far as 
x/0 ■ 1.85 for 1/15 second shutter speed. 

Figures 2A end 25 ere photographs of the re-screen of tufts for a * 

90 and p ■ A5 degree flow conditions at x/0 • 0.5, i«0, 1.5, end 2.5. The 
shutter speed Is 1/2 second in order to capture information regarding tur- 
bulence levels as well as flow direction. The camera Is positioned slight- 
ly downstream end on the centerline of the test facility, and hence an 
annular ring of distorted tufts appears around the tuft screen because of 
the reflectivity of the plexiglass tube. The variation in brightness of 
the tufts is merely caused by the angle at which the light strikes each 
particular tuft. 

Figure 2A(a) apparently shows noticeable but random asymmetry in tur- 
bulence levels at x/D ■0.5, whereas Figures 2A(b) and 25(a) seem to re- 
veal a horizontal strip of highly fluctuating tufts. However, it may be 
unreal istc to conclude from this that such a region of intense turbulence 
exists because the effect of gravity acting on tufts pointing vertically 
upward or downward is different from that on horizontal tufts. Figure 
29(b) shows that at x/D ■ 2.9, lar^ diraction fluctuations occur only 


near the centerline. This probably results from' the tendency of the cen* 
ter of the princess Ing vortex core to randomly oscillate about the test 
section center 1 1 ne \9i). 

5.2.2 Smoke-Wire 

Local details in the nonswlriing flowfields are clearly revealed 
through the visualization of streak! ines indicated from the generation of 
illuminated smoke, as discussed in section 4.3.2. In the swirling flows 
strong mixing diffuses the smoke so that streaklines are not distinguish- 
able. However, under such conditions recirculation zone outlines are 
visible, especially in the region near the smoke-generation wire. A 
selected photograph is exhibited and discussed for each of the six flow- 
fields considered in previous sections. 

The corner recirculation bubble in the nonswirl inn case with wail ex- 
pansion angle a » 90 degrees is revealed in Figure 26(a) using 1/15 second 
shutter speed. Also, the radial location of the zero velocity point with- 
in the upper and lower corner bubble is estimated to be approximately 
0.15 0 from the respective waits of the test facility. This agrees with 
the velocity measurements presented in Chapter VI. 

The moderate swirl vane angle flow with a 1/30 second exposure is 
shown in Figure 26(b). The shortened corner zone is easily identified 
because the adjacent downstream flow contains no smoke near the inlet. 

This bubble is seen to extend to approximately x/0 » 0.45. The tufts in 
Figure 22 indicate a slightly shorter zone ending at approximately x/D - 
0.4. The upstream portion of the central zone is also clearly seen, as 
low velocity fluid carries a dense mass of smoke which slowly moves up- 
stream of the inlet. Further, the precessing vortex core is seen to 
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contain no smoke. Since it exhibits negligible axial velocity, the smoke 
is essentially carried around it by the high velocity fluid outside the 
core. 

Figure 26(c) is a photograph at 1/8 second of the ^ « 70 degree case 
wherein this core is not as distinct. In this case some smoke has diffus- 
ed into the core due to a slightly longer delay before activating the cam- 
era shutter. Observe that both the corner and central zones near the 
inlet reveal that the only change from those for ^ ■ 45 degrees is a 
si ightly wider central region. 

The same swirl sequence for the gradual expansion geometry with a ■ 

45 degrees is presented in Figure 27. The nonswirling case at 1/15 second 
reveals about the same details as the corresponding case in Figure 26. 

The (j> « 45 degree flowfield at 1/60 second exhibits no evidence of a cor- 
ner bubble and a thinner central bubble at the wire position of x/D « 

0.25. The existence and location of the vortex core is very similar to 
that in Figure 26(b). The large vane angle case shown in Figure. 27(c) 
using 1/15 second reveals no evidence of a corner zone. Also, the central 
zone at x/C- * 0.25 appears considerably wider than that for 4* « 45 degrees 
shown in F'cjre 27(b). 

5.2,3 Neutral ly-Buoyant Soap Bubbles 

Soap bubbles injected into the air flow upstream of the test section 
trace pathlines clearly when illuminated. This technique is discussed in 
detail in section 4.3*1. Photographs of illuminated bubbles are taken 
with various shutter speeds and camera positions. In relatively lower 
turbulence intensity portions of the flowfield mean flow directions can 
be obtained by ensemble averaging local tangents to pathlines traced out 
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by soap bubbles. This helps define the flowfteid geometry In terms of 
the outline of recirculation regions. 

A sample flow visualization photograph is presented In Figure 28(a) 
corresponding to the zero swirl, 90 degree expansion angle flowfield. 

The photograph, taken with a relatively long time exposure (1/8 second), 
clearly shows a great number of individual pathlines. Photographs of 
this type can be used to indicate regions of highly turbulent flow such 
as that near the centerline of the flowfield which exhibits smoother, 
straighter pathlines. In addition, the outline of the corner recircula** 
tion region can be estimated from Figure 23(a) (and numerous additional 
photographs taken at the identical run condition). For this geometry the 
mean stagnation point defining the end of the recirculation zone appears 
to be at x/0 * 2.0. 

A photograph with ^ • kS degrees and 1/3 second shutter speed is 
shown in Figure 28(b), where the precessing vortex core is clearly seen 
extending from x/D ■ i.5 to the exit. Its upstream edge fluctuates ran- 
domly from approximately x/D > 1.2S to 1.75. The corner bubble is observ- 
ed in both the upper and lower portions of the flowfield, extending to 
approximately x/0 0.4, which agrees almost exactly with the smoke flow 

pattern for this flowfield in Figure 26(b). A photograph using 1/15 sec- 
ond shutter speed is presented in Figure 28(c) for the - 70 degree flow- 
field where a thicker vortex core is obvious. Also, it extends farther 
upstream to x/D ■ 1.25 in this photograph, although this fluctuates from 
about x/D - 1.0 to 1.5. The corner zone is faintly visible here, and its 
axial length seems to extend to about x/D • 0.4. 

Figure 29(a) is a 1/15 second exposure of the nonswirling case with 
a ■ 45 degrees. Essentially the same result is found as in the 
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corresponding a • 90 degree ftowfield in Figure 28(a). The intermediate 
vane angie case is visualized using 1 A second shutter speed in Figure 
29(b). As with the smoke flow pattern in Figure 27(b)* there is apparent- 
ly no evidence of a corner bubble. Otherwise* there is little difference 
between Figures 29(b) and 28(b). The ^ * 70 degree flow shown in Figure 
29(c) is photographed at 1/8 second* and agrees with the smoke pattern in 
Figure 27(c) that there is no apparent corner zone. Further* the vortex 
core extends upstream to x/0 > 1.25* as in the corresponding swirl case 
in Figure 28(c). 


5.3 Predicted Streamline Patterns 

The inlet boundary cond-tion for the flowfield prediction program 
consists of a flat velocity profile as in plug flow for both u and w, 
whereas v ■ 0. This approximation has often been employed when pertinent 
measurements are npt available (A6, 47). it was discovered near the end 
of this investigation from the five-hole pitot velocity measurements for 
the swirling flowfie Ids, that in the test facility, the inlet velocity pro- 
files are actually much different in shape and magnitude from those assum- 
ed for the general predictions. These measurements are presented and fur- 
ther discussed in Chapter VI. This difference results from: (a) imperfect 
blade efficiency, (b) the existence of a hub, and (c) the fact that the 
downstream edge of the swirl vanes of the test facility is actually locat- 
ed approximately 4 cm upstream of the flow expansion corner, where x/0 ■> 

0. This swirler location allows the central recirculation zone to begin 
upstream of x/0 • 0, thereby changing the velocity profiles there. Hence 
swirling flow streamline predictions in Figures 30 and 31, parts (b) and 
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(c), are not expected to exhibit exceilent agreement with experiments, 
particulariy near x/D ■ 0. 

Figures 30 and 31 are calcuiated and plotted by computer to show time** 
mean streamiine patterns for a ■ 90 and 45 degrees, respectiveiy. Although 
the solution of the swirling fiows is not entirely reaii Stic, the trends 
are generally correct. The nonswirl predictions do not suffer from the 
inlet boundary condition error described above. The results shown in Fig- 
ure 30(a) exhibit excellent agreement with the flow visualization results 
in Figure 20(a) and also with the streamiine pattern resulting from corre- 
sponding measurements of Chaturvedi (15). The intermediate case of ^ * 

45 degree vanes shows reasonable agreement with the results in Figure 
20(b) except at x/D * 0. Discrepancies include central recirculation zone 
lengths of 1.35 0 and 1.6 D for predictions and visualization, respective- 
ly. The corner recirculation region lengths, however, are in exact agree- 
ment. Surprisingly, for ^ * 70 degrees the predicted and photographed 
central zone lengths are x/D • 3.0 and 1.75, respectively. Further, no 
corner zone is predicted, whereas the photographs reveal one with axial 
length x/D •« 0.4. 

As with the a > 90 degree case, the nonswirling flowfield with a °> 

45 degrees, shown in Figure 31(a), is very well predicted by comparison 
with Chaturvedi *s measurements and the present photographs. No corner 
zone is predicted in this geometry for either the intermediate or large 
vane angle flows, which agrees with visualization results of Figure 21(b) 
and 21(c). For ■ 45 degrees, the predicted central zone length is x/D 
« 1.35 Whereas that photographed is approximately 1.5. The central region 
length for the strongly swirling flow from computations and visualization 
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is x/0 ■ 1.75 and 1.85» respectively. Hence the central zone is more 
accurately predicted for a • kS degrees than for a * 90 degrees. 

S»k Parametric Effects 

Streamline plots, including those deduced from flow visualization as 
well as those predicted, allow recirculation zones to be characterized 
parametrically. The parametric effects are illustrated in Figures 32 and 
33 for predictions and visualization, respectively. Figure 32 shows the 
predicted effects of a and ^ on the corner and central recirculation zone 
lengths. Several observations should be noted. First, zone lengths are 
only slightly affected by a, as found previously for the corner bubble 
under isothermal nonswirling conditions (IS, 31). Further, at these high 
Reynolds numbers, the corner region length is independent of Reynolds num- 
ber (33, 34). As a decreases from 90 to 4$ degrees, the zone lengths tend 
to decrease slightly and the inlet flow is encouraged to impinge on the 
confining walls. Second, the corner bubble length decreases upon increas- 
ing ^ from 0 to 4S degrees, and also the central region and the precessing 
vortex core appear. Increasing the vane angle to 70 degrees slightly en- 
larges the central zone and the vortex core. 

Visualization results exhibit close agreement with predictions for 
the shape of both corner and central regions for ^ • 0 and ^ * 4$ degrees. 
Further, this also occurs for ^ ■ 70 degrees with the gradual expansion 
geometry where a ■ 45 degrees. However, for the abrupt expansion case 
where o - 90 and ♦ ■ 70 degrees, the visualized comer zone is not predict- 
ed. Moreover, the central region length is seriously overpredicted at x/D 
" 3*0 as compared with 1.75. A lack of prediction accuracy for swirling 
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flows is expected because of complicated inlet conditions and/or the lack 
of realism in the turbulence model; however, it is surprising that the a • 
45 degree case exhibits close agreement while that for o * 90 degrees does 


CHAPTER VI 


VELOCITY HEASUREMENTS. PREDICTIONS. AND 
COHPARISON WITH EXPERIMENT 

Each of the six flowflelds under investigation is first character- 
ized in terms of axial and swirl tine-mean velocity measurements obtained 
using a five-hole pitot probe. Then general veiocity predictions using 
estimated inlet conditions are presented. A special test case is de- 
scribed also which serves as a vaiidity check of the computer code for 
the abrupt expansion angle a ■ 90 degrees using measured inlet conditions. 
Next, reliability of predictions and measurements is discussed, and final- 
ly, the results of another validity test for a non-expanding geometry is 
described. 


6.1 Five-Hole Pitot Probe Velocity Measurements 

Details of the five-hole pitot measurement apparatus and technique 
are discussed in section A. 4. Readings were taken at approximately 20 
positions across the test facility radius at the various axial stations 
for which velocities are presented. Before production measurements be- 
gan, readings were taken for ^ ■ 0, 45, and 70 degrees at various values 
inlet Reynolds number. Comparison of these data determined the 
Reynolds number required to ensure that the nondimensional ized measure- 
.nents are independent of Reynolds number. A sample of these data is pre- 
sented and discussed in section 6.4.2. 


Figures 34, 35* end 36 show the axtat and swirl velocity profiles 
for ^ ■ 0, 4S, and 70 degrees with a ■ 90 degrees. Note that the scale 
of dimensionless u/u^ shown at the Inlet in Figures 34(a) and 35(a) is 
different from that of w/u^ which is found In Figure 35(b). These scales 
apply to all figures through section 6.3* Radial velocities are consls" 
tently much smaller than the axial and swirl components, and they are 
indicated In the dividing streamline patterns deduced from flow visual* 
ization experiments depicted in Figures 20 and 21. Therefore, they are 
not presented here. Figure 34 shows the corner recirculation region for 
the nonswirling flow on the verge of reattachment at x/D ■ 1.3 and also 
at 2.0, whereas it occurs at x/0 ■ 2.1 for the corresponding visualiza- 
tion results shown in Figure 20(a). The corresponding measurements of 
Chaturvedi (15) yielded a value of x/D 2.3. 

Both of the inlet profiles in Figures 35 and 36 differ in magnitude 
and shape from the flat profiles used as inlet boundary conditions for 
the general predictions. This is discussed in section 5.3 where it is 
mentioned that the swirl vanes in the facility are not 100 percent effi- 
cient. It can be seen in Figure 35 that a 45 degree swirl vane angle 
produces a maximum swirl flow angle [tan ^(w/u)] of 30 degrees at the 
test chamber inlet. The maximum swirl flow angle for the 70 degree vane 
angle case shown in Figure 36 is only 34 degrees. Thus there is only a 
slight increase of swirl flow angle, although these two flows are differ- 
ent in that the inlet profiles are considerably more sharply peaked for 
the t * 70 degree case. Figures 35 and 36 show zero u and w velocity 
near the axis at the inlet, but actually the probe was insensitive to 
the very low velocities there. This is consistent with the flow 
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visualization result discussed in Chapter V that the central bubble ex* 
tends upstream of the inlet. 

Jt 

The measurements shown in Figure 35(a) provide no evidence regard- 
ing the existence of a corner zone. This is expected because fiow visu- 
alization reveals that the region only extends to x/0 ■ 0.4. However, 
there is clear evidence of a rather large central zone whose length is 
similar to that shown in Figure 20(b). Although the axial velocity pro- 
files are beginning to flatten in the downstream direction, they retain 
a zero velocity value on the axis, which is consistent with soap bubble 
flow patterns as seen In Figure 20(b). The early erratic behavior shown 
by the swirl velocity profile in Figure 35(b) quickly transforms, exhib- 
iting a sol id-body- rotat ion core with a rather fiat profile outside this 
reg ion. 

Figure 36 also reveals the discrepancy regarding inlet velocity pro- 
files. The large central recirculation region causes the downstream 
flow to be accelerated near the top wall. The early erratic behavior 
found in the swirl velocity profile at x/0 ■ 0.5 quickly develops into a 
shape similar to that seen in Figure 35(b). The processing vortex core 
motion discussed in Chapter V results in poor measurement repeatability 
which promotes the irregular behavior within this core region. 

Figures 37, 38, and 39 exhibit velocities for the same sequence of 
flowfieids with a <■ kS degrees. The inlet profiles were not measured in 
this geometry because the presence of the expansion block interferes. 
Effects of on velocities, similar to those found for the sudden expan- 
sion cases, are found in the flowfield sequence for this test section 
geometry. The major difference is that the sloping wall encourages the 


inlet flow to aecelereto near the top wall. Also, it tends to shorten 
or obliterate the corner recirculation region. 

6.2 General Velocity Predictions 

Ail results are obtained via a nonuniform grid system due to the re* 
suiting enhancement of solution accuracy. Nineteen ceils are employed 
in the r-direction, and they are clustered near the shear iayer region 
and along the wail and centerline. The cells are gradually expanding in 
the streamwise direction, and from 21 to 33 are employed as required to 
produce the desired a. The inlet profiles of axial velocity u and swirl 
velocity w are idealized as "flat" (i.e., constant-vaiued) with 100 per* 
cent efficient swirl vanes. Also, it is assumed that v ■ 0. The exit 
is always located well downstream of the central zone where the flow is 
considered para I ie I with zero radial velocity. Further details are 
briefly discussed in Chapter III. For comparison with other results, it 
should be noted that ^ and swirl number S are related by S ■ 2/3 tan 
so that vane angles A5, 60, and 70 degrees, for example, correspond to S 
values of 0.67» I.IS» and 1.83* respectively (12). 

Predicted time-mean axial and swirl velocity profiles for the side- 
wall expansion angle a ■ 90 degrees are shown in Figures 4P, 41, and 42 
for swirl vane angle ^ ■ 0, 4S, and 70 degrees, respectively. The non- 
swirling case in Figure 40 exhibits close agreement with measurements of 
Chaturvedi (IS) as shown elsewhere (S). Note the scale of nondimensional 
u/u^ shown at the inlet profile; this applies to all subsequent general 
predictions. Further, the centerline velocity exhibits little change in 
the streamwise direction. 
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Figure 41(e) reveele that with ^ • 45 degrees, the mean axial veto* 

city profiles are dramatically changed. A maximum velocity value occurs 

in an annular fashion near r/0 ■ 0.3S» although a more flattened shape 

quickly develops before x/D ■ 1.0. It should also be noted that the 

boundary layer on the outer sidewall is relatively thin and is not easily 

seen on the figures using the chosen scale. Note that the scale for w/u^ 

shown at the inlet in Figure 41(b) is different from that of.u/u . Both 

o 

of these scales are used in all subsequent figures through section 6 . 3 . 
All of the swirl velocity profiles exhibit sol id<*body* rotation behavior 
near the centerline, even near the inlet where a flat profile is a speci' 
fied inlet condition. The radial location of the station maximum for w 
tends to increase with x/D in Figure 41(b). Swirl as well as axial velo* 
city profiles appear to approach those corresponding to swirling flow in 
a pipe as x increases. 

Figure 42(a) reveals that the strong swirl case provides a much 
larger central recirculation region at x/D ■ 0.5, which is caused by 
strong centrifugal effects. This promotes a very large foward velocity 
near the well end almost eliminates the corner region, although experi- 
ments show it is unaffected. The radial extent of the central recircula* 
tion zone and the velocity near the wall quickly diminish downstream as 
swirl strength is dissipated. The swirl velocity profile in Figure 42(b) 
are similar to those of Figure 41(b), except much stronger magnitudes re- 
suit from the increased vane angle. 

The same sequence for the a ■ 45 degree geometry is found in Figures 
43 , 44, and 45. The inlet profiles are not shown here since they are the 
same as those of the corresponding cases in Figures 40, 41, and 42. The 
nonswirling flowfield shown in Figure 43 is essentially identical to the 
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nontwirling cast with o ■ 90 dogreot. For ♦ ■ AS citgroat, Ftgura AA 
thowt littia affact of a upon compariton with Figura Al. Tha eantral 
ragion indlcatad In Figura AA(a) It slightly thinnar at x/D ■ O.S. Tha 
strong twirl flowflald llluttratad In Figura AS(a) axhibits a smallar 
backflow ragion than for o • 90 dagraas In Figura A2(a)i although It dis- 
plays largar upstraam valocitlas. Also* tha paak In tha prof! la at x/0 
• 0.5 occurs in an annular fashion, wharaat this is not saan in Figura 
A2(a). Swirl valocitlas in Figura A5(b) axhibit a slightly largar magni- 
tude near tha top wall than that In Figura A2(b), but are otherwise vary 
simi lar. 


6.3 Tha Specific Test Case 

Due to the discrepancy regarding axial and swirl inlet velocity pro- 
files, tha importance of which was discovered near the end of this Inves- 
tigation, another prediction of tha a ■ 90 and ^ * AS degree fiowfield 
is presented. Tha measured velocities are specified as tha inlet boun- 
dary condition, and the results are saan in Figure A6 along with the mea- 
surements. 

This prediction shows that the large paak of inlet axial velocity 
for x/D ■ O.S is located at r/0 ■ 0.2S rather than O.AO as measured. 

Also, a much smallar central zona is predicted than for tha general pre- 
dictions. However, no affact of this discrepancy is seen beyond x/O • 
1.5. The iwasured inlet swirl velocity paak broadens very quickly by 
x/O ■ O.S. Also, a solid-body-rotation core is exhibited there, which 
is somewhat similar to tha previous predictions shown in Figura Al. Com- 
parison with maasuramants shows considerably largar velocity magnitudes 
at x/D ■ O.S, but lower values downstream. 
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This ditmpsncy iMsr Um tnist might rssMlt from bsekfloM at tha 
Inlat boundary which could not ba maasurad baeauea of transient affects 
or tha iMk of probe sensitivity for small veloeltlas. 

6.4 Reliability of Predictions and Maasurments 
6.k.l Predictions 

The finite difference mash size employed for the flowfleld predic* 
tions resulted from a compromise between the irreconcilable alternatives 
of computing cost and flcwfleld resolution. Several grid systmm xrere 
utilized to determine the effect of grid density on the solution of the 
present flowfields. The chosen grid is discussed In Chapter III and an 
example is shown in Figure 5. An increase of 60 percent in the number 
of computational cells produces a change of only 0.3 and l.l percent for 
top wail and centerline dimensionless axial velocities, respectively. 

Another important aspect of solution accuracy concerns the converge 
ence criterion. Rather than use the well established fractional cliange 
criterion, which is a necessary but not sufficient condition for converg- 
ence, a residual source criterion is en^loyed. This residual is defined 
by equations i ike 

- I “•’> 

which has been obtained from Equation (3.1/ and which measures the depar- 
ture from exactness (except for the use of current coefficients and 
source) for the variable d at the point P. When each of these quantities 
becomes smaller than a certain fraction of a reference value, the finite 
difference equations are considered to be solved. 
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Although inoccurocies resulting from mesh resoTution and degree of 
convergence are insignificant^ the predictions suffer from inherent in* 
adequacies In the turbulence model . This is tolerated presently because 
these models for complicated swirling flows are not yet available (49). 
The k-e model developed for nonswirling flows has been extended to in- 
clude swirl shear stresses; however, an isotropic viscosity remains 
while f ..i>ling flows in combustor geometries can exhibit anisotropic be- 
havior (26). However, the test case consisting of the measurements of 
Vu and Gouldin (29) exhibits very encouraging results as discussed in 
section 6.5. 

6.4.2 Measurements 

It is difficult to estimate the accuracy of the flow visualization 
results. The soap bubbles probably trace the flow very accurately be- 
cause they are neutral ly-buoyant and approximately only 0.5 to 1.0 mm in 
diameter. The tufts observed in the rx-plane of the test section do not 
respond well to the local flow direction in a low velocity region due to 
the small drag forces acting on them. They are considered very reliable 
in higher velocity regions, however. The smoke particles have much 
greater density than air and '‘consequently they are affected by gravity 
in low velocity areas and by centrifugal forces in strongly swirling re- 
gions. 

Little information is available concerning the effect of turbulence 
on pressure probes in swirling flows (12). However, it is presumed that 
the five-hole pitot probe is accurate within approximately 5 percent for 
most of the measurements. This value may increase to 10 percent as the 
velocity magnitude falls below approximately 2.0 m/s because the probe 
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is not very sensitive to such low dynamic pressures. In. addition, the 
calibration is no longer independent of probe Reynolds number. The pres 
sure transducer is estimated to have a measurement uncertainty .of approx* 
imately 0.25 percent for a typical velocity magnitude of 5 m/s and 0.75 
percent at 2 m/s. 

An experiment to determine the extent to which the five-hole probe 
disturbs the flow was conducted. A dummy probe of approximately the same 
size and shape was inserted through the lower wall of the test section, 
exactly opposite to the five-hole probe entering through the upper wall. 
Axial and swirl velocity profiles with and without the dummy probe, which 
is traversed exactly like the measuring probe, are shown in Figures kj 
and kB for x/D ■ I.O and 2.5, respectively, with « ■ 90 and ♦ * 45 de- 
grees. Clearly the presence of the dummy probe has virtually no effect. 
Further, since the dummy probe has a negligible influence, these figures 
show evidence of good measurement repeatability. 

Before beginning velocity surveys, an initial experiment was per- 
formed to investigate the asymptotic invariance of Reynolds number Re^. 
This was needed to ensure that dimensionless measurements will be appli- 
cable to combustors which operate at higher Re^ values than those desired 
for the test facility. Measurement traverses were repeated at the same 
conditions for increasing Re^ for each of the flowfields at x/0 * 2.5. 
This axial station was chosen because a probe location slightly down- 
stream of a recirculat ion zone was considered most sensitive to Re and 
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the central rather than corner zone was selected since most of the com- 
bustion and flame stability occur there. A sample of the resulting 
axial and swirl velocity profiles is found in Figures 49(a) and 49(b), 
respectively. Note that there is little change for 68,000 s Re^ s 92,000. 



Hence all measurements of the moderate ^ ■ 4$ degree flowftelds are per* 
formed at Re^ ■ 78,000. Similarly, ail data for ^."6 and ^ ■ 70 degrees 

Q 

are obtained at Re^ ■ 105,000 and 78,000, respectively. 

Figure 50 shows axial and swirl velocity profiles at x/D ■ 1.0 with 
a ■ 90 and 4; • 45 degrees. The measurements above and below the axis 
were obtained from probe entry through the upper and lower wall, respec- 
tively. Although the upper and lower profiles reveal a similar shape, 
they exhibit different magnitudes. This shows clear evidence of a slight 
asymmetry in the flowfield. 

Averaging the mass flow rate calculated from each radial traverse 
in Figure 50 results in a value within 2 percent of that measured at the 
test section inlet. Calculation of mass flow rate from each traverse 
individually yields a discrepancy of approximately 14 percent of the mea- 
sured value. This asymmetry is indicative of the lack of agreement which 
may be expected if mass flow rate is calculated from the upper traverse 
only. 

6.5 Prediction of Coannular Swirling Pipe Flow 
6.5.1 The Flowfields Considered 

Due to the rather poor agreement resulting from the special test 
case near the inlet of the test section, another validity test of the 
prediction code was undertaken. This case, experimentally investigated 
by Vu and Gouldin (29), has no expansion and the inlet velocity profiles 
do not exhibit backflow. Therefore, it is considered less difficult to 
predict. 

This test consists of both counterswi rl inr; and coswirling flows 
from a coannular jet discharging into a straight pipe. A schematic 
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diagram of the flow assen4>ly is illustrated in Figure 51 • The assembly 
consists of two major segments: an inner flow passage and an outer flow 

channel. The inner passage consists of a circular tube with inner and 
outer diameters of 3>72 cm and 3*8 cm» respectively. A swirl vane assem** 
bly containing 12 equally spaced vanes with a swirl vane angle of 68.5 
degrees is iocated 11.9 cm upstream of the discharge end of the tube. 

The outer flow channel is a radial inflow passage with 24 adjustable 
swirt vanes near its inlet. The outer flow turns from the radial to the 
axial direction whereupon it exits in an annular fashion into a l4.6 cm 
diameter pipe where the mixing of the twu streams occurs. 

6 . 5.2 Counterswi ri i ng Case 

The inner and outer flows have swirl numbers of 0.49 and -0.51 • re- 
spectively. A nonuniform fine mesh of size 30 x 24 (Nl • 30, NJ ■ 24) 
is arranged to cover a pipe region of length 0.17 m and radius 0.0727 m. 
Inlet profiles of axial and swirl velocities are taken from the measured 

values at the location x * 0.002 m (x/D ■ 0.01375). as indicated in Fig- 

3/2 

ure 52 . Inlet k and c (« k*^ /t) are specified via 

k ■ 0.03 u^ (6.2a) 

I • 0.005 u/2 (6.2b) 

where 0 • 0.1454 m. The exit flow is paraliel, with zero axial gradients 
for ali variabies and v ■ 0. 

After a preiiminary calculation of 22 iterations and a further i89 
iterations, finai convergence is obtained. The criterion is that the 
maximum normaiized residuai source sum for each of the variabies mass, u, 
V, w, and p is less than 0.004. Radial profile computations of axiai 
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and swir) velocities u and w are shown in Figure 53 at the axiai station 

x/0 ■ 0.i37 (x • 0^02 m), and in Figure 56 at the station x/D ■ 0.276 

(x • 0.06 III). These figures show velocities normalized with respect to 

a fixed quantity u^» taken as 21.6 m/s as an estimate of the outer flow 

spatial average axial velocity* and the experimental data* quoted by 

Srinivasan and Mongia (69)* have been appropriately transferred to these 

• * 

figures for comparison. The two axial stations chosen lie in a central 
zone and show dramatically the differences. 

Notice that* whereas the agreement between swirl velocity profiles 
is very satisfactory at both locations, the agreement between axial velo- 
city profiles is not. In the present work the central zone is predicted* 
but it is not as large or with such large negative velocities as found in 
the experiment (29). Finally* Figure 55 exhibits the centerline axial 
velocity development as a function of downstream location. The velocity 
predictions at the nearest gridline to the axis are taken as representa- 
tive of axis values (r/D ■ 0) since 3u/3r is zero there. The central re- 
circulation zone is rather well predicted in the present study, although 
the axial velocity magnitudes are all underpredicted. 

6.5.3 Coswirling Case 

The flowfield for this experiment concerns inner and outer flow 
swirl numbers of 0.58 and 0.56, respectively. Inlet profiles of axial u 
and swirl w velocities taken from the measured values at the location x 
0.002 m (x/0 ■ 0.01375), as indicated in Figure 56. Other boundary 
conditions are as taken previously. Radial profile confutations of axial 
and swirl velocities are shown in Figure 57 at the axial station x/0 * 
0.0688 (x • 0.01 m) * and in Figure 58 at the axial station x/0 ■ 0.2063 
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(x 0.03 m). These may be compared directly with experimental evidence 
given in the corresponding figures. Velocities are normalized as in the 
previous counterswirl case. The present predictions are very encourag* 
ing when compared to the experimental data. Finally, Figure 59 illus* 
trates the centerline axial velocity development as a function of down- 
stream distance. Again, velocity predictions are taken from the nearest 
gridline to the axis. Notice that predictions of the present study com- 
pare quite well with the experimental data, especially for x/D < 0.5. 




CHAPTER VII 
CLOSURE 

7.1 Conclusions 

The present research is concerned with a specific sub'problem of 
the complete three-dimensional problem and is restricted to axisymmetric 
geometries under low speed and nonreacting conditions. This is an area 
in which there is a need for more fundamental research, particularly re- 
garding recirculation zones where most of the burning occurs in reacting 
flows. 

Many factors affect the existence, size, and shape of both the cor- 
ner recirculation zone and the central toroidal recirculation zone. A 
major outcome of the current study is the experimental and comoutational 
characterization of these for six basic flowfields. Parameter variations 
which define these flow conditions are side-wall angle a * 90 and 1>5 de- 
grees with swirl vane angle ^ ■ 0, AS, and 70 degrees. The size and 
shape of the recirculation bubbles for each flowfield <s illustrated as 
an artistic impression deduced from a collection of flow visualization 
photographs of tufts, smoke, or neutral ly-buoyant soap bubbles responding 
to the flow. 

Photographed results show that increasing ^ from 0 to kS degrees 
produces a shortened corner region and the appearance of a central bubble 
typically extending downstream to approximately x/0 ■ 1.7. Also, a pro- 
cessing vortex core is observed which stretches from the end of the 
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centra] region to the test section exit. A further Increase in ^ to 70 
degrees enlarges the central zone and vortex core with negligible effect 
on the corner region in those flowfields where It occurs. 

The effect of a on the nonswirling flows is negligible; however, a 
decrease from 90 to kS degrees apparently eliminates the corner bubble 
in each swirling flow. This decrease in a also causes the inlet flow to 
impinge more severely on the top wall, where larger axial velocities 
occur. 

A more detailed experiment consists of the measurement of time-mean 
veiocity components in the axial, radial, and azimuthal directions using 
a five-hole pitot probe. These measurements generally agree with the 
flow visualization results and provide a more complete understanding of 
each flowfield. Also, they constitute a data base allowing validity 
testing of flowfield predictions. Further, inlet velocity measurements 
reveal that the estimated inlet velocity profiles employed in the present 
swirling flowfield predictions are unrealistic. This arises from blade 
inefficiency, the presence of a hub, and the fact that the swirl vane 
exit station in the test facility is actually located 4 cm upstream of 
the test section inlet. This allows the central recirculation zone to 
begin upstream of the inlet where x/0 ■ 0. 

The flowfield prediction computer program is employed in a paramet- 
ric study to computationally determine trends of interest to combustor 
designers. The predicted effects of a and ^ on swirling flows are almost 
the same as those determined experimentally, even though the inlet veio- 
city estimates employed consist of flat profiles with 100 percent effi- 
cient swirl vanes located at the inlet with no recirculation. Further 
predictions using measured inlet velocities as boundary conditions were 


70 


obtained as a special test case to allow appropriate comparison with inea> 
surements for' the a ■ 90 and ^ ■ 45 degree flowfield. Poor agreement re- 
sults at x/D “ 0.5* although it is excellent for x/0 2 1.5< This dis- 
crepancy may be due to the lack of realism inherent in the k-e turbulence 
model (extended to include swirl shear stresses), inasmuch as the empiri- 
cal constants pertain to nonswirling flows. Another possibility may be 
an inaccurate specification of velocity and turbulence distributions at 
the inlet to the test section. 

7.2 Recommendations for Further Work 

Fundamental research should be continued in several areas. First, 
time-mean velocity measurements should be obtained at additional axial 
stations, especially in the upstream region, to provide more details and 
allow more extensive prediction validity tests. Second, movies of the 
flow visualization experiments should be obtained and analyzed to charac- 
terize the flow instabilities and temporal behavior of the recirculation 
zones and the trailing vortex core. Third, further details of the flow- 
fields are needed to allow a more complete fundamental understanding. 

This should consist of turbulence quantity measurements, including all 
of the Reynolds stresses. This would allow the deduction of more sophis- 
ticated turbulence models relating turbulent shear stresses with time- 
mean velocity gradients. 
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TABLE I 

SOURCE TERM IN THE GENERAL EQUATIONS 
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<(• 


- 

8x 


i£ pw 2^ ^ V 

3r r 2 

r 


. ^ A (,,) * s" 

r 


G - CpPC 


(C,cG - C^pc )/k 

















p. w ETC. 
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V 

\ 




O 



♦ s 



CONTROL VOLUMES C. U, V FACE 

AREAS An. A3, Aq AND A^^ FOR C. SIMILAR FOR U AND V 


The Three Control Voluries Associated With 
Points of the Three Grids 
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View of trie Swirl Vane Assembly 
From Upstream 






Figure 8. 


V I ev*. of the Sift/irl Vane Assennbly 
Frur Downstream With the In- 
strument ^or Measuring the 
Vane Annie 
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Figure 17. Velocity Coefficient Calibration Characteristic for the Five-Hole Pitot 
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Flow Vi sud I i zat ion Photographs of S'^^^ke-Wire 
Streaklines With Wall Expansion Angie a = 
90° for Swirl Vane Angles: (a) x/D = 0.5 

and (b) >'/D = 1.0 
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OF PObR QUALITY. 


Flow Visualization Photographs of Tu^ts 
in the r6-Plane Wit. Wall Expansion 
.^ngle a = 90® and Swi rl Vane Angie 1> 

= 45® for Axial Stations: (a) x/D = 

1 . 5 and (b) x/D = 2.5 
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Flow Visualization Fhotogrophs of Smoke 
Wire Streaklir.es With Wall Exoansicr 
Angle i = ^5° for Swirl Vene Angles: 
(a) C ^ 0°. fb) i- = AA-, and (c) c = 


Fiaure 27 
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Flow Visualization Photographs of Pathlines 
Indicated by Illuminated Neutral ly-Buoyant 
Soap Bubbles for Wall Expansion Angle^.i = 
90° and Swirl Vane Angles: (a) ; = 0° , 

(b) $ = 45' , and (c) 6 = 70° 
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Flow Visualisation Photographs of Path 
lines Indicated by Illuminated 
Neutral ly-Buoyant Soap Bubbles for 
Wall Expansion Angle n » ^*5° and 
Swirl Vane Angles: (a) b = 0", (b) 

? ^5’, and (c) b “ 70“ •> 
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DO 3.00 


Axial Position x/D 

Figure 30. Predicted Streamline Plots With 
Wall Expansion Angle a * 90° 
for Swirl Vane Angles: (a) 4- 

» 0°. (b) 4 = ^5°, and (c) 4 


2.0 

) « » 70 * 

Axial Posi tion x/0 

Figure 31- Predicted Streamline Plots With 
Wi>ll Expansion Angle a » ^5° 
for Swirl Vane Angles: (a) 

■ 0°. (b) (f - 45*. and (c) $ 
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(a) Corner Recirculation Zone 


Figure 33 Experimental Effect of Wall Ex- 
pansion Angle a and Swirl 
Vane Angle ^ on: (a) Corner 

Recirculation Zone Length and 
(b) Central Toroidal Recircu- 
lation Zone Length 















Figure 37- Measured Velocity 
and Swi rl Vane 



gure 38. Measured Velocity Profiles 






gure >9. Measured Veloctty Profiles 
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-Q2 0 Q2 0.4 0.6 0.8 | 

(a)u/Ug I 

Figure h 7* Probe Interference Effect on Measured 
Velocity Profiles With Wall Expan- 
sion Angle a « 90° and Swirl Vane 
Angle (J “ at Axial Station x/D 
= 1.0 for (a) Axial Velocity and 

(b) Swirl Velocity j 
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0.30 


(a) u/u^ 

Figure ^8. Probe Interference Effect on Measured 
Velocity Profiles With Wall Expan- 
sion Angle a • 90** and Swirl Vane 
Angle ^ « Ii5“ at Axial Station x/0 
• 2.5 for (a) Axial Velocity and 
(b) Swirl Velocity 








Fi 9 ure k$. Effect of Inlet Reynolds Number on 
Measured Velocity Profiles With 
Wall Expansion Angle a ■ 90" and 
Swirl Vane Angle ♦ • ii5* at Axial 
Station x/0 ■ 2.5 for (a) Axial 
Velocity and (b) Swirl Velocity 
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(a) u/u 


Figure SO. treasured Velocity Profiles From an Entire 
Uiameter Traverse With Wall Expansion 
Angle a ■ 90* and Swirl Vane Angle c> “ 
l>5* at Axial Station x/0 “ 1.0 for (a) 
Axial Velocity and (L) Swirl Velocity 
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(a)u/u„ 

Kigure 52. Inlet Velocity Profiles at x/D - 
0.01375 (x - 0.002 m) for the 
Counterswirling Coannular Pipe 
Flow Test Case (a) Axial Velo- 
city and (b) Swirl Velocity 



txx Reference (29), 


Experiment 

X 


X 


X 


y X ^ 
X ^ 

X 

X 

X 

X X 

X ^ 

» 1 L 

X 

1 L_J L- 











-Q5 0 0.5 1.0 1.5 
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Figure 5^. Velocity Profiles at x/D ■ 0.27^ (x ■ 4 cm) 
for the Counterswirling Coannul ar Pipe 
Flow Test Case (a) Axial Velocity and (b) 
Swirl Velocity 
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Figure 56. (Continued) 



















APPENDIX C 

DATA REDUCTION COHPUTER PROGRAM LISTING FOR 
FIVE-HOLE PITOT PROBE MEASUREMENTS 


159 


160 


C A COHPUTER PROORAH FOR DATA REDUCTION OF FTVE-HOLC PITOT 
C NEASURENENTS IN TURDULENTt SUIRLINOp RECIRCUL»>TIN(S » FLOU 

C IN COHDU8TOR OEOHETRXES 

C 

C UER8X0N OF NAY» 1981 

C 

C D L RHODE 

C NECHANXCAL AND AEROSPACE ENGINEER XNO 

C OKLAHONA STATE UNIVERSITY 

C STILLUATERp OK 74078 

C 
C 

CHAPTER 00000000 PRELIHINARIES 00000000 
C 

D I HEN810N HEDN < 9 > » HEDUNN < 9 > t HEDNN8 ( 9 > t HEDCNU ( 9 > » 

»HEDU( 9 ) tHEDVC 9 > t HEDU( 9 ) t HEDVT ( 9 ) pHEDUST < 9 > t 
«HEDVST(9> tHEDySTC?) fHEDVTS(9) tHEDDEL(9) »HEDDET(9) p 
♦HEDHNF ( 9 ) p HEDNX V ( 9 ) p HEDHIP ( 9 ) p HEDNPP ( 9 ) p HEDAN ( 9 > 

C 

CONNON 

• / C AL I D/CP I TCH ( 26 ) p C DEL T A ( 2 6 > p C VELCF < 26 > 

•/HE A8UR/RDETA < 48 p 24 ) p RPNHPS ( 48 p 24 > p RPCHPU (4ep24)«NDATA<48)p HAX JPT p 


• RDNPR8(48) 

• /OEON/X ( 48 > P R ( 24 ) P XND ( 48 ) P RND ( 24 > p DYPS 04 ) p 0 YNP ( 24 ) p 

• SNS ( 24 ) p N8T ATN p X I NCHS ( 48 ) p R I NCHS < 24 ) 
•/CALC/VT0TAL(48p24>>**^ *8p24> p V(48p24> pU<48p 24) t 

• VTSTAR(48p24)p. «R( 48p 24 ) t VSTAR( 48 p 24 ) p NSTAR( 48p 24 > . 

• P ZCHCF ( 48 P 24 ) P ^ rLCF ( 48 P 24 ) P DEL T A < 48 p 24 > f PET A < 48 p 2 4 ) p 
4 ANGHOH < 48 ) p UHEAN < 48 ) p HASS < 48 ) p HASFLO ( 48 > p U 1 N < 48 > 

REAL HASS p HASFLO 
LOGICAL IURITEpDIAGNS 

C SET IWRITE-pTRUE, FOR URITINO SOLN. ON DISK STORAGE 

C SET OIAONS-.TRUE, TO ACTIVATE DIAGNOSTIC WRITE STATEHENT3 

C N3TATN IS NO, OF AXIAL STATIONS FOR WHICH DATA IS 

C SUPPLIED 

C- HAXJPT IS HAXp no* OF RADIAL POSITIONS FOR ALL AXIAL 

C STATIONS 

C CPITCH IS CALIBRATION PITCH COEFF, 

C CDELTA IS CALIBRATION PITCH ANOLE(DEO.) 

C CVELCF IS CALIBRATION VELOCITY COEFF, 

C NCAL is NO, OF CALIBRATION POINTS 

C XINCHS IS AXIAL POSITION! IN, ) OF EACH PARTICULAR 

r RADIAL TRAVERSE 

C NOATA IS NO* OF RADIAL LOCATIONS FOR WHICH DATA IS 

C supplied for each PARTICULAR RADIAL TRAVERSE 

C RDNPRS IS HEASURED OYNAHIC PRESbURE ENTERING THE 

C SWIRLER 

C RINCHS IS RADIAL POSl TION( IN . ) OF EACH HEASUKEHENT 

C LOCATION IN EACH PARTICULAR RADIAL TRAVERSE 

c rbeta is neasured flow anol£(DEG,) where yaw 

C ANGLE BETA « .160,0 ’ RBETA 

C RDNPRS IS HE ASURED VOLTS FOR PNORTH - PSOUTM 

C DIFF, PRESSURE 

C RPCHPW IS HEASURED VOLTS FOR PCENTER • PWEST 

C DIFF, PRESSURE 

C PICMCF IS REDUCED VALUE FOR PITCH COEFF. 

C DELTA IS REDUCED VALUE FOR PITCH AN0L£(DE0,) 

C- VELCF IS REDUCED VALUE FOR VELOCITY COEFF. 

C seta is REDUCED VALUE FOR PROBE YAW ANOLE(DEO.> 

C vtotal is total vector velocity haonitude<m/sec> 

C U IS AXIAL velocity <H/SEC) 

c radial velocity (H/SEC) 

c u IS SWIRL velocity (H/SEC) 

C --VTSTAR IS DIHENStONLESS TOTAL VELOCITY 

C USTAR IS DIHENSTONLESS AXIAL VELOCITY 

r VSTAR IS DIHCNSI0NLES8 RADIAL VELOCITY 

C-— W8TAR IS OIHENSIONLESS SWIRL VELOCITY 

C ALL PRIHARY USER INPUTS ARE LOCATED HL*.E 


tWRITC«,TRUC* 

DIA0NS»,TRUE« 

ALPHA-90, 

PH I -45 • 


16 ) 


VISC0S«1.8E-5 

NSTATN-A 

HAXJPT-20 

PATH-73,8 

TATH-28. 

ITM8 

JT-24 

RL AROE- U , 75/ < 2 * 0*39 * 37 ) 

R8HALL-RLARQE/2,0 

READ<Sf205) HEDHtHEDUHNf HEDU*HEDM«HEDyf 
• HCDVT » HCDUST » HEDMST t HEDU8T f HEDMTS « HEDDEL t HEDBE T » 

•HEDNH8 » HEDCHU » HEDHHP » HEOH I » HEDN I P » HEDHPP t HEDAH 
205 F0RHAT<9A4> 

C INITIALIZE VARIABLES TO ZERO 

CALL INIT 

c read five-hole pitot calibration data 

NCAL-25 

DO 10 I-1»NCAL 

READ<5»210 ) CPXTCH( I ) # CDELTAC I > »CVELCF< I ) 
to CONTINUE 
210 F0RHAT<3F10,5) 

IFvDIAQNSJ «RITE<6f400) (CP^*^CH< I ) f 1*1 1 25> 

IF(0IA0N8) WRITE(4»400) it 0*:\ ^M^'Dtl-l.^S) 
ir(DIAONS) URITE(4f400) (CVcLCh < I ) » I-l f 25) 

400 F0RHAT(///flX»13<Fe«4»lX)f//fSXf 12<F8,4>) 

C READ RAW HEASURED DATA TO BE REDUCED 

DO 30 I-lfNSTATN 

READ<5f230) XINCHS( I ) t NDATA( I > t RDNPRS( I ) 
JPTS-NDATA<I> 


20 

30 


32 

33 
220 
230 


37 

450 

500 


DO 20 J-lfJPTS 

RE AO < 5 » 220 ) R I NCHS U ) t RBET A ( 1 1 J > f RPNHPS < I » J > « RPCHPUl < 1 1 J ) 

CONTINUE 

CONTINUE 

CONVERT X AND R TROH INCHES TO HETERS 

DO 35 I^ltNSTATN 

X<l>«XINCHS(I)«0,02S4 

JPTS-NDATA(I) 

R?J?^Rl/NiH^U?*0.0354 * 

CONTINUE 
CONTINUE 
FORMAT (4F 10. S> 

FORMAT! IFlO.StlllOtlFlO.S) 

IF<DIAONS> URITE(Af470> (NOATA< I > > I«ItNSTATN> 

<X<I>fI-lfNSTATN) 

(R(J)tJ-liJPTS) 


IF<DIAONS) URITE(6>4S0) 
IF<PIAGNS> URITE(AtSOO) 

DO 37 I-t»NSTATN 
IF<DIAONS> URITE(6fS00) 
IF<DIAONS) URITE(6>S00) 
IF<DIA0N8> URITE<6fS00) 
CONTINUE 

F0RMAT<///»4v>>1(F8.4»1X) ) 
F0RMAT<///.20X.I0<F8.4)> 


<RBETA(ItJ)tJ>lf JPTS) 

< RPNMRS < I . J ) t J- 1 . JPTS ) 
<RPCMPW<It J).J-l»JPTS» 


C 

CHAPTER 1 1 1 I I DATA REDUCTION I I I 1 1 1 

C 


470 FORMAT!///. 40X«1(I8»1X)) 


c calc pichcf and interpolate for delta from 

C PITOT CALIBRATION CURVE 


IDID>0 

DO SO I>1.N8TATN 
JPTS-NOATA!!) 

DO 40 J-I.JPTS 

IFv !RPCMPW!I. J> .EO. 0.0) .AND. !RPNMPS!I»J) .EQ. 0.0»' CO TO 30 


P ICHCF < I f J )«RPNMPS< I . J > / < RPCMPH! 1 1 J ) f 1 . E -4 > 

IF! iPICHCF!l.J).0T.3.399) .OR. (PICHCF! I . J) .LT. -3. 759> ) 00 TO 38 
IFUDID .EQ. 0) DCLTA!I.J)>SPLINE!CPITCH. 
f CDELTA.NCAL.PICHCFil.J)) 

IF!IOID .OT. 0) DELTA! l.J>-SP<CPITCH.CDELTA. 

• NCAL. PICHCF! I tJ>> 

IDID-1 
GO TO 40 
38 CONTINUE 


n n 


ORNSHIAL PAGE IS 
OP POOR QUALITY 


DELTA(ItJ>-0.0 
UR1TE(A»850> ItJ 

8S0 E0RHAT(20Xr'PlCHCF IS OUT OF RANGE OF CALIBRATION AT 
•'>I3>' AND J>'tI3> 

40 CONTINUE 
50 CONTINUE 

C INTERPOLATE FOR UELCF FROM PITOT CALIBRATION DATA 

IDID>0 

DO 80 Ial>NSTATN 
JPTS-NDATA(I) 

DO 70 J>1»JPT8 

lF((RPCMPU<ItJ) .EQ. 0.0) .AND. <RPNHPS(IfJ) .EQ. 0.0) > 80 TO AS 
IF((AB8(DELTA(I»J))> .GT. S8.0) GO TO AS 
IFCIDID .EQ. 0> UELCF<I.J)-SPLINE(COELTAf 
♦CVELCF , NCAL t DELTA ( I » J ) > 

IFdDID .QT. 0> UELCF(I>J)«8P(C0ELTAtCUELCF* 

*NCAL>DELTA(I»J>) 

IDlO-1 
60 TO 70 
AS CONTINUE 

UELCF(I>J)>0.0 
URITE<Af890> ItJ 

890 FORMAT <20X> 'DELTA IS OUT OF RANGE OF CALIBRATION DA^A 
♦AT l-'fI3»' AND J»'.I3) 

70 CONTINUE 
80 CONTINUE 

DO as I*1»NSTATN 

IF<DIAONS> URITE<AfSOO> (PICHCF<I>J)rJ>lf JPTS) 

IF(DIAGNS) WRITE(A»SOO) (DELTA( I . J> > J-lt JPTS) 

IF<DIAONS> URITE<A>SOO) <UELCF( It J> f J>lt JPTS) 

35 CONTINUE 

CALC MAGNITUDE OF TOTAL MEAN VELOCITY VECTOR AND 

U* V. » U COMPONENTS 

RHO-PATM* a 3A . 0/0 . 102 ) / < 287 . 0« < TATM4273 . 0 ) ) 

PI-3.14159 
DO 100 1-lrNSTATN 
JPTS-NDATAd) 

DO 90 J> If JPTS 
BETA( 1 f J)-3A0.-RBETA(If J) 

IF( (RPCHPUdf J).EO.O.O) .AND. (RPNMP8( I • J) .EQ.0 . 0) ) BETAdfJ)«0.0 
VTOT AL (I f J > -SORT ( ABS < 2 . 0/RH0» VELCF d . J ) 4RPCMPU dfj)«133.9)) 

Udf J)-VTOTAL<lf J) « COS(DELTA(If J)»PI/180.0) » 

♦ COS(BETA<If J)«PI/180.0) 

Vdf J)«VTOTALdf J) * SlN(DELTAdfJ)«PI/180.0) 

Udf J>>VTOTAL(I«J) « COS(DELTAdf J)«PI/180.0) * 

♦ SIN<BETA(If JXPl/180.0) 

90 CONTINUE 

100 CONTINUE 

IF(DIAGN8) URITE(Af500)(VT0TfL<I. J). J-lf JPTS) 

IF(DIAONS) URITE(Af500XU<If J.'f J-lf JPTS) 

IF(DIAGNS) URITE(Af500><V<If J)f J«lf JPTS) 

IFCDIAONS) WRITE<Af500)(U(If J)f J«lf JPTS) 

CHAPTER 222222 AUXILIARY CALCULATIONS 22222 
C 

DO 130 I-lfNSTATN 

C CALC GEOMETRIC QUANTITIES 

JPTS-NDATAd) 

JPTSM1-JPT9-1 

OYPS<1)-0.0 

DYNP ( JPTS ) -2 . 0« ( RLARGE-R < JPTS ) ) 

DO 110 J-lfJPTSMl 
OYNP< J)-R(J+1)-R< J) 

DYPS(J+1)-DYNP(J) 

110 CONTINUE 

DO 115 J-lfJPTS 

SNS( J)>0.5«(DYNP(J)'fDYPS<J)) 

115 CONTINUE 

IF(DIAONS) URITE<AfSOO) (DYNP< J) . J-l f JPTS) 

IF(DIAGNS) WRITE<6fS00) <SNS( J) t J-lt JPTS) 

FLOU-0.0 
UMOM-0.0 
DO 120 J-I.JPTS 
ARDEN-RHOBR ( J > »SNS ( J ) 


ORIQtNAL PAQE 18 
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UHOM-UMOH-t-ARDEN«U < I » J > «U ( I » J ) «R ( J > 
rLOU«FLOU'l’ARDeN»U (I t J > 

tF(DIAONS) URITE(6>900) J>FLOU>ARDEN»UHOMrRHO 
900 F0RHAT(///>IS»4<ri2.4.2X>) 

120 CONTINUE 

ANOHOM(I)-UHOH 

NAS8<I>-2.O»PI«PL0U 

UNEAN < I ) -HASS ( 1 ) / < RH09R X«RL AR0E»«2 ) 

130 CONTINUE 

IF(0IA6NS) URITE(4»4S0> <UHEAN<I>*I-1»NSTATN) 

IF<OIAQNS> URXTE(A»4S0) (HASS(I>»I-lrN8TATN) 

C N0N0IHEN8I0NALIZE VELOCITIES 

00 ISO I-1»NSTATN 
XN0<I>-X(I)/<2.0«RLAR6E) 

JPTS-NDATA< I ) 

UIN<I)-(SQRT(2.0/RH0«RDNRRS(I>«249.08))«(A.312/S.94>««2 
NASFL0< I )-2.0«PI»RH0«UIN( I >«R8HALL.«»2/2.0 
DO 140 J-lfJPTS 

VTSTAR< I f J)-VTOTAL( I* J) /UIN<I ) 

USTAR(I»J>-U(ItJ)/UIN<X) 

VSTAR<I>J)-V<IfJ)/UIN(I) 

U8TAR<I>J)-U(If J)/UIN(I) 

140 CONTINUE 
ISO CONTINUE 

DO 140 J-lfHAXJPT 
RNO<J)-R<J)/<2.0»RLAR6E) 

160 CONTINUE 
C 

CHAPTER 3 3 3 3 3 OUTPUT 3 3 3 3 3 3 3 

C 

IF(.NOT. lURITE) 00 TO 16S 
URITE<11) X 
URITEdl) R 
URITE<11> U 
URITEdl) V 
URITEdl) U 
16S CONTINUE 

URITE(6t3ll) 

URITE(6t32S) ALPHA 
URITE(6>330) PHI 
URITE(6r33S) RSHALL 
URITE(6t340) RLARGE 
URITE(6>3SS) VISCOS 
URITE(6»360) RHO 
C 

CALL WRITE < 1 1 1 > NSTATN >1 f IT > JT » X » R f MASFLO f HEDHHF ) 

CALL WRITE <1 » 1 > NSTATN »lfIT>JT>X»Rf HASS t HEOH ) 

CALL URITE<ld»N8TATNf lfITrJT>X>R>UINfHEDHIV) 

CALL URITEdfl»NSTATN«lrIT>JT>XtR> UHEAN t HEDUHN ) 

CALL WRITEd>l>N8TATN»lfIT>JT»X»R>AN0H0HrHEDAH) 

CALL PRINT <1 r 1 > NSTATN » HAX JPT rIT>JTrX>RtUf HEDU ) 

CALL PR I NT < 1 . 1 » NSTATN » HAX JPT »ITfJTfXfR»V. HEOV ) 

CALL PR I NT <1 f 1 • NSTATN r HAX JPT >ITfJT>XfRfUrHEDU) 

CALL PR INT (1 f 1 > NSTATN r HAX JPT > I T > JT » X » R r DELTA t HEODEL ) 

CALL PRINT d » 1 » NSTATN » HAX JPT » IT f JT » X »R * SETA»HEDBET ) 

CALL PRINT d f 1 » NSTATN # HAX JPT f I T » JT f X » R . VTOTAL »HEDVT ) 

CALL PR INT <1 « 1 r NSTATN » HAX JPT f I T f JT r XNO f RND « USTAR f HEOUST > 

CALL PR INT d f 1 > NSTATN t HAX JPT 1 1 T f JT > XND f RND t VSTAR f HEDVST ) 

CALL PRINTd f 1 fNSTATNf HAXJPTf IT» JTrXNDf RNDrUSTARf HEDUST) 

CALL PR INT d r 1 fNSTATN f MAX JPT f IT f JT r XNO • RND » VTSTAR . HEDVTS ) 

C ALL PR I NT d > 1 » NSTATN > MAX JPT • I T f JT f X INCHS t R I NCHS • RPNMPS , HEDNHS ) 
CALL PRINT < 1 f 1 f NSTATN t MAX JPT 1 1 T > JT t X INCHS > RINCHS • RPCMPW r HEDCHU ) 
CALL URITEdtlfNSTATNf If IT.JT.XINCHStRINCHSrRDNPRSfHEDMIP) 

CALL PRINT<lflfNSTATNfMAXJPTfITfJTfX INCHS t R INCHS # P I CHCF t HEDMPP ) 
STOP 

C FORMAT STATEMENTS 

311 F0RMAT<lHlfT37f 'AXISYHMETRICfISOTHERHALf GT COMBUSTOR FLOUFIELD 
•MEASUREMENTS 'f//fTS3f 'USING A FIVE-HOLE PITOT PROBE') 

32S FORMAT(////f T40f 'EXPANSION ANGLEIOEG.) « ' f T77f 1PE13. 3) 

330 FORMAT!//, T40f'SUIRL VANE ANGLEIPEG.) - ' , T77, IPE13. 3) 

33S FORMAT (//. T40, ' INLET RADIUS(M) -' ,T7r, 1PE13.3) 

3A0 F0RHAT(//,T40f 'COMBUSTOR RAOIUS(H) , T77, 1PE13.3) 

3S3 FORMAT <//,T40f 'LAMINAR VISCOSITY<KG/N/SEC) *' «T77, IPE13.3) 
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3A0 F0RMAT<//»T40# 'HENSITY<KG/CU. M> •' *T77f lPEi3. i*//// > 

END 

r • 

SUBROUTINE INIT 

C**************************************************************** 

c 

COHNON 

*/NEASUR/RBeTA<48f24) •RPNHPS<48r24) rRPCHPU<48t24) fNDATA<48) fNAXJPT* 

# RDNPRS ^ 48 ) 

♦/0£0H/X(48) fR<24) f XND<48) tRND<24> 8DYP3(24) t DYNP(24) f 

♦ SN8(24)tN6TATNfXXNCH8(48)fRINCH6(24) 
♦/CALC/VT0TAL<48f24)»U(48f24)#g(48»24)»M(48f24>f 

« gT8TAR<48»24>»U8TAR<48f24)fU8TAR(48f24>»U8TAR(4ef24)f 

♦ PICHCF<48»24)»VELCF<48t24)»DELTA<48»24)fBETA<48»24>9 

# ANOHON( 48 > »UHEAN<48 > t HASS ( 48) »HA8FL0( 40 ) tUIN < 48 ) 

REAL HASSfHASFLO 

C 

DO 20 I»lrNSTATN 
HA8FL0(I>»040 
MASS(I)*040 
AH6HOH(I)«040 
UHEAH(I)«0.0 
UINOoOaO 
DO 10 J-1»HAXJPT 
OTOTALdf J>-OaO 
U<If J)*0*0 
0<lf J)a0«0 
U(If J>«0.0 
VTSTARdf J)-0.0 
USTARdf J)>0«0 
MSTARdf J)«0«0 
USTARdf J)»040 
RBETAdf J)s0.0 
PETAdf J)«0.0 
PPNHPS(IpJ)- 040 
RPCMPU<If J)-0*0 
PICHCFdf J)»OaO 
VELCFdpJ>-0.0 
DELTAd»J)«040 
10 CONTINUE 
* 20 CONTINUE 
RETURN 
END 
C 

FUNCTION SPLINE(Xf FXp Nt XI) 

C CUBIC SPLINE CURUE FITTING IN 2 DIHEN8I0NAL DATA PLANE 
C INPUT VALUES I 

C Xf FX DATA ARRAYS t ONE DIHENSIONALf X IN INCREASING ORDER 

C N NUHBER OF DATA POINTS IN X» HAX 24 

C XI POINT OF INTERESTf WHERE F<X1) IS TO BE FOUND 

C 

C RETURN VALUE I 

C SPLINE OR SP • F<Xl) 

C THIS ROUTINE ACTIVATES ROUTINE ABUILD» Ht AND GAUSS. 

C FOR INTERPOLATION OF A LARGE NUHBER OF DATA POINTS p FUNCTION 

C SPLINE MAY BE CALLED ONLY ONCE # AND SUBSEQUENT CALLS MAY USE 
C ENTRY POINT 8P. 

DIHEN8I0N X(l)f FX<l)f A<24t27) 

C CONSTRUCT SPLINE hATRlX 

Ni«NiX 

DO 10 I*lf N 
DO 10 J«lt Hi 
10 Adfj)*0. 

Hl*N-l 

DO 20 I«2f HI 
CALL ABUlLDCXf FXt At Nt I) 

Aafl>«H<Xp2) 

A(lf2)«-H<Xf 1)-M<X*2) 

A<lf3)-H<X8l) 

H2-N-2 

A(NfH2}«H(X8Hl) 


20 
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A(Nfrtl)— H(XrM2)-H<XrHl) 

A<N»N)*H<XfM2> 

C FIND SECOND DERIVATIVES 

CALL 0AU88<Af Nt Nl) 

ENTRY SP(Xr FXr N» XI) 

C FIND F<X1> 

DO 40 Ul9 HI 

n-i+i 

IF<X1 •EQ« X<I>) 00 TO SO 

IF(X1 .LT. X(l) .ANO« XI *0T« XdD) 00 TO 41 
IF(X1 .OT. X(I) *AND* XI *LT. X(I1> ) 00 TO 41 

40 CONTINUE 

IF<X1 .EQ. X(N>> 60 TO 60 
UR1TE(6» 42) XI 

42 FORHAT<" Xl-'» 014. 7t " OUT OF INTERPOLATION RANOEr RETURNED VALUE 
*-0^ > 

SP«0. 

SPLINE «0. 

STOP 

41 CONTINUE 
11*1+1 
H1*H(X»I) 

HX«X<I1)-X1 

HX2«X1-X(I) 

FX1-HX«43/HI-HI*HX 

FX1*FX1«A(I#N1) 

ST0*MX2**3/HI - HI4HX2 
FXl-(FXl+8T0«A<IlrNl) )/6. 

SPL INE* ( FX ( I ) *HX+FX < 1 1 ) «HX2 ) /HI +FX 1 
SP*SPLINE 
RETURN 
C 

SO CONTINUE 

SPLINE*FX(I> 

SP»SPLINE 

RETURN 

C 

60 CONTINUE 

SPLINE*FX(N) 

SP-SPLINE 

RETURN 

END 

FUNCTION H<X>1) 

C*»««*.»*XX.«««*««***..*«******************.*****.****«.***.«*«********* 

C CALCULATE DELTA X UHICH IS USUALLY CALLED AS H. 

C««.**»*«*«««««««««*»K«)ll«)ll«*««««*«*«*««**«*«»*«X******»«)l>****««***««.**. 

DIMENSION X(l> 
ll-I+l 

H-X(I1)-X(I> 

RETURN 

END 

C 

SUBROUTINE ABUILD(X> Ft Af Nr I) 

C CONSTRUCT SPLINE MATRIX FOR FINDING 2ND OERIUATIVES. 

DIMENSION X(l>f F<l)r A<26»27) 

IMl-I-1 

Il-I+l 

N1-N4-1 

STO«H<X,I) 

HIMl«H<XtIMl> 

AdflMD-HIMI 

A( Irl )>2.«(HIHI.ST0) 

A(IrIl)>STO 

A<I»Nl)-< (F<I1)-F<I>)/ST0 - <F<I)-F<IMI))/HIM1 >*6. 

RETURN 

END 

C 

SUBROUTINE OAUSS<At Kr M) 

C OAUSS-JOROAN ELIMINATION 
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DIMENSION A<26f27> 

Ml-M-l 

K1«K-1 

DO 3 L»t» K1 
Ll-t+1 

DO 3 I«Llt K 
C0N8T«A<IfL)/A(LfL) 

DO 3 J-Lf M 

3 A<tf J)-A(I»J)-C0N8T«A(L»J) 

DO 6 1*1 f K1 
Il-Z-fl 

DO 6 L-Ilf Ml 
CONST»A(IfL)/A<LrL) 

DO 6 J«Ir M 

6 A(If J)*A<I»J>-C0N8T«A<Lf J) 

DO 10 I»l» K 
A(I»M)-A(I»M)/A<Irl) 

10 A(IfI)«l« 

RETUAN 

END 

C 

SUBROUTINE PRINT < 1ST ART t J8TART f NI t NJ . I T f JT r X t Y # PHI r MEAD > 

C 

DIMENSION PHKITf JT)#X<IT)fY(JT)fHEAD<9>fST0RE(48) 

ISKIP-1 

JSKIP*1 

URITE<6fllO)HEAD 

ISTA«ISTART-12 

100 CONTINUE 
ISTA-ISTA+12 
XEND«ISTAfll 
IF(NI.LT,IEND)IEND«NI 
URITE<6flll)<ItI>ISTA»IEND»ISKIP) 

URtTE<6»114)(X<I)f I»ISTAf IEN0*ISK1P> 

URITE(4fU2> 

DO 101 JJ«JSTARTtNJ»J8KIP 
J-JSTART4NJ-JJ 
DO 120 Z«ISTAfIEND 
A»PHKlf J> 

ir(AB8(A>*LT«l*E-20) A«0*0 
120 STORE(I)«A 

101 URITE<4f 113)J»Y( J>*(STORE<I)f I«ISTA»IEND»1SKIP) 
IF(IEND*LT«NI>00 TO 100 

RETURN 

110 FORMAT ( IHO » 1 7 ( 2H«- ) # 7X » 9A4 1 7X # 1 7 < 2H-» > ) 

111 F0RMAT(lH0rl3H I » fI2fllI9) 

112 F0RMAT(8H0 J Y> 

113 F0RMAT<13f0PF8*5flXflP12E9.2) 

114 FORMATdlH X « tFS.St 1 1F9 *5 ) 

END 

C 

SUBROUTINE WRITE ( ISTART t JSTART f NI r N J » IT f JT • X # Y f PHI r HEAD > 
C 

DIMENSION PHiaT)fX(IT)»Y< JT)#HEAD(9)fST0RE(48) 

ISKIP-1 

JSKIP-1 

URITE<6f 110>HEAD 
ISTA-ISTART-12 
100 CONTINUE 

ISTA-I3TA412 

IEND-ISTA411 

IF(NI .LT* lEND) lEND-NX 

URITE(6flllHI»I-ISTA»XENDf ISKIP) 

URITE(A»114>(X<I)d-ISTAtIENDf ISKIP) 

DO 101 JJ-J8TARTfNJtJ8KIP 
J-JSTART4NJ-JJ 
DO 120 I«ISTA»IEN0 
A-PHKI> 

XF(AB3iA).LT.l.E-20) A«0,0 
120 STORE(I)«A 


101 URlTe(6>113) (8T0Re(l>fI>I8TA»IEN0>t8KIP> 
IF<XEN0.LT.NI)00 TO 100 
RETURN 

110 f0RNAT(lH0»17(2H*-)«7X>9A4f7Xtl7<2H-«l) 

111 RORNATdHOflSH I - »I2tlll9) 

113 P0RNAT<12XflPl2E9.2> 

114 F0RMAT<11H X • *F8.S*11F9.S) 
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